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The Growth of Ice Crystals on Freshly Cleaved Covellite Surfaces 


By J. Hauietr 
Department of Meteorology, Imperial College, Londont 


[Received January 23, 1961] 


ABSTRACT 


Oriented ice crystals growing as thin plates on freshly cleaved surfaces of 
covellite, when examined by reflected light, show interference colours which 
give a measure of the crystal thickness. Some growing crystals remain con- 
stant in colour, failing to thicken in a direction normal to the basal plane. 
Growth can be initiated on the basal plane by contact with a thicker crystal, 
by increase of the excess vapour density, or by causing mechanical damage 
to the surface. 

When two crystals of different thickness touch, layers grow from the point 
of contact across the thinner crystal, with velocity inversely proportional 
to the thickness, for thickness 200 to 1000 &. This implies that the layers are 
growing primarily by surface diffusion of adsorbed molecules to their edges. 
This surface diffusion is characterized by a distance x, from which a growing 
layer collects molecules; measurements indicate that x has a maximum value 
of 7-5 4 at —6°c. The shape of layers and crystals are circular, polygonal, 
or dendritic depending on their thickness, and the excess vapour density of 
the environment. The relevance of the results to the variation of ice crystal 
habit with temperature is examined. 


§ 1, INTRODUCTION 


In a previous investigation of the epitaxial growth of ice crystals on single 
crystalline substrates (Bryant et al. 1959), it was discovered that crystals 
growing on freshly cleaved surfaces of covellite (natural cupric sulphide) 
were often very thin, and gave rise to brilliant interference colours when 
viewed by a microscope with vertical illumination. As each crystal grows 
vertically the colour changes, enabling its rate of growth to be found. In 
the present investigation this technique has been exploited to select 
erystals which do not grow perpendicular to the basal plane, and to investi- 
gate the relation between the rate of growth of discrete layers across the 
crystal, and the layer thickness. Previous observations had shown that 
layers growing on p-toluidene (Kowarski 1935) and cadmium (McNutt and 
Mehl 1955) crystals grown from the vapour, and on KH,PO, crystals grown 
from solution (Bunn and Emmett 1949) propagated with a velocity 
inversely related to their thickness. Other substances, however, forexample, 
erystals of sucrose growing from solution (Albon and Dunning 1959) had 
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layers whose growth velocity was independent of thickness. By observing: 
the growth rate, shape, thickness, and spacing of the layers on ice crystals 
we hoped to learn something about the growth mechanism which would 
throw some light on the changes of crystal habit with temperature exhibited 
by ice crystals growing from the vapour (Hallett and Mason 1958). 


§ 2, EXPERIMENTAL TECHNIQUE 


The principle of the apparatus is similar to that used in the previous. 
investigation. The substrate crystal rested on the flat end of a brass rod 
which dipped into a vacuum flask containing liquid nitrogen (fig. 17). 
Its temperature was controlled by a resistance heater soldered just below 
the surface. The environmental conditions were controlled by the 
temperature of an ice-lined outer brass chamber, also cooled by liquid 
nitrogen. Its temperature was controlled by a second heater, and could be 
kept constant to within 0-2°c. It was measured by a copper—constantan 
thermocouple in contact with the ice surface. An estimate of the sub- 
strate crystal temperature was obtained from a thermocouple 40s.w.g. 
copper—constantan, which lay alongside the crystal. This gave a value 
which was too low, since it had a higher thermal conductivity than the 
substrate crystal and both were heated by the warmer environment. A 
correction was made empirically by observing the temperature difference 
when ice crystals neither grew nor evaporated, which should have occurred 
when the crystal and environment temperatures were just equal. The 
maximum correction was about 4°c. This effect was partly caused by the 
proximity of the warm objective. In practice this limited the working 
distance to a minimum of 4mm. The vapour density excess or deficiency 
Ap, was calculated from the difference of the equilibrium vapour density 
over ice at the environment and substrate temperatures. An error of 
0-2°c in either temperature would give rise to an error of 10~* g em- in the 
excess vapour density near 0°o, and 10-°gcm-? near —20°c. The whole 
cold unit rested on the stage of the microscope and could be moved with the 
usual vernier adjustments, about $cm in each direction. 

The optical system employed a 16mm bloomed metallurgical objective. 
The object was illuminated by a tungsten lamp through a vertical illumina- 
tor, viewed by a x 20 eyepiece, and recorded simultaneously via a part 
silvered mirror on 35 or 16mm colour film, Superanscochrome, Type 
PN 1 226 (tungsten). 

In performing an experiment, the current in each heater was adjusted 
to give any desired value of crystal and environment temperature, and hence 
vapour density difference. A freshly cleaved specimen of covellite about. 
3x3x0-5mm? was selected, placed in position with a pair of tweezers, 
where it came into thermal equilibrium with its surroundings in a few 
seconds. The microscope was then moved into place for viewing and. 
photography. 


} Figures 2, 5, 8, 9, 10 and 15 are shown as plates. 
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Fig. 1 
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§ 3. MEASUREMENT oF Ick CrysTaL THICKNESS 


Events took place too quickly for accurate visual observation, and all 
measurements were made from frames of 16mm film, taken each second. 
A calibration slide was made of Newton’s rings produced by freezing a 
lens of known radius of curvature on a surface of covellite whose surface 
had been tested against an optically flat surface, using white (tungsten) 
light for illumination. Under these conditions, each order of interference 
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had a characteristic colour and with a little practice the actual thickness of 
any crystal could be estimated to an accuracy of 150 to 2004. As crystals 
became thicker, the interference colours were fainter, and could not be 
distinguished when the crystal thickness exceeded about 2 1, beyond which 
all crystals became a uniform pink colour. 


§ 4. EXPERIMENTAL OBSERVATIONS 


4.1 
The first crystals to appear on the covellite invariably grew along those 
steps in the substrate which exceeded about 0-1 in height (see fig. 10). 
All crystals assumed a constant orientation related only to the substrate 
crystal and not to the direction of the step, which varied continuously. 
The presence of steps less than 0-1 in height could be seen when crystals, 
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nucleated elsewhere, grew across them, and showed a change of hue. This. 
is clearly shown by the large crystal in the upper part of fig. 2, which has. 
grown across steps of height 800 to 300 4. The total thickness of this erystal 
was about 50004. Crystals subsequently grew on flat areas of the surface, 
particularly when the excess vapour density was large. 


Fig. 4 
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growing vertically. 


All crystals which originated at steps, and most others, grew uniformly 
in thickness until they became pink. At low values of excess vapour 
density, however, a few isolated crystals increased only in diameter, the 
colour and thickness remaining constant. The growth rate of a typical 
crystal is shown in fig. 3. Measurement made on about 50 such crystals 
showed that the crystal diameter increases linearly with time between 15 
and 100p. Crystals which thickened invariably grew more slowly in 
diameter. Figure 4 shows that the rate of approach of two plates of 
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constant, nearly equal thickness does not decrease as they become close 
together and finally join. Curve (b) shows that the rate of approach does 
decrease when the crystals are thickening in the vertical direction. : 

The rate of vertical growth depended also on temperature, and prism- 
type crystals often appeared below — 25°c and above — 8°C. 


4.2. Crystal Shape 


The shape of the ice crystals depended on the crystal thickness, changing 
from circular to polygonal as the thickness increased from 3000 to 5000 a 
This is shown by the lower crystal in fig. 2, which was originally circular, 


Fig. 6 
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and is beginning to take polygonal form as it thickens. At high values of 
excess vapour density, the crystal shape became very irregular, as in fig. 5, 
and finally took a roughly dendritic form with branches lying along (1120) 
directions. These changes are summarized in fig. 6. There was not any 
dependence of shape on crystal diameter from 15 to 100 . 


4.3. Ageing of the Covellite Surface 


If a piece of covellite which had not been freshly cleaved was used as the 
substrate crystal, the type of growth was quite different from that 
described above. As the ice crystal diameter increased beyond about 
20, it began to grow upwards from the surface, forming a hollow hexagonal 
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‘cup. The covellite behaved in this way after about 10min exposure 
to the (London) atmosphere. 


4.4. Growth of Ice Layers 

Ice iayers were occasionally seen spreading over flat ice crystals already 
growing on the substrate. These layers were often initiated as two crystals 
of different thickness grew laterally, and joined. If one crystal was much 
thicker than the other, say 2000 4 and 10, a succession of discrete layers of 
different thickness often spread out from the point of contact, and grew 
across the thinner crystal (fig. 5). The shape of these layers depended upon 
their thickness and the excess vapour density of the environment as shown 
in fig. 7. Thin layers, less than 300A, were invariably circular (fig. 5, 
crystal in upper left), becoming polygonal with increasing thickness. 


Fig. 7 
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Layers between 800 A and 2000 A spread first around the edge of the crystal, 
and then inwards towards the centre (fig. 8). This effect appears to be 
related to the increasing importance of diffusion of water vapour through the 
air as the step height exceeds the mean free path of air molecules, about 
0-ly. Layers of thickness greater than 20004 spread with a sloping 
front, and usually have a greater velocity of growth in the (1120) directions, 
giving in effect a three dimensional dendrite (see fig. 2, lower crystal, and. 
figs. 9,10). When the touching crystals were of almost the same thickness, 
the thickening appeared to take place uniformly over the whole crystal. 
~ As soon as the difference exceeded 800 A, however, the thickening took place 
first around the edge as described above. 
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Layers were also occasionally initiated remote from the edge of a crystal, 
possibly by the inclusion of a foreign particle as the crystal grew. In 
this case layers of thickness up to 20004 spread out in circular form. As 
soon as any layer greater than 800 A reached the crystal edge, however, it 
spread around the edge, as if it had been initiated there and then filled in 
toward the centre. 

Observations of the rate of growth of layers were made from time lapse 
films. For layers whose thickness lay between 200 and 1000 A, the velocity 


Fig. 11 
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was found to be inversely proportional to the layer thickness. A typical 
set of measurements are shown in fig. 11. The velocity of layer growth 
increased with the excess vapour density, until a value of 0-2 x 10-* g em-* 
was reached, when the growth rate became nearly constant. If a crystal 
were thickening uniformly however, its rate of growth was approximately 
proportional to the excess vapour density. 

It was found that the rate of growth of circular layers was constant for 
radii from 10 to at least 80. A typical example of the growth of such a 
layer is shown in fig. 12. As in the case of complete ice crystals, the 
growth rate is reduced when the layer is also increasing in thickness. The 
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resolution was not sufficient to make measurements on layers with radius. 
less than about 10 p. 


4.5. Influence of Temperature on Velocity of Ice Layers 
In view of the strong dependence of ice crystal habit on temperature 
(Hallett and Mason 1958), an investigation of the influence of temperature- 
on the rate of layer propagation was made. The initiation of layers was a 
very uncertain event, and many growths were examined before a layer: 


Fig. 12 
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suitable for measurement was obtained. Hence, to facilitate this investi- 
gation, use was made of the above relationship between layer thickness and. 
velocity. The velocity of any layers of thickness between 200 and 1000 4 
was measured, and the velocity of a layer of 250A thickness deduced. 
All measurements were made with an excess vapour pressure of 3 
0-25 +0-05x 10-*gem-8, It was very difficult to obtain observations of 
layers at temperatures greater than about —5°c, as crystals would not 
nucleate on the covellite surface. It was necessary to begin the experiment 


1082 J. Hallett on the 


with the substrate at about —10°c, and increase the temperature to the 
required value after nucleation had occurred. This took about one minute 
to reach a steady value, and usually any plates present had grown vertically 
in the high excess vapour density. It was not practicable to change the 
environment temperature because of the large thermal inertia of the outer 
walls. The results are shown in fig. 13, with the dashed part of the curve 
representing the less certain results at higher temperatures. 


Fig. 13 
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The initiation of layers described above usually resulted in the distance 
between layers of order 30. Occasionally, however, when a thin layer 
formed after a thick layer, it would spread with a corresponding greater 
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velocity and eventually catch up the first to give a thicker, more slowly 
moving layer. This event happened very rarely, and only two cases 
have been recorded on film. It was found that as the layers approached 
each other, at a definite distance apart, each began to move more slowly, 
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until the two merged (fig. 14) and the resulting thicker layer again moved 
with a constant velocity. The interaction distance was about 5y 
at—17-5°c, and about 10, at —10°c with excess vapour density 
0-2 x 10-*g em-*. ' 
4.6. Evaporation 

Evaporation studies of ice crystals on covellite were difficult because of - 
the inability to change the temperatures sufficiently quickly, and observa- 
tions could only be made while crystals were increasing in temperature. 
Evaporation first began at corners and edges and eventually the crystal 
took the shape of a spherical cap. Layers present when evaporation 
began, usually evaporated as discrete layers before the remainder of the 
crystal. Occasionally crystals evaporated by the initiation of layers at the 
edge. 
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4.7. Mechanical Damage of Crystals 


A piece of 36s.w.g. stainless steel wire was cemented to the base of the 
objective so that the tip was just in the field of view. This was then 
racked down to a suitable crystal which was not thickening, and a scratch 
made in the surface by moving the crystal beneath the wire. This caused 
the crystal to grow vertically over its whole surface, the line of the scratch 
filling in slowly. 


§ 5. Discussion. EVIDENCE FOR SURFACE MIGRATION OF WATER 
MOLECULES 


5.1. The Variation of the Velocity of Ice Layers with Thickness 


Consider a crystal layer growing only parallel to the basal plane both by 
direct deposition of molecules from the vapour, and also by the collection 
of mobile molecules from the surface layer within a distance x of the growing 
edge. 
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If the density of surface molecules is small, then x represents the mean 
free path ofa molecule on the surface. The velocity of advance ofa straight 
step, U, of constant height h, is given by 

hs 2 es 0 Ben ae 

pik 

where A =mass flux of vapour cm~*sec, p;,=density of ice. This result 
is in agreement with the experimental result that U=const/h for values. 
of h between 200 and 1000 A when 2ah. We can therefore conclude that 
surface diffusion is important for layer growth under these conditions. The 
condition that 2x h is shown to be valid by an independent determination. 
of the value of x (see next section). 


5.2. Rate of Growth of Layers 


Consider the rate of growth of a circular layer, radius r, of constant 
thickness, by direct deposition and by surface migration of molecules from 


a distance w on either side. We take h<z, neglecting growth by direct. 
deposition. 
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There are two cases: 


1. r<a. Ifm is the mass of the layer, 


dm dr 2 
an = 2mpihr a =An(r+x)*. when = const ao : (2) 
2 pn 
d ; 
= =4Anrx when es = constant. Ee ano) 


Experimentally it is found that dr/dt=constant for layers r>10p, 
accurate measurements not being possible for smaller radii. This sets an 
upper limit of 10 u to the migration distance. 

It may be pointed out that as the whole crystal is subject to the ambient 
supersaturation, those parts remote from the growing edge must be 
returning as many molecules to the vapour as leave it—and must have a 
vapour pressure larger than the value for the ice step at this particular 
temperature. 

The lack of variation of layer velocity beyond a critical value of the 
excess vapour pressure would imply that the surface layer had reached an 
equilibrium and become saturated. 

Direct evidence for the absolute value of x for the ice surface comes from 
the two observations of the interaction of approaching layers: 5p at 
—17-5°c and 10 at —12°c at Ap=0-2 x 10-° g cm-3, which gives values of 
« of 2-5 and 5p respectively. 

In the case of growth of crystals directly on the substrate, we suppose 
that molecules will be collected from a distance x on the ice and a distance 
Xe on the covellite. In the case r>z2z, this again leads to the result 
dr/dt=constant. This is observed for crystals whose radius is greater 
than about 7p. It is also found that two approaching crystals do not 
interact until they are separated by at least 2, so that we conclude that 
to< lp. 


5.3. The Influence of Temperature 


If we consider the growth rate of a layer of constant thickness and at 
constant flux, then eqn. (1) gives U=const. x. Hence it is possible to 
interpret the variation of growth rate with temperature (fig. 13) as a varia- 
tion of x with temperature. The two values of # obtained in § 5.2 above, 
enable a scale to be placed on this curve, to give a maximum value of x 
of about 7-5 at a temperature of —6+$°C. 

This temperature dependence differs significantly from the rapid increase 
with decrease of temperature suggested by Burton e¢ al. (1951) for growth 
directly from the vapour in the absence of air. In view of the importance 
of this result it would be desirable to obtain further information by direct 
measurements of x from layer interaction over a wide range of tempera- 
ture and also in the absence of air. 
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§ 6. Tur INFLUENCE OF THE SUBSTRATE ON CRYSTAL 
GROWTH 


If we intend to apply the results obtained in the previous section to 
crystals growing freely in the atmosphere or in a diffusion chamber, it will 
be necessary to investigate the effect, if any, of the substrate. 

As was shown in a previous paper (Bryant et al. 1959), the substrate 
reduces the supersaturation required for nucleation and determines the 
orientation of the crystal. In the case of freshly cleaved covellite, however, 
some crystals grow which do not thicken in a direction perpendicular to. 
the basal plane. We can interpret this as being caused by an absence of 
suitable sites for growth in this direction—a hypothesis which is sup- 
ported by the sudden growth of these crystals when large excess vapour 
densities are produced, suggesting that homogeneous nucleation of layers 
is taking place—and also by the initiation of growth (at low excess vapour 
density) when the crystal is subject to mechanical damage. It appears that 
exposure to the atmosphere leads either to a deposit of particulate matter 
or an adsorption of gaseous impurities on the covellite surface. This 
causes regions of strain in the overlying ice crystal, which ultimately lead 
to the formation of imperfections in the basal plane, and consequent growth 
in the prism direction. Crystals perfect in the basal plane have not been 
observed on old covellite surfaces, or on other substrates, and although 
crystals are occasionally sufficiently thin to show interference colours, they 
quickly grow in thickness and become colourless. 


§ 7. CRYSTALS AND LAYER SHAPE 


It is difficult to interpret the changes of crystal and layer shape shown in: 
figs. 6 and 7 without a knowledge of the structure of the surface layer. The 
transition from polygonal and circular growth to dendritic growth with 
increase of excess vapour density takes place much more completely for 
layers on an ice surface than for crystals growing on covellite, and suggests. 
that the surface layer on ice is quite different from that on covellite. This. 
is confirmed by the lack of interaction between approaching crystals when 
growing on a covellite surface. The growth of dendritic layers on ice 
resembles more closely that which has been observed on freshly cleaved. 
muscovite mica (fig. 15) and also on biotite mica (Kleber and Weis 1958), 
indicating that the migration distance of water molecules on mica is 
comparable with that on ice, a point which requires further experimental 
investigation. 

§ 8. Ick Crystan GRowTH 


These results show that ice crystals can grow on a covellite substrate 
without growth taking place perpendicular to the basal plane, irrespective: 
of the temperature, from — 4 to — 40°0, providing the vapour density excess. 
is not too large. This is in striking contrast to the habit variation with. 
temperature observed for crystals growing in the atmosphere and the 
diffusion chamber, and also for the majority of crystals growing on the. 
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covellite surface. In all these cases, sites for growth parallel to the c axis 
must be produced, possibly by stresses of a solid impurity built into the 
lattice, or by thermal stresses occurring during growth. 

In all cases, there are present sites for growth parallel to the basal plane. 
These may be imperfections in the ice itself, or the step between the ice and 
the covellite surface. The habit of a crystal containing growth sites in 
both prism and basal plane may therefore depend upon the relative rates of 
surface diffusion on the two faces, providing the growth site concentration is 
not large. The dependence of habit on a surface parameter of this kind 
is to be expected in view of the immediate change of habit.as the crystal 
temperature is changed during growth in the diffusion chamber—which. 
does not occur for the ‘ perfect’ crystals growing on the covellite surface. 

The present experimental results only give the variation with tempera- 
ture of x on the basal plane, so that a complete explanation of the complex 
habit temperature relation in ice must await an experimental determination 
of the corresponding parameter on the prism face. 
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Steep-sided Trigons on Diamonds 


By Emzren M. Wiixs 
Clarendon Laboratory, Oxford 


[Received February 6, 1961] 


ABSTRAOT 


The natural octahedron surfaces of diamond exhibit many small triangular- 
shaped cavities or trigons. The majority are low-angle pyramids, but there is 
another vartety which has steep sides and a flat base. Contour measurements 
on the pyramidal form have already been reported; similar measurements 
have now been made on the steep-sided variety. It is found that the sides of 
these trigons are octahedron planes. Thus the steep-sided trigons must arise 
by some quite different mechanism from the etch process responsible for the 
low-angle pyramids. 


§ 1. INTRODUCTION 


TRIANGULAR-SHAPED cavities, commonly called trigons, have long been 
observed on the natural octahedron faces of diamonds. They are oppo- 
sitely oriented to the octahedron face, and the majority have sides which 
are inclined to the surrounding octahedron surface at angles of a few degrees. 
Not all of these low-angle trigons are simple pyramids, some are truncated 
with flat bases, and some have slightly concave sides (Tolansky 1955). 
The morphology of trigons is of interest in connection with theories of 
crystal growth and dissolution, and we have described the contours ot these 
low-angle trigons in a previous paper (Frank ef al. 1958). In addition, 
there is another type of trigon, found on natural octahedron surfaces, 
which has steep sides and a flat base; it is with this type that we are here 
concerned. Examples of the two forms are shown on fig. 17, trigon X is a 
low-angle pyramid, while trigon Y has steep sides and a flat base. The 
latter type is not as common as the pyramidal form, although it has been 
observed on all the diamonds we have studied. For example, on one typical 
stone there were more than 100 pyramidal trigons deeper than 100 4, but 
only eight steep-sided ones. (If trigons have depths less than about 100 4, 
it becomes difficult to make an accurate estimate of their shape.) 

The differences in outline between the two types of trigon may be seen 
from multiple-beam interferograms. In fig. 2 the fringes cross a steep- 
sided trigon ; there is an abrupt change of depth at the edge. Figure 3 shows 
fringes crossing a low-angle trigon; the sides slope down to the apex at the 


} Figures 1, 2, 3 and 5 are shown as plates, 
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centre. As mentioned above, contour measurements have already been 
made on the pyramids; we have now made similar measurements on the 
steep-sided trigons, and this paper describes results obtained on one of the 
stones previously used in studying the pyramidal trigons; a good quality 
‘blue white’ regular octahedron of about | carat. 


§ 2, EXPERIMENTAL METHOD 


To determine the inclinations to the octahedron surface of the sides of: 
a steep-sided trigon, two measurements are required: the depth ¢ and the- 
lateral dimension a as illustrated in fig. 4. 


Fig. 4 


Plan and side views of a steep-sided trigon showing the dimensions a and f 
required to determine the inclination of the sides to the octahedron surface. 


It was not possible to determine these dimensions on all the available. 
trigons owing to several limitations imposed by the method of measure- 
ment. The dimension a could be found from photomicrographs at a 
magnification of 3000 to within +0-3 micron, and by taking some 20: 
readings on each photomicrograph we obtain a value of a accurate to + 0-1 
micron. ‘This sets the lower limit to our measurement and we have only 
investigated trigons with values of a greater than 1-2micron. Additional. 
restrictions were also introduced by the interference method of determining 
the depth. To determine the correspondence between orders, the trigons 
had to have a large-enough surface area so that at least two fringes could 
be observed within the trigon, and also the method becomes difficult for- 
depths greater than about 7microns.. (For trigons deeper than about 
15 microns ordinary microscopy was used.) Again many trigons are located. 
near to the edges of the octahedron surface, and some of these were obscured. 
by the opposite octahedron surface when viewed by transmitted light in 
order to determine the correspondence between orders. As a result of these- 
restrictions the number of available trigons on any given face was rather- 
limited. Thus on the face shown in fig. 5, there are about 18 trigons, but. 
measurements could only be made on six of them: and on the other face: 
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studied, measurements could be made on three out of a total of 12 trigons.. 
The results for these nine trigons are given in the table below. 

§ 3. RESULTS 


a 


Face | Trigon number | Depth (microns) | Side length (microns) Pel selone) 


(degrees) 
A 1 4-5 52 75 
A 2 5-4 62 70 
A 3 17-5 109 75 
A 4 6-7 41 68 
A 5 5-9 64 70 
A 6 25-0 114 68 
B 7 57 100 70 
B 8 5-4 81 74 
B 9 5:7 90 71 


The accuracy of these results is estimated to be about 6°. As the angle 
between two intersecting octahedron planes is 70°39’; the results show 
that to within the accuracy of measurement the steep-sided trigons have 
sides which correspond to octahedron planes. These results are in marked 
contrast to those found for the pyramidal trigons where the sides are inclined 
to the octahedron surface by angles no greater than about 2°. It was also 
apparent from the Fizeau fringes that the base of a steep-sided trigon is. 
parallel to the octahedron surface. We have already noted that our method 
of measurement severely restricts the size and depth of trigons available 
for inspection. Therefore no significance is to be attached to the narrow 
range of depths in our results. There is a suggestion that the deeper steep- 
sided trigons have also larger surface areas but this again is probably not 
significant. 


§ 4, DIscUSSION 


The present results together with those previously obtained (Frank et al.) 
show that the trigons observed on diamond consist of two distinctly 
different types: the pyramidal form with shallow non-crystallographic 
sides and a flat-based form whose sides are octahedron planes. This 
distinction is relevant to recent discussions as to whether natural trigons 
are the result of etch or growth processes. Tolansky (1955) argued that all 
natural trigons are due to growth because trigons produced in etching 
experiments (Fersmann and Goldschmidt 1911, Wilhams 1932, Omar ef al. 
1954) are oppositely oriented to the natural ones. In addition he pointed 
out that edges of etched trigons in the octahedron surfaces are generally 
rounded in contrast to the straight sides of natural trigons (although trigons 
with straight edges have since been obtained by etching (Omar and Kenawi 
1957 and Tolansky and Patel 1957.) Other support for a growth mechanism 
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was deduced by Tolansky and Wilcock (1947) from a study of a flat-based 
trigon, the base of which lay very exactly in the same plane as did the 
surface of the stone some distance away. However, we have recently 
presented arguments, based on surface features associated with pyramidal 
trigons in contact, which imply that this type of trigon is due to etch (Frank 
et al. 1958). Subsequently, pyramidal trigons oriented in the same way as 
the natural ones have been obtained by etching (Frank and Puttick 
1958). Thus, although it seems that pyramidal trigons with non-crystallo- 
graphic sides are due to dissolution, the flat-based trigons with octahedron 
sides must arise by a different mechanism. 
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ABSTRACT 


Electrical resistance, at temperatures from —196 to 700°c, was used as an 
indicator of crystallographic order to test the reality of certain anomalies in 
the long-range order in solid solutions near the composition Fe,Al, which had 
previously been detected by x-ray diffraction. | Resistance anomalies were 
found in three of the five alloys tested, containing 22-9, 23-5 poe 24-8 at.% 
aluminium, but not in the alloys containing 22-2 and 25-5 at. % aluminium. 
Their significance is briefly discussed. 


§ 1. IyTRODUCTION 


THE experiments to be described here were carried out to confirm the 
reality of an anomaly in the ordering behaviour of iron—aluminium solid 
solutions, originally discovered during an investigation of these alloys by 
x-ray diffraction (Lawley and Cahn 1961). It was found that from several 
of the alloys, the intensity of superlattice lines was higher if the alloys were 
slow-cooled all the way from above the critical ordering temperature 7’, 
than if they were rapidly cooled through a range of temperatures near 7c, 
and then slow-cooled. Figure 1 shows an example of this kind of behaviour. 
It implies that the rapid cooling through 7, (termed transcritical fast 
cooling) works some not readily reversible structural change in the alloy 
which prevents it from reaching full equilibrium order. Only a subsequent. 
anneal well above 600°c restores to the alloy the capacity to attain equili- 
brium order, corresponding to the upper curve in fig. 1. For this reason 
the alloy is said to possess a thermal memory. 

Electrical resistance is a good index of atomic order and therefore, if this 
effect is genuine, then the resistance-temperature relation should be 
sensitive to thermal history. The investigation included alloys containing 
25-5, 24-8, 23-5, 22-9 and 22-2at.% aluminium. X-ray diffraction studies 
of the first three of these had revealed a thermal memory effect, while the 
last two, which gave no superlattice lines and are subject to short-range 
ordering only, were included to test whether these alloys can also show 
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signs of thermal memory. The results of studies of an alloy containing 
19-4 at. % aluminium, which shows highly unusual resistivity properties, 
will be published separately (Feder et al. 1962). 


Fig. 1 


O lOO” +200 3300 4007 S00 ous 
TEMPERATURE CF MEASUREMENT °C 


The intensity ratio of the (111) superlattice line and the (220) fundamental line, 
as a function of temperature, for the 24-8 at. % alloy. @, equilibrium; 
O, after transcritical fast cool. (From Lawley and Cahn 1961.) 


§ 2. EXPERIMENTAL METHODS 


Samples were milled in zigzag form from vacuum-cast and hot-rolled 
sheet about 1mm thick, as for previous resistance studies (Feder and Cahn 
1960, Cahn and Feder 1960). This shape provided an adequate total 
resistance (0-2—0-3 ohm at room temperature) but kept the sample compact 
and thus made it easy to hold its temperature uniform. This was further 
assisted by enclosing the sample, insulated by thin mica sheets, in a slotted 
copper furnace block. 

For resistance measurements, the sample was heated or cooled con- 
tinuously in an electric furnace, in air, at a rate of 1-2°c/min, which was 
slow enough to maintain equilibrium at temperatures down to about 
300°c (Lawley and Cahn 1961). At 700°C, oxidation of the samples was 
negligible; anneals at 800°C were normally carried out in a vacuum, but 
occasional anneals in air at this temperature had no permanent effect on the 
resistance (and therefore on the composition) of the samples. 

Resistance was measured by a double potentiometer method, with due 
precautions to eliminate thermoelectric e.m.f.’s, to an accuracy of 
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~0-2mohm, at intervals of 10-15°c. In view of the complicated shape 
of the samples, no attempt was made to calculate absolute resistivities. 

Since the test furnace was not adapted for rapid cooling, the effect of a 
rapid transcritical cool through 7’, (= 525-550°c according to composition, 
from x-ray data for those alloys which possess long-range order) was 
achieved by first water-quenching the sample from an auxiliary furnace 
held at 600°c or higher, and then placing it in the test furnace, which had 
been preheated to about 500°C, so that it rapidly attained this temperature. 

The programme of heat treatments given to each sample is indicated by 
the inset diagram in the corresponding figure. A sloping line here indicates 
heating or cooling at 1-2°c/min, vertical lines indicate rapid heating or 
cooling. Continuous lines indicate heating or cooling in air, dashed lines 
indicate heating or cooling in vacuum. For each sample some provision 
was made in the design of the experiment for testing the reproducibility of 
the resistance-temperature curves. 


§$ 3. RESULTS 


Figures 2-6 show the results obtained with the five experimental samples. 
Individual experimental points are not shown for lack of space; experi- 
mental scatter was of the order of the thickness of the curves as drawn. 
Reproducibility of similar runs was good throughout. 

Figure 2 shows quite clear evidence of the influence of a transcritical 
fast cool, and of the existence of a thermal memory. The resistance after 
a transcritical fast cool} (curve A) is higher at all temperatures than it is in 
equilibrium (curve B). Moreover, anneals at a series of temperatures up to 
620°c do not remove the effect of the transcritical fast cool; that is to say, 
a thermal memory of the transcritical fast cool persists in the alloy. A 
heat treatment at 800°c however suffices to remove the thermal memory. 
The sequence of experiments was carefully designed to make sure that any 
irreversible changes in the state of the specimen (e.g. oxidation) would have 
been detected. The degree of order corresponding to equilibrium (curve B) 
is higher than after the transcritical fast cool (curve A), in accord with the 
conclusions drawn from the x-ray diffraction data. 

The slow heating of the quenched sample (curve C) shows that ordering is 
sufficiently rapid above about 300°c to maintain equilibrium. The slight 
dip of curve C below curve B is not thought to imply the existence of an 
abnormally high degree of order, for theoretical analysis of the resistance— 
temperature curve (Cahn and Feder 1960) shows that at high temperatures 
there is no assurance than an increase of order will lead to a reduction of 
resistance. . 

In fig. 3, the results for the 25-5at.% alloy show no evidence that 
transcritical fast cooling has any effect, and there is no thermal memory. 
The magnitude of the effect detected by x-raysin this alloy was rather smaller 


+ The term is used to represent the equivalent heat treatment referred to 
above. 


1096 R. Feder et al. on the 


than in the 24:8 at. % alloy (except in the range 400-500°c) and the con- 
sequent effect upon the resistivity might be below the limit of detection. 
It is also possible that the filings used for the x-ray work were slightly poorer 
in aluminium than the solid alloy, from which both the filings and the 
resistance specimens were made. This would help to account for the null 
effect in fig. 3, since the order anomaly rapidly decreases with increasing 
aluminium content (Lawley and Cahn 1961). 


Fig. 2 
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Resistance-temperature curves for 24-8 at. % alloy. 7, is the critical ordering 
temperature and 6, is the ferromagnetic Curie temperatures. 


The 23-5at.% alloy again showed a marked influence arising from a 
transcritical fast cool. Curve a (fig. 4) defines the equilibrium behaviour, 
and curve 6 the perturbation due to transcritical fast cooling. Here, the 
effect of the transcritical fast cool on resistivity shows most just below the 
critical temperatures, and fades away towards room temperature. Since 
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the effect of order on resistivity is always most pronounced at low tempera- 
tures, this behaviour suggests that the difference in state of order induced 
by the transcritical fast cool rapidly diminishes on continued coolingt. 
Nevertheless, its latent influence (i.e. thermal memory) persists, for on 
reheating (curve c) the path of curve b is almost exactly retraced. The re- 
heat to 680°o (end of curve c) suffices to erase the thermal memory, and 
curve d is once again close to curve a. 

It is striking that in this alloy, the resistance during slow heating of the 
disordered sample (curve e) follows the path of the anomalous state of the 

Fig. 3 
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Resistance-temperature curves for 25:5 at. % alloy. 
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alloy rather than the equilibrium path. This peculiar behaviour was also 
found to a pronounced degree during an examination of an alloy containing 
19-4at. % aluminium (Feder et al. 1962). 

The 22-9 at. °% alloy had showed no detectable long-range ordering when 
examined by x-ray diffraction, but the 8°% difference in room temperature 
resistance between quenched and slow-cooled samples (fig. 5) suggest that 
very slight long-range ordering may in fact exist, and have been missed in 
the earlier study. However, it is equally plausible that the resistance 
difference can be explained entirely in terms of short-range order. The 
inversion of resistance as between curves a and d in the range 170—270°c is 
presumably merely a trivial consequence of the smaller temperature 
coefficient of resistivity of disordered alloys (Cahn and Feder 1960), coupled 
with the sluggishness of diffusion, and therefore of ordering, below ~ 250°c. 
It is noteworthy that the ordering, when it comes, leads to a sharp increase 
of the resistance, which lends support to the view that we are concerned 
with short-range rather than long-range order. 


Fig. 4 
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Irrespective of whether the ordering in question here is long range or 
short range, the heat treatment programme used did lead to very pro- 
nounced resistance anomalies. The fact that these anomalies more or less 
‘disappear at 500°c suggests, by analogy with fig. 4, that if there is indeed 
any long-range order, then the critical temperature is near 500°C. 

The resistance of the 22-2 at. % alloy (fig. 6) is virtually unaffected by 
heat treatment; in the quenched condition it has a very slightly lower 
resistance than in the slow-cooled condition. This alloy had shown only 
short-range order when examined by x-rays (Lawley and Cahn 1961). 
However, this short-range order was quite pronounced and it is surprising 
that it has so little effect on the resistance. 
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§ 4, Discussion 

The experiments described above confirm the reality of the transcritical 
fast cool effect, originally discovered during the x-ray diffraction experi- 
ments, and also confirm the persistence of a thermal memory in the samples. 
In detail, there are some points of disagreement with the x-ray results, 
but in evaluating there it has to be remembered that we have little know- 
ledge of the relation between order and resistance for the various alloys ; 
only for the stoichiometric alloy have the details been worked out (Cahn 
and Feder 1960). Indeed, one cannot even be sure whether in the tem- 
perature range just below 7’, a rise in resistance corresponds to an increase: 
or a decrease of the degree of order. For this reason no firm conclusions 
can be drawn, from the resistance data alone, about the extent and even the 
sense of the order anomalies produced by the heat treatments used. It is. 
beyond doubt, however, that the order anomalies exist. Moreover, in 
spite of the slow heating and cooling rates used, the various resistance 
curves did not merge; this confirms the important conclusion, drawn from 
the results of the x-ray work, that the anomalous states of order resulting 
from transcritical fast cooling were metastable and not simply a result of 
sluggish diffusion. Indeed, near 400°C isothermal ordering is complete 
in a matter of minutes (Lawley and Cahn 1960, Feder and Cahn 1960, 
McQueen and Kuezynski 1959), so that slowness of diffusion cannot 
account for the separation of the various curves at 400°C, in figs. 4 and 5. 

A tentative interpretation of the transcritical fast cooling effect and 
of thermal memory, in terms of competition between long or short-range 
order and the supplementary directional order caused by ferromagnetic 
influences, has been proposed elsewhere (Lawley and Cahn 1961), and the 
reader is referred to that paper. 

The present resistance-temperature curves have been compared with 
those published by Bennett (1952), although the compositions studied by 
him were not identical with ours, and he was concerned only with simple 
heat treatments. Allowing for the small differences in composition, the 
forms of the equilibrium resistance curves in the two studies agreed quite 
well. It is of interest, however, that Bennett’s 21-5at. % alloy showed 
a much greater difference in resistance between the quenched and the 
slow cooled conditions than did our 22-2 at. % alloy. Masumoto and Saito. 
(1952) also found a considerable resistance difference between the quenched 
and the slow-cooled conditions of a 22-0 at. % alloy. It can only be pointed 
out, in explanation of these differences, that the fine structure of the alloys 
varies rapidly with composition in this range (Lawley and Cahn 1961), 
and that moreover it is not easy to get reliable chemical analyses of these. 
alloys, so that published compositions should be treated with some reserve. 


ACKNOWLEDGMENTS 
We are indebted to Professor G. V. Raynor, F.R.S. for the provision of 
laboratory facilities and for his support during this work, to Professor J. D. 
Fast of Philips Research Laboratories for casting some of the alloys used, 


High T emperature Electrical Resistance of Iron—Aluminium Alloys 1101 


and to the Secretary of the U.S. Army for a Fellowship awarded to one of 
us (R.F.), during the tenure of which most of this work was carried out. 


REFERENCES 


BENNETT, W. D., 1952, J. Iron St. Inst., 171, 372. 

Cann, R. W., and Frepmr, R., 1960, Phil. Mag., 5, 451. 

Ferper, R., and Cann, R. W., 1960, Phil. Mag., 5, 348. 

Fereper, R., Kramer, K., Levin, E., and Cann, R. W., 1962 (to be submitted to 
Phil. Mag.). 

Lawtey, A., and Cann, R. W., 1961, J. Phys. Chem. Solids, 20, 204. 

Masvumoro, H., and Sarro, H., 1952, Sci. Rep. Res. Inst. Téhoku Univ. A, 4, 321. 

McQueen, H. J., and Kuczynsk1, G. C., 1959, Trans. Met. Soc., AJ.M.EH., 215, 
619. 


ni : 


tA ins vee 
Piaaee ons 


PeetLOs, | 


Observations of Precipitation in Thin Foils of 
Aluminium + 4%, Copper Alloy+ 


By G. THomast 
Metallurgy Department, Cambridge University 


and M. J. WHELAN 
Crystallographic Laboratory, Cavendish Laboratory, Cambridge 


[Received December 12, 1960] 


ABSTRACT 


A high temperature specimen stage has been used to study precipitation in 
thin foils of Al+ 4% Cu alloy directly in the electron microscope. Precipitates. 
a few hundred Angstroms in size are observed to form in a few minutes at 
temperatures ~250° to 300°c. This initial precipitation has been identified. 
by electron diffraction as the 6’ phase, and is found to occur in a well-defined 
single orientation with respect to the matrix. The orientation observed is. 
consistent with one of those of @’ formed in bulk material. The growth at 
higher temperatures of larger @ precipitates at the expense of the 0’ precipitates 
has also been studied. The results show that precipitation in thin foils differs 
from that in bulk material due to the fact that the surfaces act as preferential 
nucleation sites. Most of the precipitates in the thin foils appear to be at the 
surfaces. In aged bulk material preferential precipitation occurs on helical 
dislocations formed from screw dislocations by vacancy climb. The pre- 
cipitates are observed to go into solution at ~500°c. Measurement of the 
rate of dissolution by ciné techniques has enabled an estimate of the diffusion 
coefficient of copper in aluminium to be made. Other observations include 
the formation of regions denuded of small precipitates near grain boundaries 
and large precipitates, and the growth and dissolution of precipitates under a 
fluctuating temperature. 


§ 1. INTRODUCTION 


THE stages of precipitation in the Al+Cu system (low copper content) 
have been extensively investigated in recent years, and much is now 
known about the various types of clusters (Guinier—Preston zones) and 
precipitates formed by ageing at different temperatures. Low temperature 
ageing (~ 100°c) of a quenched alloy generally produces ‘zones’ or plate- 
like clusters of copper atoms on the {100} planes of the matrix (Guinier - 
1942, 1956, Nicholson and Nutting 1958, Nicholson et al. 1959), while 
ageing at higher temperatures (~200°c and above) produces large 
precipitates of #’ CuAl, also on {100} planes, or the equilibrium @ CuAl, 
structure. For a review of the work on Al+ Cu see Silcock e¢ al. (1954), 
Hardy and Heal (1956) and Nicholson ef al. (1959). All of the information 
obtained so far comes from experiments on bulk specimens examined 
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so that only the aged state of the bulk material is studied. A notable 
exception is the work of Castaing and Laborie (1954), who studied 
precipitation both in bulk specimens and in thin foils aged outside the 
electron microscope. These workers have emphasized the differences in 
precipitation behaviour expected by use of the two methods. 

It was evident to the authors that much useful information might be 
obtained from a study of precipitation in thin foils aged inside the electron 
microscope during observation. A comparison of the ageing behaviour 
of thin foils (< ~1, thick) and bulk material can therefore be made. 
The nature of the sites of nucleation of precipitates in thin foils is also 
open to investigation. The experiments reported in this paper have been 
carried out with this object in view. 


§ 2. EXPERIMENTAL TECHNIQUE 

Foils of Al+ 3-7 wt. °% Cu alloy of total impurity 0-003°% and 0-005 in. 
thick were prepared from thicker sheet by several stages of cold rolling 
with intermediate anneals. The foils were solution heat treated by water 
quenching from 540°c. Thin foils suitable for transmission electron 
microscopy were prepared by electropolishing in a phosphoric—chromic 
acid bath (Nicholson ef al. 1958). The specimens were examined in a 
Siemens Elmiskop 1 electron microscope operated at 100 kv using the 
high temperature stage described by Whelan (1960). The accuracy of 
temperature measurement in these experiments is roughly estimated as 
+30°c. For dynamic experiments photographs were recorded on 16 mm 
ciné-film in a camera mounted outside the viewing window of the electron 
microscope. 

§ 3. EXPERIMENTAL RESULTS 
3.1. Hxamples of Precipitation Patterns Observed 

The experiments on the thin foils were carried out by ageing at 
temperatures ~ 250°c and above, where the @’ and @ phases are expected 
to form in the bulk material. Selected area electron diffraction patterns 
indicated that both these phases are also formed in the thin foils (see 
§3.3). Figures 1 to 7 and fig. 13} are typical micrographs obtained from 
thin foils of the alloy aged at various temperatures showing the precipitation 
patterns observed. Figure 1 (a) and (b) is a micrograph of an area where 
the foil surface is approximately (112). Figure 1(a) taken after ageing 
for a few minutes at ~ 270°C, shows rounded or irregular-shaped precipi- 
tates of average diameter ~ 500A with an areal density of ~ 101° em-, 
while fig. 1(b) shows the same area a few minutes later at a temperature 
of ~340°c. A marked coarsening of the precipitates is apparent. The 
density has decreased by a factor of three or more while the average 
diameter is greater than ~ 10004. This coarsening is the result of growth 
of larger precipitates at the expense of smaller ones which disappear. 
Many of the small precipitates on fig. 1 (a) have shrunk and disappeared 
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on fig. 1 (6), while others (e.g. A, B and C on fig. 1 (a) and (b)) have grown 
and assumed angular shapes. The small round or irregular-shaped 
precipitates first observed are hereafter referred to as the ‘initial 
precipitation’. Electron diffraction patterns and dark field experiments 
($3.3) have shown that the initial precipitation is the 6’ phase of CuAl.,. 
The larger angular-shaped precipitates on fig. 1(6) have been similarly 
identified as the @ phase. Figures 3 and 4 are micrographs of areas of 
approximately (110) orientation where the contrasting morphology of 6’ 
and @ precipitates is evident. 

The change from the small @’ precipitates to the angular 0 precipitates 
occurs very slowly below about 300°c. Electron diffraction patterns 
($3.3) show only @’ spots below this temperature for ageing times of 
several minutes. However, many of the @ precipitates in fig. 1(b) can 
be recognized as rather large precipitates on fig. 1(a). In some cases 
dark field micrographs (§3.3) have shown that even at temperatures 
~250°c these larger precipitates are not the 6’ phase. It appears 
therefore that the @ phase is present even at the lower temperatures 
although the quantity is insufficient to give appreciable diffraction. It 
is also unknown whether 6 is formed by transformation from @’ or whether 
it nucleates independently. The growth of @ precipitates at the expense 
of & precipitates occurs fairly rapidly above ~350°c. 

A considerable amount of evidence has accumulated which suggests 
that much of the initial precipitation and of the subsequently formed 
4 phase is nucleated preferentially at the foil surfaces. The facts are as 
follows : 

(a) Features of the initial precipitation (e.g. areal density, size and shape 
of precipitates) depend markedly on orientation, i.e. on the crystallographic 
nature of the foil surface, and vary from grain to grain in the same 
specimen. When the foil surface is near (100) the initial precipitation is 
usually slow to appear, very irregular, and is difficult to see clearly in the 
early stages. Figure 2(a) is a micrograph of a foil with a (100) surface 
aged for about 10 min at ~270°c. Two types of precipitate are visible. 
The first includes the large 0’ platelets on {100} viewed edge-on; in this 
case they are known to have nucleated in the interior on dislocations 
previously observed but out of contrast in fig. 2 (a) (see § 3.2). The second 
type (also 6’ phase) includes the irregular precipitates of poorer visibility 
which are thought to lie at the surface. This type is always characteristic 
of (100) areas not containing internal centres of nucleation such as 
dislocations. Regions of near (112) and (110) orientations are shown 
in fig. 1 (a) and (6) and in figs. 3 and 4 respectively. In these micrographs 
the 0’ precipitates assume a more rounded shape. coe ane 

(b) A study has been made of the early stages of the initial precipitation 
in an attempt to detect any precipitate density variation with foil thickness, 
as expected if the precipitation is a volume effect. No noticeable effect 
has been observed except in thin regions (see below). This suggests that 
the density of precipitates depends on the surface area of the foil. 
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(c) Preferential precipitation has been observed at inclined grain 
boundaries running through the foil. Figure 5(a) and (b) shows such a 
boundary at A. The precipitates are observed to form along the edges. 
of the boundary B and C corresponding to the top and bottom surfaces 
of the foil. 

(d) In areas of (100) orientation an extended type of precipitate pattern 
is produced at temperatures ~330°c. Figure 2(b) shows an example 
formed in a different area of the same grain as in fig. 2(a). This extended 
pattern forms suddenly from the irregular precipitates of fig. 2(a). It 
has been identified as the 6 phase. These extended patterns are thought 
to be produced by coalescence of 6’ plates lying on the (100) surface plane. 
Detailed inspection shows that the pattern consists of two overlapping 
patterns, one of which shows bend extinction contours (Heidenreich. 
1949) while the other shows no extinction contours. These probably 
represent extended precipitates on opposite surfaces of the foil. Further- 
more dislocations have been observed to move in an unhindered manner 
through such areas. 

(ec) When @ precipitates dissolve at ~500°C (§3.4), marks are often 
observed on the image covering the areas previously occupied by 
precipitates. These marks are particularly noticeable if the foil has. 
been held at high temperature for several minutes before the solution 
treatment. An example of these marks is visible in fig. 13 (e) around the 
large precipitate. The marks may be produced by increased oxidation 
of the foil at the site of the precipitate. Whatever the origin of the marks 
they are much more likely to be a surface than a volume effect since the 
interior of the metal would be homogeneous above the solution temperature. 

_ It therefore seems probable that in thin foils the surfaces provide sites. 
for preferential nucleation of precipitates and that many of the precipitates 
observed are at the surfaces. This is expected on the hypothesis of Seitz 
(1952) by which solute atoms are associated with vacancies and move 
with them. Quenched-in vacancies would be expected to diffuse to the 
surfaces of the thin foils on ageing and deposit copper at the surfaces, 
thereby providing nuclei for further precipitation. Surface diffusion of 
copper atoms would also account qualitatively for the rapid rate of changes. 
which take place suddenly at comparatively low temperatures, e.g. for 
the change between fig. 2 (a) and (6) following a small rise of temperature. 
Precipitation in the interior of the foil is sometimes observed; e.g. the 
6’ plates viewed edge-on in fig. 2(a) are undoubtedly nucleated in the 
interior. However, it is believed that these precipitates are characteristic 
of nucleation at dislocations (§3.2) and would be absent if dislocations. 
were not present. Dislocations are usually observed in such cases. 

In thin regions of the specimen near the edges of the foil, precipitates. 
are often observed which run right through the foil. These precipitates, 
are usually separated from the surface precipitation patterns of thicker 
regions by a zone denuded of surface precipitation. An example is shown 
in fig. 6(a); the precipitates have been identified as the @ phase. The. 
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boundary surfaces between matrix and precipitates are usually nearly 
perpendicular to the foil surface, so that the projected widths of the 
boundary planes are small. However, dark field images in the strong 
matrix reflections have shown that the areas occupied by these precipitates 
behave like holes in the foil, i.e. they scatter a negligible number of electrons 
into the matrix reflections. This proves that the precipitates run right 
through the foil. Figure 6(b) shows a similar 6 precipitate in a thicker 
region where the inclined boundary plane with thickness interference 
fringes is visible. Several different orientations of these precipitates 
occur in the same region but no systematic study has been made. 
Regions denuded of precipitation are frequently observed near grain 
boundaries. Figure 5(a) and (6) is a typical example at a triple grain 
boundary junction. Preferential precipitation has occurred at the grain 
boundaries leaving a denuded zone ~4y wide. Typically the denuded 
zone is observed to increase in width with ageing time, until ultimately 
most of the initial precipitation has disappeared leaving behind large 
precipitates at the grain boundaries and in the interior of the grains. 
Denuded regions also form around large precipitates in the interior of 
grains. Figure 7 (a) and (b) is a micrograph of the same area of a foil of 
unknown orientation. Figure 7 (a) taken after 30 min at ~250°c shows 
the initial precipitation. Figure 7(b) shows the rather striking effect 
observed in this area after raising the temperature a few degrees. The 
initial precipitation coarsened rapidly and larger precipitates such as: A 
and B suddenly formed leaving denuded regions around them. The large 
precipitates appear to run right through the foil with almost perpendicular 
boundaries. Ultimately all the initial precipitation disappeared leaving 
only a few of the large precipitates. Although the large precipitates were 
not identified in this case they are undoubtedly the 6 phase. All these 
changes can be easily followed visually in the electron microscope. 


3.2. Comparison of Precipitation in Thin Foils and Bulk Material 


To test whether the precipitation observed in thin foils is a surface 
phenomenon, the following experiment was performed. A thin foil was. 
prepared from quenched bulk material. The electropolished foil is 
generally only suitable for transmission microscopy in local areas near the 
edge of the foil; eisewhere it may be several microns thick and can be 
considered as bulk material. The specimen (i.e. thin parts and thick parts 
together) was then aged for a few minutes in high vacuum at 270°C to 
simulate conditions occurring on the electron microscope high temperature 
stage. Thin specimens were cut from the edge and were examined in the 
electron microscope. The micrographs were similar to those obtained 
from thin foils aged in the microscope at about this temperature (e.g. 
compare fig. 8 with fig. 1(a)). The specimen was then electropolished 
again to remove several microns, i.e. a thin foil was prepared from what 
may be considered as aged bulk material. The micrographs obtained 
from this specimen were noticeably different and showed much larger 


4B2 


1108 G. Thomas and M. J. Whelan on 


precipitates of 6’ phase lying on {100} planes in Widmanstatten-type 
patterns similar to those previously observed by Thomas and Nutting 
(1955) using oxide replica techniques (see, for example, fig. 9). This 
experiment is fairly convincing proof that much of the initial precipi- 
tation in thin foils is characteristic of the surfaces of the specimen, and 
that thin foils differ appreciably from bulk material in ageing behaviour. 

The precipitation patterns observed in aged thick specimens subsequently 
thinned as in fig. 9 are probably produced by preferential precipitation at 
dislocations. Itisnow known (Thomas and Whelan 1959) that quenched -in 
vacancies in Al + 4°/, Cu alloys migrate to screw dislocations and cause them 
to climb into helices. The segments of a helical dislocation have consider- 
able edge component and it seems probable that if copper atoms migrate 
to the dislocations with vacancies they will segregate at the compression 
side of the dislocation, thereby relieving the strain and forming nuclei 
for further precipitation. Therefore, on ageing, a helical dislocation may 
give rise to the arrays of precipitates lying on {100} planes visible in fig. 9. 
Figure, 10 shows an example of precipitation at helical dislocations in a 
bulk specimen aged for 2 min. at 270°c. This micrograph shows an earlier 
stage of precipitation than that visible in fig. 9. In this case the small 
precipitates A formed at the helical dislocations seem to be partially 
coherent and give rise to diffraction contrast effects due to strain in the 
surrounding matrix similar to those observed by Nicholson and Nutting 
(1958). The contrast effects are very similar to those observed at dislocation 
loops in quenched aluminium (Hirsch ef al. 1958). The dislocation line B 
connecting the precipitates is clearly visible in fig. 10, but elsewhere (e.g. 
at C-C) it appears to be out of contrast (Hirsch ef al. 1960). 


3.3. Identification of Precipitates 


The selected area diffraction technique has been used to identify the 
precipitates at various temperatures. The structure of the @ phase of 
CuAl, has been studied by Friauf (1927) and by Bradley and Jones (1933). 
The Bravais lattice is body-centred tetragonal with a= 6-054 A, c=4-8644 
(Bradley and Jones 1933). These lattice parameters have been used to 
compute d-spacings which were compared with observed spacings measured 
from electron diffraction patterns. The structure of the 0’ phase is also 
based on a body-centred tetragonal lattice with a=4-:044, c=5-8A 
(Silcock et al. 1954, Silcock and Heal 1956). The a parameter is practically 
identical with that of the matrix. It was found in all cases that the electron 
diffraction spots could be interpreted by assuming that one or other (or 
both) of the above phases was present. 

The dark field technique has been used to distinguish between the two 
phases when present together. A dark field micrograph taken with one 
of the known 6’ or @ reflections serves to show up the location of many of 
the corresponding precipitates in a given area. In this way it has been 
possible to show that the large or angular shaped precipitates in figs. 1 (b), 
2 eh 3 and 4 are the @ phase and that the initial precipitation is the 6’ 
phase. 
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Selected area diffraction patterns have shown that the initial precipitation. 
formed on the first ageing invariably consists of precipitates of the same 
orientation. The orientation of 6’ observed is that in which the c-axis 
makes the smallest angle with the normal to the foil. The orientation 
relationship usually holds over the whole grain. Figures 11 and 12 (a) 
are typical examples of electron diffraction patterns of the initial precipi- 
tation for foils with near (011) and (001) surfaces respectively, showing 


Orientations of precipitates 


Number of 


ae Approximate A i ‘i 
Precipitate pp : er te reas pproximate orientation 

P foil surface | “* ice of precipitate 

6’ (001) several [001],- parallel to [001 ] yt 
[100] common to both. 

0 (O11) 3 [011],- parallel to [011], 
[100] common to both. 

g’ (112) 3 [111],, parallel to [112], 
[110] common to both. 

0 (O01) several [001], parallel to [001], 
[100], parallel to [110]q, 

6 (O11) 1 [113], parallel to [O11] y_ 
[110], parallel to [100], 

6 (O11) 1 [111], parallel to [011],, 


[110], parallel to [100], 


+t M denotes the matrix. 


single crystal type patterns from 0’ (selected area diameter ~1y). The 
table summarizes the results for three different foil surfaces. The 
orientations are approximate; the first parallelism quoted in each case 
may be in error by as much as ~10°. All the orientations of @’ in the 
table are consistent within the error with the well-known orientation of 6” 
in the bulk material, i.e. the tetragonal axes of 6’ are parallel to the cubic 
axes of the matrix. The orientation of the 6 phase formed after the 
disappearance of 6’ has not been extensively studied except for foils 
with near (100) surfaces. In this case the observed orientation (see table) 
is the same as that reported previously by Guinier (1942) (Guinier’s group 1 
orientation) for 6 formed from 6’ in the bulk material. It is possible that. 
the last two cases of the table for foils with near (011) surfaces represent. 
this case also within the experimental error. Figure 12(b) shows the 
diffraction pattern of the 6 phase formed from the 6’ phase of fig. 12 (a). 
Both these patterns were obtained from areas similar to those shown in 


fig. 2(a) and (6). 
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Little can be said at present about the habit plane of the initial 
precipitates. In the case of foils of (100) orientation (fig. 2 (a) and (0)), 
the habit plane is believed to be the (100) surface. In some other 
orientations the 6’ precipitates appear to have grown into the foil on an 
inclined plane which might be {100} after nucleating at the surface. The 
precipitates in fig. 1(b) at D, E and other places lend support to this view. 
These rounded precipitates appear elongated in the [110] direction as if 
they were roughly circular discs viewed at an angle. In this case this 
direction is parallel to the trace of the (001) plane perpendicular to the 
c-axis of 6’ on the foil (112) surface One edge of each precipitate (parallel 
to this direction) is rather irregular suggesting that this is the intersection of 
the (001)-habit plane with the irregular surface of the foil. Similarly in 
fig. 4 the large @ precipitates probably have the inclined (100) plane as 
habit plane and run right through the foil. However, in many other cases 
there is no clear indication of a habit plane. 


3.4. Dissolution of Precipitates at High Temperatures 


When the temperature of a thin foil containing precipitates is raised 
above ~400°C coarsening occurs rapidly, and ultimately at ~500°c 
all the precipitates go into solution. Figure 13 is a typical sequence of 
micrographs showing how the originally angular-shaped @ precipitates 
become rounded-off and shrink as the temperature is raised and eventually 
disappear. The foil surface is (100) and the precipitates in this case are 
probably thin plates lying at the foil surface. If the temperature is 
suddenly raised and held constant, the rounded precipitates are observed 
to shrink slowly at first when their diameters are large and then more 
rapidly as their diameters decrease. 

The rates of solution of the last few remaining precipitates have been 
studied using ciné techniques. It has been found in three cases out of 
twelve examined that d(r)/dt is approximately constant over the range 
of r (r< ~ 20004) (r is the precipitate radius at time t). Figure 14 shows 
a plot of r? versus t in such a case. In the other cases however, similar 
plots are irregular or show a roughly linear variation for large particle 
radius but have negative values of d?(r?)/dt? for small particle radius. The 
linear variation can be understood on the basis of a diffusion-limited 
dissolution process. Several authors (Zener 1949, Wert and Zener 1950, 
Frank 1950, Ham 1958) have considered the reverse case of growth. For 
a precipitate of negligible initial dimensions in an initially uniform super- 
saturation, where the growth rate is determined by diffusion and not by 
any interfacial phenomenon, these authors have shown that the solution 
of the time-dependent diffusion equation leads to r? = kDt, where kis the rate 
constant and D is the coefficient for diffusion of material to the precipitate. 
Boundary conditions at large distances from the precipitate are assumed. 
Ham (1958) has also considered growth of a spheroidal precipitate (rod-like 
or plate-like) and has shown that a similar solution exists, the eccentricity 
of the precipitate remaining constant. The growth constant k depends on 


Precipitation in Thin Foils of Al+4%, Ou Alloy 1 


the shape of the precipitate. Frank (1950) has given graphs from which & 
may be determined for the cases of spherical, cylindrical and one-dimen- 
sional diffusion in terms of the ratio (Ps—Ppe)/(pe—ps) Where pe, ps and pe 
are the concentrations of diffusing material inside, at the surface, and at 
a large distance from the precipitate. For small values of this ratio in the 
spherical case the rate constant k is 2(ps—pe)/(pc—ps), a result which can be 
obtained by an approximation using steady-state diffusion theory. 


Fig. 14 
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Plot of (radius)? versus time for the last few seconds of the lifetime of a shrinking 
precipitate. Measurements taken from ciné-film. 


As far as the authors are aware, no exact solution of the dissolution 
problem has been given. When the rate constant is small however, it is 
expected that it will be the reverse of that in the growth case, i.e. 

d(r?) 
dt 


=-—kD, eens Oe eo atm (1) 


ease G. Thomas and M. J. Whelan on 


The applicability of the above theory to precipitates in thin foils is: 
questionable. The effects of the surfaces have been neglected. Even if 
the foil is thick compared with the precipitate diameter, the correct 
value of k is obscure since it depends on the precipitate shape. The influence 
of surface diffusion on dissolution of precipitates at the surface of the foil 
is also unknown. It is not surprising therefore that a linear variation is 
the exception rather than the rule. However, if it is assumed that in the 
cases where a linear variation is observed the precipitates are small spheres. 
in the interior of a relatively thick foil it is possible to estimate the diffusion 
coefficient of Cu in Al from eqn. (1). Spherical precipitates satisfying this 
condition will probably be formed from @ precipitates which initially run 
right through the foil es in fig. 6(b). For a sphere we have 

Daas y Pere tlt aa 

(ps— pe) dt 

Pc; pe and ps can be taken approximately as the atomic per cent concentra - 
tions. Since nearly all the copper is in solution pe is about equal to the 
concentration of copper in the alloy (1-6 at. %). ps can be neglected in 
comparison with pe in the numerator of eqn. (2). pe is about 32-4 at. % 
ps—pe can be obtained from the slope of the phase boundary on the 
equilibrium diagram if A7’ (the temperature increase above the equilibrium 
point at 1-6 at. % copper) is known. ps—pe is about 1-3 x 10-AT at. %. 
Putting d(r?)/dt equal to ~—5x 10-1 cm?sec from the slope of 
fig. 14, we obtain D~6 x 10-8/AT' cm? sect. For AT ~30°%c, D is 
~2x10-*%cm?sec, which is slightly larger than, but of the same order 
of magnitude as the values of the diffusion coefficients for copper in 
aluminium above 500°c reported in the literature (see Smithells 1955). 
The observed rate can therefore be understood on the basis of a diffusion - 
limited process. With further refinement the electron microscope method 
might provide a suitable method of measuring diffusion coefficients in | 
other alloy systems. However, it appears to be limited at present by 
uncertainties in the theoretical interpretation, besides the purely practical 
difficulty of estimating the temperature accurately. 


3.5. Some Further General Observations 


Several other interesting observations have been made on the dissolution 
process. An experiment was performed in which precipitates were observed 
while the temperature was successively raised and lowered by about 20°c 
from ~500°c. In this way precipitates could be made to diminish in size 
(raising the temperature) or grow (lowering the temperature) at will, 
i.e. there is no difficulty in growth provided a site for precipitation is still 
available. However, if the precipitates are allowed to dissolve completely 
at 500°0, re-precipitation does not occur when the temperature is lowered 
by ~20°c. In order to cause precipitation after dissolution the temperature 
had to be lowered by about 100°c below the limit of solid solubility. The 
precipitates then formed were quite different morphologically, consisting 
mainly of extended patterns probably lying at the foil surfaces. 
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Most of the effects described in this paper have referred to the first 
ageing treatment. It has also been possible to quench-in the solid 
solution in thin foils from above the solubility limit}. From successive 
precipitation—dissolution—quenching—+re-precipitation sequences carried 
out on the same specimen in the microscope it was found that the second 
and subsequent precipitation differed considerably from the first precipi- 
tation. In general much less @’ is produced by low temperature ageing 
during the second precipitation. Angular precipitates characteristic of 
the @ phase are often formed during the second precipitation. Figure 5 (a) 
and (b) is a micrograph (at low magnification) of the same region at ~ 300°C 
during the second and third precipitations. The density of the precipitates 
is lower than that during the first precipitation while their size is greater. 
Diffraction patterns of typical regions during the second precipitation 
showed spots belonging to 6’ and @ of several orientations in contrast to. 
the well-defined single orientation observed during the first precipitation. 
Jt was also found that the specimen has no memory regarding nucleation 
sites, i.e. by observing precipitates formed at a given temperature, 
dissolving them, quenching and by re-precipitating at the same tempera- 
ture, keeping the same region in view, it was found that the new 
precipitates appeared to nucleate quite randomly (e.g. compare figs. 5 (a). 
and (b)). The differences between the first and subsequent precipitations. 
are not understood at present and further work is required on this point. 
Finally, it may be stated that the observations described in this paper are- 
most conveniently demonstrated by ciné-photography. 


§ 4. CONCLUSIONS 

1. The experiments have shown that it is possible to observe directly 
precipitation in thin foils of Al+ 4% Cu alloy. 

2. Precipitation in thin foils generally appears to differ from that in. 
bulk material. Precipitation seems to occur preferentially at the foil 
surfaces, and this is probably due to the fact that solute atoms are associated 
with and migrate to the surfaces with quenched-in lattice vacancies. 
In contrast to this, precipitation in bulk specimens occurs preferentially 
at helical dislocations formed from screw dislocations by vacancy climb. 
Again this can be interpreted in terms of migration of solute atoms with 
vacancies. 

3. During the first ageing treatment the first precipitates to appear 
at temperatures ~ 250°-300°c are of the 6’ phase. Usually only one. 
orientation of 6’ is observed over a whole grain. The orientation observed 
is consistent with one of the three possible orientations of #’ in bulk material. 

4. At temperatures ~350°c the small 6’ precipitates disappear, while- 
larger angular precipitates of the @ phase are formed. In a few cases. 
examined the orientation of @ observed is consistent with one of the 
orientations of @ formed from 4’ in bulk material. 


pee ee ee ES 
+ Switching off the heating current to the high temperature stage generally 
gave a sufficiently fast cooling rate to effect this. 
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5. At ~500°c the precipitates go into solution. In a few cases d(r?)/dt 
has been found to be constant. This result can be interpreted in terms of 
diffusion-limited dissolution and the diffusion coefficient can be roughly 


estimated. 
6. Quenching and re-precipitation sequences have shown that subsequent 
precipitation differs considerably from the first precipitation. 
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ABSTRACT 


Previous measurements of the amplitude-dependent internal friction of 
NaCl at 90 ke/s and room temperature (Whitworth 1960) have been extended 
to study the corresponding changes in Young’s modulus. It is concluded that 
the motion of dislocations giving rise to the amplitude-dependent damping 
lags in phase behind the stress by about 30° to 60°. Previously described 
changes in damping produced by high amplitude vibration are associated with 
corresponding changes in the modulus defect. 


AN account of some experiments on the internal friction of single crystals 
of NaCl at 90kc/s and room temperature has already been published 
(Whitworth 1960, hereafter called I), and that paper described how the 
damping could be changed by vibrating the specimens at high amplitudes. 
There is always some uncertainty in the interpretation of such dislocation 
damping measurements unless the change in elastic modulus associated 
with the internal friction is known, but the apparatus described in I was 
not suitable for measuring these modulus changes. The present work was 
therefore undertaken to clarify this point. 

The composite resonator described in I was again used, but in this case 
it was driven by a crystal stabilized oscillator whose frequency could be 
controlled to 1 part in 10°. . Both the internal friction of the specimen, 64, 
-and the resonant frequency of the system were measured as a function of the 
maximum strain amplitude, ¢, in the crystal; from the frequency measure- 
ments (AH/E),, the relative decrease of Young’s modulus below its value 
at the lowest amplitudes, wasdetermined. Inall specimens this amplitude- 
dependent modulus defect was found to be closely related to the amplitude- 
-dependent part of the decrement 6,,= (5 —59), where 5, is the decrement at 
low amplitudes. The modulus defect, (AH/E),, associated with 6) was 
not measured in these experiments. 

The figure shows the results of measurements of 6, and (AH/E), as 
functions of ¢ for two specimens deformed in compression by 0-2% and 1% 

respectively. 6, and (AH/Z), are plotted on a common logarithmic 
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scale, and the values of 8, are given below the figure. The precise nature of 
the curves obtained varies somewhat from one specimen to another, but the 
principal point to note is that 7=6,/(AH/Z), varies only slightly for quite 
large changes in §,, and (AH/Z),. The table shows some typical values of r 


exiGge 


Amplitude-dependent decrement and modulus defect for two crystals. 
Specimen deformed 0:2%: O, 8,; @, (AH/E),; 89=1:7x 10-+. 
Specimen deformed 1% : A, 8,; A, (AH/E),; 8)=6-1 x 10-. 


for three specimens deformed by differing amounts; the values quoted are 
accurate to about + 25%, and show that while r is effectively constant in 
the most heavily deformed crystal, it rises with ¢ in the more lightly 
deformed specimens. These observations are consistent with the very 
approximate results quoted in §3.1 of I, and confirm the approximate. 
values for the amplitude of motion of the dislocations calculated in that 
paper. 

A knowledge of both AZ/E and 5 leads to information about the motion _ 
of the dislocations under the oscillatory stress, and if this motion is assumed. 
to be sinusoidal, it can be shown that it lags behind the stress by a phase- 
angle equal to tan-t [(3/7)/(AH/E)] (cf. Nowick 1953). At very low ampli- 
tudes the results of Gordon and Nowick (1956) and of Bauer and Gordon. 
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(1960) show that 5,/(AH/E), ~ 0-2 in their specimens, and this indicates that 
the motion of the dislocations giving rise to the amplitude-independent 
damping is only a few degrees out of phase with the stress. In the ampli- 
tude-dependent range, however, r is an order of magnitude larger, being 
~ 7, and the additional motion of the dislocations giving rise to this damping 
therefore lags in phase behind the stress by 30° to 60°. 

In any model in which the dislocation displacement is a single-valued 
function of the stress there will be no phase lag, but the introduction of a 
relatively small dislocation damping term can lead to small phase differences 


Typical values of r=5,/(AZ/EZ),, 


Deformation of r ? 
specimen (°%) at e=10= athe Oe? 


0-2 2-0 6-0 
1-0 2-0 3:4 
35 They, 2°3 


such as those observed at low amplitudes. At the other extreme, if the 
dislocation velocity is dependent only on the applied stress there will be no 
modulus defect, and the motion of the dislocations will lag 90° behind the 
stress. The present results show that the model appropriate to the 
amplitude-dependent damping of NaCl lies somewhere mid-way between 
these two simple cases. 

It was shown in I that if a specimen is deformed and then annealed at 
100 to 200°c its internal friction is lower than its value determined 
immediately after the deformation; this was attributed to pinning of the 
dislocations during the anneal. It was further observed that vibration of 
such an annealed specimen at fairly high amplitudes increased the values 
of §,, subsequently observed at lower amplitudes, and this was explained 
as being due to the pulling of dislocations away from pinned sites to which 
they do not return. Corresponding effects on the modulus defect would be. 
expected, and they have now been observed, the values of (AH /E), being 
reduced by annealing and in part restored by high amplitude vibration. 

Another effect described in I was a decrease in the internal friction of 
deformed specimens vibrated at strain amplitudes above about 3 x 10~°. 
Measurements of the modulus defect before and after this high amplitude 
vibration treatment have shown that it too is decreased, thus confirming 
that the effect is due to some process which prevents the dislocations from 
moving under the oscillating stress. This * pinning * is associated with 
the production of ‘debris’ by the oscillating dislocations (Davidge and 


Whitworth 1961). 
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ABSTRACT 


In Al-1-7 at. % Cu the critical resolved shear stresses, tension stress—strain: 
curves, and the dependence of the flow stress on strain rate were determined: 
from 4-2 to 273 or 373°K for the following structures: supersaturated « solid 
solution, «+GPI, «+GPII, «+ 6’, and «+6. The difference between the- 
CRSS at 4-2 and 273°x is large for the first two structures, but small with 
GP II, & or @ precipitates. With «, «+GPI1, and «+GPII, a single slip 
system is operative over a significant portion of the stress-strain curve, and the 
rate of work hardening is about the same as that for pure Al. With 6 or 6’ 
multiple slip is immediately detected beyond the yield stress, the work-harden- 
ing rates are initially extremely rapid, and fall to small values after some per 
cent strain so that a plateau occurs in the stress-strain curves. 

For supersaturated «, the strength is attributed to clustering in the solid 
solution. 

In the cases of «+GPI and «+GPII, the zones must be sheared 
during plastic deformation and the strengthening arises from the atomic 
scrambling which occurs and the local strains near a zone which must be 
overcome when a dislocation moves through a zone. The variation of 
critical resolved shear stress with temperature and the strain-rate data are 
discussed in terms of the theory of thermally activated slip. Activation 
energies computed from the data are larger than those predicted on the basis 
of atomic scrambling alone. 

With @ or @ precipitate dislocations bend between the particles. 


§ 1. IyTRODUCTION 


Four distinct stages of precipitation occur in Al base Cu alloys: 
Guinier—Preston zones of the first and second kinds, 6’, and @, the last. 
being the stable precipitate. These have been extensively studied by 
x-ray diffraction (for reviews see Heal and Hardy 1954 and Guinier 1959) 
and the electron microscope (Nicholson and Nutting 1958, Nicholson et al. 
1958-59). Commercial alloys based upon this system are heat treated to. 


een ne eee 
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contain substantially GP I or GP IL; the presence of substantial amounts 
of 6 or 6’ corresponds to overageing. While these alloys have found 
important uses for many years, many uncertainties remain concerning 
the influence of the various precipitates on mechanical properties and the 
sources of strengthening are not well understood. Several years ago a 
programme was undertaken to determine the fundamental strength 
properties of Al base Cu alloys with various precipitates. The initial 
phase (Kelly and Chiou 1958) reported, for Al-1-8at.% Cu with GPI 
zones and for the supersaturated solid solution, the variation of the 
critical resolved shear stress, and the flow stress with temperature over 
the range 77 to 373°K. This work has now been extended to include all 
four precipitates with an alloy of 1-7at.% Cu. The temperature range 
of study has been lowered to 4:2°x. The effect of strain rate (10~° to 
10-8 sec) has also been investigated. 

Concurrently and independently Dew-Hughes and Robertson (1960) 
determined the room temperature plastic properties of single crystals of 
Al base Cu alloys with various precipitates as functions of Cu content. 
The experimental situation is now much clearer. 

The purpose of the present paper is to report the experimental findings 
concerning plasticity of single crystals as function of structure, tempera- 
ture and strain rate and to discuss these in terms of existing theory in an 
effort to understand the causes of the strength. 


§2. EXPERIMENTAL 
2.1. Sample Preparation 


An alloy of 1-7at.°% Cu in the form of 0-260 in. diameter forged rod 
was provided by the Aluminum Company of America. The impurities 
reported by Alcoa are 0-002°%% Fe, 0-002°% Si and 0-004°% Mg. The forged 
rod was rolled to a thickness of about 0-060 in. and machined into tensile 
specimens of about lin. gauge length and about 0-1 in. wide in the gauge 
length. Single crystals were grown by first annealing at 550°c for 5 min, 
then plastically straining about 1-5% in tension and finally lowering at 
1 cm per hour into molten salt (Houghton 430) at 550°c. The orientations 
of the crystals were determined to +2° by the back reflection Laue 
technique. 

Heat treatment consisted of holding the specimens 15 hours at 550°c 
in molten salt, quenching into water at 25°c, followed immediately by an 
ageing treatment. Ageing treatments to 220°c were done in a circulating 
hot-air furnace and those above 220°c in molten salt. Resolution of GP I 
was accomplished by holding specimens in a bath of boiling benzyl aleohol 
(205°c) for 6 min. 

2.2. Plastic Deformation 


Testing was carried out with an Instron tensile testing machine at a 
strain rate of 10-4sec~' except where noted. During deformation the 
temperature of the crystal was controlled by immersing it in liquid helium 
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for 4-2°K, liquid nitrogen for 77°K, dry ice and acetone for 198°x, ice and 
water for 273°x, and boiling water for 373°R. 

In measuring the effect of strain rate on the flow stress at a given 
temperature, the cross-head motion was stopped, the drive gears changed 
to those for the desired speed, and the cross-head motion was then resumed. 
During the period when the test was stopped the stress was maintained 
manually. 

For deformation at 4:2°K a liquid helium cryostat similar to that of 
Blewitt et al. (1956) was constructed. A photograph of the parts is shown 
in fig. 1. Holes drilled through the wall of the compression tube permitted 
liquid helium to be in direct contact with the specimen. 


2.3. Shear Modulus 
The shear modulus, @, was determined by the resonant dynamic method, 
utilizing a cylindrical quartz crystal vibrating in the torsional mode. 


The values of G at various temperatures for specimens with GP I, GP II 
and 6’ precipitates are given in fig. 2. 


Fig. 2 


o 


2.9 x an 
e GPI 


2.8 Belg 4 GPL 


Dat, x10" 


2.6 x10" 


RIGIDITY MODULUS in DYNES/CM® 


2.5 x10"+*—105 150 200 250 300m 350 


TEMPERATURE in °K 


Shear modulus as function of temperature for Al-1-7 at. % Cu with various 
precipitates. | 


2.4, Structure Determinations with X -rays 


For the identification of the precipitates after ageing, specimens of the 
same thickness as the tensile specimens were solution treated, quenched 
and aged simultaneously with the tensile specimens. They were then 
etched to about 0-010 in. thickness in a 10% NaOH solution and cleaned 
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+ Figure 1 is shown as a plate. 
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with concentrated HNO,. These specimens were oriented in a three-circle 
goniometer so that {100} planes were normal to the incident beam. 
Transmission Laue photographs were obtained using Mo radiation (40 kv, 
10 ma, 8 hrs). 

Most of our ageing treatments were based upon the comprehensive X-ray 
study of Silcock et al. (1953), to which reference should be made for details 
of the structures of the precipitates. Our x-ray results for the most part 
confirmed theirs. Guinier—Preston zones of the first kind (GPI) were 
produced by ageing for two days at 130°c, which gives zones of about 
1004 diameter. For Guinier—Preston zones of the second kind (GP II) 
a number of ageing treatments was used. One day at 190°c results in a 
mixture of GPII and some 6’, the amount of @’ being variously estimated 
as 15% (Guinier 1952) and 30% (Silcock et al. 1953) of the maximum 
possible. In the present work, photometry of the 6’ spots indicated 13%, 
of the maximum amount of 6’. The thickness of the 6’ plates is 80 A (2bid.). 
1000 hours at 130°c gives GP IT plus 5 to 10% of the possible 6’; the zones 
are 20-30 4 thick and 500 A in diameter while the 0’ platelets are morethan 
1004 thick (ibid.). Five hours at 190°C gives only GPIT but the 
precipitation is not complete. In the present research some GP II was 
found after 4092 hours at 130°c. Relative intensity measurements 
indicated the presence of only about 35°, of the maximum amount of 
&’. Sileock et al. (1953) report no GPII beyond 2600 hours at this 
temperature. 

6’ was produced by ageing one day at 240°c. This gives a platelet. 
thickness of 300A; no other precipitate is present (ibid.). @ was formed 
by ageing for one day at 465°c followed by two days at 305°c; this gives 
a particle spacing of 3-8 x 10-4 cm (Shaw ef al. 1953). 


§ 3. RESULTS 
3.1. Critical Resolved Shear Stress 


The critical resolved shear stress (CRSS) of a large number of crystals. 
has been determined in this study (Byrne 1960). Values are plotted as 
a function of temperature for most of the aged states and for crystals of 
approximately the same orientation (fig. 3). The variation of CRSS with 
temperature is small for all states above about 150°K. The reverted 
crystals and those containing GPI show a very rapid temperature 
dependence of the flow stress below 150°K. The other aged states produce 
only a small variation in CRSS with temperature relative to the reverted 
and GPI states. Thus the reverted crystals or those containing GP I 
lose a considerable proportion of their strength on heating to room 
temperature. It is generally supposed that the strengthening effect of 
GPII precipitates is greater than that of GPI. This is not true at low 
temperatures, as can be seen from fig, 3, by comparison of the CRSS at 
4:2°K of crystals containing GPI and GPII produced at the same 
temperature of ageing, 
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Fig. 3 
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Critical resolved shear stress of Al-1-7 at. % Cu as function of temperature. 
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In fig. 4 the CRSS at 77 and 273°K are plotted against ageing time at 
130°c. The hardening at early times is due to GPI and at later times 
GPII. The difference in the temperature dependence for crystals with 
mainly GPI and mainly GPII is again evident because the 77 and 273°K 
curves come much closer together at longer times. 
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3.2. Stress—Strain Curves 


Stress-strain curves for GPI and reverted crystals are shown in fig. 5, 
GPII in fig. 6, and 6 and 6’ in fig. 7. Kelly and Chiou (1958) verified 
that the slip system is {111} (110) for the Al-1-8% Cu alloy containing 
GPI. In the case of GPII, slip lines are not observed with the light 
microscope, nor with the electron microscope, except for large strains 
(Carlson and Honeycombe 1954-55, Thomas and Nutting 1957-58). 
This has been verified. In the present work a marked change occurred 
in the orientation of the crystal axes with respect to the tensile axis. 
The change in orientation was always found, within an error of 2°, to be 
consistent with slip on the {111} (110) system having the maximum 
resolved shear stress. The Laue spots show only a small amount of 
asterism. In figs. 5 and 6 the resolved stress may thus be plotted versus 
glide strain. Slip on a single system has been assumed to occur throughout 
the test. 

Fig. 5 
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Resolved shear stress-glide strain curves for reverted crystals and crystals 
containing GP I. 


In the case of 6 and 6’ being the major precipitate, deformation developed 
extensive asterism rather than a change in orientation. Thus slip must be 
occurring on several slip systems which is consistent with the change in 
form of the specimen during deformation. For this reason the nominal 
tensile stress is plotted in fig. 7. 

In the testing at 4-2°K (not shown on the curves in figs. 5 to 7) sudden 
yielding was observed in the ‘pre-plastic’ and plastic ranges giving 
serrations in the stress-strain curves. These were accompanied by audible 
sounds. Serrations in the plastic region have been studied extensively 


Fig. 6 
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Resolved shear stress—glide strain curves for crystals containing GP IP. 
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by Basinski (1957) who attributed them to local heating from plastic 
deformation arising because the specific heat is so small at 4:2°K. In 
accord with Basinski, the magnitude of the drops in load and the spacing 
between them increased with elongation. The ‘pre-plastic’ serrations 
occurred only with CRSS values greater than 4:3 kg/mm? and therefore, 
did not occur with 6. Similar ‘ pre-plastic’ serrations were shown in curves 
published for 2024 alloy by Basinski. The effect is not due to frozen air 
on the sample since the system was well purged with the He gas prior to 
the transfer of liquid He, and no condensation was ever noted inside the 
apparatus. It is believed that all the serrations observed here are due to 
the local heating effect. 

The fractures at 4-2°K were similar to that previously recorded by 
Basinski (1957) in an Al-Mg alloy. They occur by shear through a highly 
localized region and as noted by Basinski are associated with the 
discontinuous yielding phenomenon. Fracture is usually preceded by 
several large discontinuities and results when the specimen can no longer 
arrest the highly localized deformation. The fracture surfaces are almost 
plane. Laue x-ray photographs of the fracture surface show very large 
amounts of asterism. 

In reverted crystals and crystals air-cooled from the solution treating 
temperature and immediately tested, precipitation occurs during deforma- 
tion (Carlsen and Honeycombe 1954-55, Kelly and Chiou 1958). The rate 
of strain hardening is then anomalously high. For this reason stress-strain 
curves for the reverted crystals are not shown in fig. 5 for test temperatures 
above 77°K. 

It is at once clear from figs. 5, 6 and 7, that the stress-strain curves for 
crystals containing GP I and GP II are quite different from those of crystals 
containing 6’ and @ precipitates. The stress-strain curves fall into two 
groups, those of crystals containing GP Land GP I are, at low temperatures, 
similar to those of pure aluminium single crystals (Staubwasser 1959, 
Hosford et al. 1960). It has been customary to define three stages of 
work hardening in face-centred cubic metals: Stage I is the region in the 
stress—strain curve where the stress increases only slightly with deformation, 
Stage II is the linear hardening range, and Stage III is where the curve is 
concave to the strain axis. 

The slopes of Stage II of the low-temperature curves do not depend 
markedly on whether GPI or GPII is present, and their magnitude, 
about G/200, is roughly the same as that found in pure aluminium. The 
total increase in flow stress due to work hardening is about the same in 
the alloys as in pure aluminium at the same temperature. 

When the principal precipitate is 6’ or 6, then the curves of nominal 
tensile stress against elongation bear no similarity to those of an aluminium 
crystal of the same orientation. The initial rate of work hardening is 
extremely rapid for both precipitates, being about @/10. This is similar 
to that of aluminium crystals deforming on a number of slip systems 
(Hosford e¢ al. 1960). After a certain amount of elongation, the curve 
levels off and becomes quite flat; no further work hardening occurs. 
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The complete flattening offis not apparent in Dew-Hughes and Robertson’s 
(1960) curves but the other features are the same for crystals deformed at 

room temperature. From fig. 7 it is seen that the maximum stress reached 

depends markedly on temperature and the total amount of work hardening 

is very large. For crystals containing 6’ the maximum nominal tensile 

pie depends logarithmically upon temperature obeying a relation of the 
orm 


In Omax =A = BT 


where 4=3-72 and B=1-7x 10-3/°K for stresses in units of kg/mm?. 
For the crystals containing @ the dependence of g,,,.. on temperature was 
still greater, particularly at low temperatures. 


Fig. 8 
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Ratio of the flow stress at 77 (a7) to flow stress at 273°K (o4,3) plotted against 
elongation for crystals containing GP I, 6’ and 8. 


3.3. Temperature Transitions During Deformation 


Following the procedure of Cottrell and Stokes (1955) a large number of 
measurements were made of the change in flow stress accompanying a 
change in temperature during the deformation of these single crystals. 
In polycrystals and single crystals of pure metals the ratio of the flow 
stress at a temperature 7' (c) to that at some standard temperature, 
e.g. 273°K (0973), is found to be independent of deformation at large 
deformations. In aluminium containing 1-8 at. °% copper aged to contain 
GPI, Kelly and Chiou (1958) found that o7,/o973 initially decreased with 
deformation and then was roughly independent of further deformation. 
Later Kelly (1959) pointed out that o,,;—o;; was independent of 
deformation. In fig. 8, 077/073 is plotted against elongation, and in 
figs. 9 and 10, o77— 073 is plotted against elongation and 0473, respectively, 
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for representative specimens of the various states of precipitation studied 
in the present work. For the states containing principally GP II, o77/c97¢ 
is roughly independent of elongation, and O77 — Oy73 increases slightly with 
deformation. In the case of @ and 6’, 977/73 aNd 07, — O73 always increase 
with elongation, the changes being less for #6. For 6 and 6’, o is nominal 
tensile stress. 
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Ratio of the flow stress at 7' (a7) to flow stress at 273°K (73) plotted against 
temperature for crystals containing GP I. In curve with triangle points. 
op/Oo73 has been multiplied by G5,5/G7 where G is shear modulus. 
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6 p[O73 aNd (7/0273) (Go73/Gp) are plotted against temperature down to 
4-2°K for GP I precipitate (fig. 11), GP II precipitate (fig. 12), 6’ precipitate 
(fig. 13) and 0 precipitate (fig. 14). When the value varied with deformation, 
average values for each temperature are plotted. For GPI the data of 
Kelly and Chiou (1958) are used for the higher and intermediate 
temperatures since essentially identical results for o7;/973 were obtained 
‘with the present alloy. 


Fig. 13 
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o7[Fo73 ANA o7G'y73/0973G'7 plotted against temperature for crystals containing 0’. 


Fig. 14 
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G77 ANA oplty75/o7367 plotted against temperature for crystals containing @. 


No data are plotted for the reverted condition since alternate cycling 
between 77 and 273°K during deformation results in rapid GPI formation 
(Kelly and Chiou 1958). 
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3.4. Strain-rate Transitions During Deformation 


Information complementary to that obtained from changing the 
temperature of deformation was obtained by changing the strain rate 
during isothermal tensile deformation and measuring the resulting change 
in flow stress. Changes in flow stress (c, — og) corresponding to a strain-rate 
transition between €, and ¢, were measured for both €,<€, and €,> 5. 
Three strain rates were used, 3-3 x 10-5, 1-65 x 10 and 1-65 x 10-3 per sec. 
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In figs. 15 to 19, values of o,, —o, (where a, indicates the flow stress for the 
lowest strain rate) are plotted against the common logarithm of the strain 
rate. For GPII, fig. 16, transitions between only the two higher strain 
rates were performed; the zero point was established by extrapolation. 

In crystals containing GPI, and in these only, yield drop effects were 
found on increasing the speed of deformation at temperatures above 
198°k. The maximum values of Ac were then used. 

As previously noted with crystals containing @ and 6”, slip is not on a 
single system. Therefore to obtain Ac the tensile stress changes were 
multiplied by 0-5, an average value for cos¢cosA. For the other states 
the values refer to resolved shear stresses. 
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The increase in flow stress on increasing the strain rate for crystals containing 
GP II. See caption for fig. 15 for explanation. 


The vertical lines in figs. 15 to 19 represent the spread of values. The 
uncertainty is shown to increase with strain rate because of the way the 
origin of the curve is defined. 

While the strain-rate tests showed considerable scatter, the parameter 


O,— 0; 
log P| €y 

remained roughly constant throughout deformation at a given temperature ; 

a, —0, in the figures is roughly a linear function of log €,. 

The slopes of the lines in figs. 15 to 19 are observed to increase as the 
temperature is raised from 4:2°K to 77°K and then decrease again as the 
temperature of deformation is further raised. The 6’ state is an exception. 
in that the slope increases rapidly on heating above 273°K. 
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The large slope at 77°K of the o,—0, vs. logé, line compared to 
4-2 and 273°K is also observed in pure metals (Wiedersich 1960, private 


communication) and Ag—Al solid solution alloys (Hendrickson and Fine 
1961). 
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The increase in flow stress on increasing the strain rate for specimen 
containing 6’. See caption for fig. 15 for explanation. 
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Fig. 18 
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The increase in flow stress on increasing the strain rate for specimen containing 6. 
See caption for fig. 15 for explanation. 
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The increase in flow stress on increasing the strain rate for reverted crystals. 
See caption for fig. 15 for explanation. 
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§ 4. DiscussIon 
4.1. The Reverted State 


At room temperature the CRSS is low, ~3 kg/mm2, and depends little 
on temperature (fig. 3). Below 100°K it rises rapidly to ~8 kg/mm? at 
4:2°x. This effect appears characteristic of solid solutions whether 
supersaturated or not, see, for example, Seeger (1958). Recent studies of 
Cu base—Al and Ag base—Al alloys (Hendrickson and Fine 1961, Kopenaal 
and Fine 1961) also fit this rule. Al base—Zn alloys have been studied by 
Dash and Fine (1960); the temperature effect as well as the CRSS at 
273°K is much smaller in the Al-Zn than in the Al-Cu alloys, suggesting 
that an important factor may be the mis-match in atomic size between 
solute and solvent. This is much greater for Al-Cu than Al—Zn. 

At low temperatures the stress-strain curve is very similar to that of 
pure Al. Thus, the solute element affects mainly the CRSS and not the 
rate of work hardening provided the crystal is deformed at temperatures 
too low for precipitation or clustering to occur during deformation. 

To account quantitatively for the magnitude of the CRSS and its 
variation with temperature is difficult. Friedel (1957) has proposed a 
theory for the CRSS and its variation with temperature based on the 
interaction of solute atoms with dislocations due to the stress field around 
the solute atoms. This theory gives for the CRSS at 0°K 


o=4Gef fumes OES et eae GLY 


where G is the rigidity modulus, f the volume fraction of solute, and ¢ the 
misfit parameter. The value of ¢f can be obtained from the measured 
variation of lattice parameter with composition. The theory predicts that 
the CRSS should vary with temperature in the manner found experi- 
mentally, as shown in fig. 3, through the effect of thermal assistance in 
overcoming the stress fields. However, the value given for o from eqn. (1) 
is only 1 kg/mm?, nearly an order of magnitude less than the observed 
value, and, hence, fails to account quantitatively for the magnitude of 
the CRSS. 

Since the solid solution is supersaturated and tends to cluster, it seems 
likely that the CRSS, particularly at low temperatures, is due to the 
cutting of small clusters of copper atoms by dislocations. In analogy with 
GP I it is expected that the larger clusters would form thin discs on {100} ~ 
planes and these may be important strengthening sources at all 
temperatures. 


4.2. Stress—Strain Ourves of Crystals Containing Precipitates 


From the stress-strain curves at low temperature, shown in figs. 5, 6 
and 7, it is clear that the presence of GPI and GPII affect mainly the 
magnitude of the critical resolved shear stress and do not greatly alter the 
work-hardening characteristics from those of pure aluminium or aluminium 
containing small amounts of other elements in solid solution. 
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We define the term ‘friction stress’ as the stress to move a single 
dislocation unaffected by the presence of other glide dislocations. The 
presence of either GPI or GPII then increases the strength of the alloy 
by introducing an additional friction stress opposing the motion of 
dislocations in the crystals. This additional friction stress is roughly 
constant throughout deformation. On the other hand, the precipitates 6” 
and 6, while affecting the CRSS less than do GPI and GP II, completely 
alter the work-hardening characteristics of the material. It has been 
shown that the deformation of crystals containing 6’ and @ does not take 
place on a single system, but that a number of glide systems must be 
operative. The single crystals deform, as does a polycrystal with no external 
change of shape, while crystals containing GP I or GP II deform throughout, 
as do single crystals. These facts lead to an immediate qualitative 
explanation of the findings reported in § 3.3 of the effect of temperature 
transitions during deformation on the magnitude of the flow stress. For 
crystals containing GPI and GPII the increase in flow stress during 
deformation is not large compared to the values of the CRSS at room 
temperature. When the temperature dependence of the flow stress is 
measured as a function of strain during deformation, one measures the 
change of flow stress due to the temperature dependence of the friction 
stress plus that due to the temperature-dependent part of the increase in 
flow stress due to deformation. In the case of GPI the temperature 
dependence of the friction stress due to zones is large and masks the effect 
of the dislocations introduced during deformation; thus o,,— 0,3 changes 
little during deformation. 

For GPII the temperature dependence of the friction stress is much 
smaller and so the effect of the increase in dislocation density during 
deformation can be detected as an increase in o7,— 75 (see figs. 9 and 10) 
during deformation. Some 6’ is always present particularly with GP II 
produced at 190°C; this will lead to an increase in o77— 0973. 

The differences in behaviour among the crystals containing GP I, 
GPII, 6’ and @ are strikingly seen if one plots the quantity o,,—o 73 (the 
temperature-dependent part of the flow stress) not against strain but 
against the flow stress at 273°K. This is done in fig. 10. If the increase 
in the temperature-dependent part of the flow stress during deformation 
is taken as a measure of the increase in the number of ‘forest’ dislocations 
which must be cut by gliding dislocations, as it is in all current theories 
of work hardening, then we see from fig. 10 that the values for GP I and 
GP II are much the same as for pure aluminium. This fact, in conjunction 
with the similarity of the stress-strain curves, is convincing evidence that 
the essential mechanism of work hardening in pure aluminium and in 
crystals containing GPI and GPII is the same. This is also true for 
aluminium-silver alloys containing zones and thus at least in these 
aluminium alloys it can be said that zone-type precipitates alter the plastic 
properties mainly through the introduction of an additional friction stress. 

Crystals containing 6’ or 6 deform on a number of slip systems. This 
leads to an increase in the temperature dependence of the flow stress 
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throughout deformation (figs. 8, 9 and 10) due to the ever greater number 
of dislocations which must be cut by a dislocation moving on any operative 
glide system. As shown in fig. 10 for 6’, the increase in O77 — Oo75 attains 
a value of more than 4 kg/mm?, indicating that the increase in density of 
the ‘forest’ is an order of magnitude greater than in pure aluminium or 
in crystals containing GPI or GPITI alone. 

Since the temperature dependence of the flow stress in crystals 
containing GPI and GPII is almost completely controlled by the 
temperature dependence of the CRSS, this means that tests involving 
change of strain rate during deformation can be used to investigate the 
strain-rate dependence of the CRSS, whereas in 6’ and @ containing 
crystals the strain-rate dependence of the CRSS cannot be simply 
investigated by changes of strain rate during deformation due to the large 
change introduced by the deformation itself. Further, several activated 
processes are involved; the o,/o,,; curves turn downward at higher 
temperatures and with 0’, do/dé is large for 373°K. 


4.3. Values of the CRSS of GP I and GP II 


In the case of crystals containing GPI and GPII the separation of the 
centres of precipitates cutting any slip plane, deduced from the x-ray 
results and confirmed by electron microscopy, is less than 1504 and less 
than 300A respectively. The stress required to bend a dislocation to a 
radius of curvature pis o~Gb/p. To avoid precipitates, dislocations would 
have to be bent to a radius of curvature of about one-half the mean 
separation of precipitates. This requires, in these cases, stresses of the 
order of magnitude of the theoretical shear strength, and an order of 
magnitude greater than the observed flow stress. The dislocations are 
thus forced to cut the zones (Kelly and Fine 1957). This has been confirmed 
experimentally in GP II containing crystals (Nicholson e¢ al. 1960) and in 
other aluminium alloys containing coherent precipitates of the spacing of 
GPI (Jan 1955, Livingston and Becker 1958, Sato and Kelly 1961). 

Comparison of the CRSS data in fig. 3 for the structures containing 
GPI and GPII clearly shows that increasing the size of the obstacles, 
which must be cut, reduces the temperature dependence of the CRSS. 

The thermally activated motion of dislocations through obstacles in _ 
their glide plane has been studied by a number of workers. Assuming 
the barrier is cut in a single activated event and following the method due 
to Seeger (1954) and to Mott (1956), we write the strain rate 


é=€,exp [— U(a)/kT] ia eel as ae) 


where é, is a constant, & and 7 have their usual meanings and U(c) is the 
activation energy for the cutting of a single obstacle under the influence 
of a shear stress c. Solving (2) for the shear stress gives the CRSS as a 
function of temperature. Mott has shown that for a barrier with a rigid 
energy profile the dependence of U on a is given by 


U(c) = U4(1 — 2/05)?” a RS ae -h io) 


P.M. 4D 
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for the case when the effect of thermal fluctuations is very small, and 
hence o is always approximately equal to oo, the stress required to shear 
an obstacle in the absence of thermal fluctuations : 


feeeie 
Co 401,Xo . 
Here U, is the increase in potential energy of the system if an obstacle 
is cut at 0°K in the absence of an applied stress. J, is the separation of the 
obstacles in the glide plane and X, depends on the energy—distance profile 
of the barrier. 2X, represents the distance moved by a dislocation into 
the barrier during the process of cutting. 

Equation (3) should be used if the effect of thermal fluctuations is small. 
In the case when the barrier is thin and the effect of thermal fluctuations 
is large, we can use the approximation (Cottrell 1953, Seeger 1954) 


U =U, — 401, X oo. 5, ere 
Substituting either (3) or (4) in (2) and defining the activation volume 
v=4b1,X, one obtains 
U, CE arn ae 
o= —8/ 1 7, ela} ] 3 (5) 
or 
U. fe 
a= <4] 1- FIn(ae) | sib irgimtod: “eg anton 


By differentiating (5) and (6) to obtain expressions for the strain-rate 
dependence of o one obtains 
do (EEA UN 


= —————_—_. . . ww. a . (7 
diné  3/2v(In €)/€)*8 (7) 


and 
do af kT 


diné v~ 


(8) 


When a dislocation cuts a zone it must overcome any misfit stresses 
associated with the zone, and it must change the arrangement of the atoms. 
Both of these must be evaluated and taken into account in any quantitative 
theory of zone hardening and its variation with temperature. 

The zones in Al—Cu present two types of barrier. One-third of the zones 
lie parallel to the Burgers vector of any moving dislocation. Consideration 
of the structures of the zones shows that for GPI and GP II a dislocation 
passing through such a zone will only alter the arrangement of nearest 
neighbours on entering and on leaving the zone. The other zones present 
a much larger energy barrier and when these are cut the type of nearest 
neighbour is altered over the whole area of intersection of the zone and 
slip plane. 

Equation (5) predicts a linear relation between the CRSS and 723; 
eqn. (6) predicts a linear relation between the CRSS and 7. In fig. 20 
the CRSS is plotted versus 77. A straight line has been drawn through 
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the GPII data for 130°C; the fit appears reasonably good. For GPII 
(190°c) containing an appreciable amount of 6’ there is positive curvature 
at the higher temperatures. In the case of GPI the three lowest 
temperature points appear to fit on a straight line. 


Fig. 20 
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Critical resolved shear stress plotted against 7?’ for 
crystals containing GP I and GP II. 


Considering the data for GPI we assume that there are two kinds of 
barrier, one with a large U, and one with a small U,. Then in fig. 20 two 
intersecting straight lines have been drawn through the data. We assume 
eqn. (5) to hold for both, but with different values of Uy. For the low 
temperature branch we see from eqn. (5) that at a temperature 7’, where 

U ,=kT ln &/€ 
the flow stress due to this barrier becomes zero. The critical temperature . 
is about 150°K, see fig. 3, whence, evaluating the log term, see below, 
U,~ 0-26 ev. 

The high temperature branch in fig. 20 gives a value of the flow stress 

at 0°K of 9kg/mm?. Whence from eqn. (5) 
Coes 9 kg/mm? 
v 
and from the slope of the line 
(OF 
v 


(ln é,/€)28=4 x 10-2 kg/mm? (°K)? 
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We can estimate €,: 

é,=NAbry 
where N is the number of points per unit volume at which dislocations 
are held up at obstacles. A is the area of slip plane swept out by a 
dislocation per activated event, b the Burgers vector and v) the vibration 
frequency of a dislocation against an obstacle. We take 


Vqg= 0 seer. 


A may be calculated from the known distribution of the zones. In this 
case it is 5 x 10-14 cm?. We assume a value of NV of 5 x 1018 corresponding 
to a dislocation density of 108 cm~?, which gives a value of the In term 
of 22. We then find U,~6ev and v~ 10-* cm?. 

For the low-temperature branch the slope is 18-6 x 10-* kg/mm? (°K)? 
and the appropriate value of the flow stress is obtained by subtracting 
9 kg/mm? from the 0°K intercept in fig. 3 to give 


oes 4-9 ko/mm?, 
v 


Taking the same value of the In term we then find U,~0-2ev and 
v~7x10-22cem3. The value of U, agrees roughly with that estimated 
from the critical temperature. 

Additional estimates of v may be obtained from plots of Ac against 
Alogé (fig. 15) using eqn. (7), since we can estimate U, from the 
temperature dependence of the flow stress and the term involving (log é)1/? 
varies negligibly over the range of strain rates used. Taking U,=0-2ev 
for temperatures of 4-2°K and 77°K, but 6 ev for 273°K gives values of v of 
4x 10-*2, 5 x 10-*% and 2 x 10-?° cm? for these temperatures respectively. 
The agreement with the values found from the temperature dependence 
of the flow stress is good, being always within a factor of two, and lends 
great support to the correctness of the analysis, when it is borne in mind 
that only average values can be obtained for what must in reality correspond 
to a range of values of both U, and v. 

We now consider the data for GP II. The CRSS varies with temperature 
approximately as 7" (fig. 20), and this is consistent with eqn. (5) assuming 
the occurrence of a single activated process. However, the strain-rate data 
(fig. 16) are not consistent with this analysis since the strain-rate 
sensitivity does not increase continuously with rise in temperature. We 
shall assume that there are again two types of obstacle to flow due to the 
different orientation of the zones, as for GP I. We assume that the obstacle 
responsible for the increase in flow stress as the temperature is lowered 
below room temperature becomes ineffective at around room temperature, 
and thus from eqn. (5) corresponds to an activation energy of about 0-6 ev, 
taking the value of the log term as 25. The increase in flow stress due to 
this obstacle is then the difference in flow stress between room temperature 
and 0°K, i.e. for GPIT produced at 130°c, 2:5kg/mm?. We then have 


= = 2-5 kg/mm? 
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and from the slope, which is 7-2 x 10-2 kg/mm? (°K)?3, we obtain. 
U,=0-5 ev and v=3x 10-*! em’. The value of the log term is 25 (using 
A =2-87 x 10-"! cm? and a dislocation density of 108 em-). 

We can also compute values of v from the strain-rate data as we did 
for the GP I containing crystals. We then obtain, taking U,=0-5ev and 
using eqn. (7), v=2 x 10-7! em? at 4:2°K and 3 x 10-21em3 at 77°K. These 
values show good agreement with one another and with the value found 
from the temperature dependence of the flow stress. If we use U,=0-5 ev 
and evaluate the strain-rate data at 273°K we obtain a value of 
v=2x10-°°cm%, which does not agree with the values obtained at the 
other two temperatures. It is for this reason that we suggest that there 
are two quite different types of barrier corresponding to the differently 
oriented zones in the case of GP II as in the case of GP I. 

We now consider the values obtained for activation energies and 
activation volumes, using the above analysis. For GPI the value of U, 
of 0-2 ev we associate with the cutting of zones where b is parallel to the 
plane of the zone and 6ev with cutting of zones where b is oblique to the 
plane of the zone; the U, for oblique cutting is 30 times the value for 
parallel cutting. The length of a zone is about 100A and this gives a value 
of about 30 atoms for the average trace of a zone intersecting a {111} 
plane. The agreement between the two estimates is thus very good. 

For GP II the results indicate that the U, for b parallel to the plane of 
the zone is about 0-5ev; it is thicker than GPI. For b oblique to the 
plane of the zone U, appears to be so high that the contribution to the 
CRSS is little affected by variation in temperature. 

The consistency of the above analysis lends great support to the idea 
that the major contribution to the flow stress of these alloys, containing 
GPI or GPII is due to the work done by the applied stress in forcing 
dislocations through the zones. The values of v, the activation volume, 
are also reasonable. A crude estimate of the activation volume for zones 
parallel to the Burgers vector is given by v=, where /, is the average 
separation of the zones and 6 the distance of closest approach of atoms. 
For GPI produced at 130°c, J,~ 1004 (Nicholson and Nutting 1958) and 
thus v~7 x 10-22em3. For GPII, J, ~ 3004 (ibid.) and so v ~ 2 x 10-7 cm’. 
These values agree with those deduced from the temperature dependence 
of the flow stress and independently from the change of strain rate during . 
deformation. 

For the zones which do not lie parallel to the Burgers vector of a 
dislocation the appropriate value of v will be v~b1,”, which gives a value 
of 2 x 10-2° cm3, again in agreement with the value found from the results 
at room temperature. 

We may now consider the value of Uo, the energy required to shear a. 
zone. As stated previously, this will include (a) any stress fields around 
and in the zone which must be overcome and (5) the energy required to 
change the arrangement of the atoms. 

A crude estimate of that part of Uy required to change the arrangement 
of atoms can be made from the heat of zone formation or the heat of solution 
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(reversion) of the zone (Kelly and Fine 1957). For a GPI zone whose 
trace on the slip plane is parallel to b of the dislocation involved, Uo is 
roughly Ae/2 where Ae is the heat of formation of the zone per atom of 
copper. Three values of Ae have been reported. These are 0-08, 0-23 
and 0:14 ev per atom of copper. The first, Suzuki (1949), and second, 
Beaton and Rollason (1957), are heats of solution estimated respectively 
from differential thermal analysis on heating a sample containing GPI 
and from the solvus line with respect to GPI determined from hardness 
data. The third and most reliable (De Sorbo ef al. 1958) is an isothermal 
calorimetric measurement of the heat evolved during formation of GP I. 
Thus a U, of about 0-lev is predicted or about half the value indicated 
from the measurements reported here. The estimate of U, from the 
calorimetric measurements is a lower limit since it is assumed that all 
copper atoms are in the zones. 

For a GPII zone about 244 thick and with b of the slip dislocation 
contained in the face of the zone, U, is roughly 2Ae where Ae is the heat 
of formation of GPII per atom of copper. Two values of Ae for GP II 
have been reported, 0:3lev (Beaton and Rollason 1957) and 0-14ev 
(Suzuki 1949). These are rougly 1-5 times the corresponding values for 
GPI; GPII is a more stable phase. Assuming Ae to be 0-2ev then 
a U, of 0-4 ev is predicted or again somewhat less than the measured value. 

Additional contributions will come from the strain fields of the zones 
which are known to exist from x-ray examination (Gerold 1958) and from 
direct observation under the electron microscope (Nicholson and Nutting 
1958). There will be a short-range elastic interaction between zones and 
a moving dislocation exactly analogous to the elastic interaction between 
a moving dislocation and another stationary one threading the glide 
plane. The magnitude of the elastic interaction is difficult to estimate 
without a very detailed knowledge of the stress field around and in the 
zones, but the low values of U, predicted from the change in arrangement 
of atoms suggests that the contribution from the elastic interaction may 
be of about the same order of magnitude. 

Kelly (1958) and Dew-Hughes and Robertson (1960) attempted to 
account for the CRSS of these alloys in terms of the theory of Mott and 
Nabarro (1940). Dew-Hughes and Robertson claim good agreement 
with experiment at room temperature, by adding together various 
contributions to the flow stress. They conclude that at room temperature 
the CRSS for GP I containing crystals is satisfactorily accounted for from 
the theory of Mott and Nabarro, plus a small contribution from the solid 
solution. The theory of Mott and Nabarro cannot apply in this case without 
considerable modification. Firstly, the formula derived is for spherical 
precipitates and secondly the flow stress is identified with the mean 
internal stress. This is the only case when the radius of curvature to which 
a dislocation can be bent by the mean internal stress is equal to the 
separation of the precipitates. We have shown above that the measured 
flow stress, which is equal to the mean internal stress in this theory, is far 
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too small to bend a dislocation to a radius of curvature equal to the 
separation of those precipitates cutting the glide plane. The dislocations 
actually cut through the precipitates in the glide plane and hence each 
dislocation is subject not to the mean internal stress but to the stress due 
to a zone at small distances from the zone. Finally, the theory of Mott 
and Nabarro predicts a temperature dependence of the flow stress equal 
to that of the elastic constants. This is not the case. 

The inapplicability of the theory of Mott and Nabarro in its current form 
does not mean that the stresses around a zone are unimportant, but that 
it is the short-range elastic interaction between dislocations and zones 
which are important here and not the mean internal stress. 


4.4. Alloys Containing 6’ and 6 Precipitates 


The CRSS of alloys containing the 6 and 6’ precipitates appears to be 
understood at room temperature. The precipitates do not deform with the 
matrix and the CRSS is the stress necessary to expand a loop of dislocation 
between them. This is given by (Orowan 1948): 


o=ab/L ae St Bi een aoe CS) 


where L is the separation of the precipitates. « is the relevant constant 
in the theoretical expression for the line tension of a dislocation. This 
formula applies well over a range of spacings (Dew-Hughes and Robertson 
1960). In our case for the 6 precipitate L is 3-8 x 10-4 cm, hence o, with 
a=4, using Cottrell’s (1953) estimate of the line tension, is 0-8 kg/mm?. 
The measured value is 1-2kg/mm?. The CRSS at 4:2°K is 1-8 kg/mm?. 
The change can be attributed to the increase in G with temperature, 
the increase in strength of the equilibrium solid solution and perhaps a 
small increase due to the density of ‘forest’ dislocations. 

For alloys containing 6’ the same appears to be true. The separation 
of the precipitates is 0-4x10-4cm (Nicholson e al. 1958-59), hence 
o=5kg/mm?; the measured value is 4-3 kg/mm?. On cooling to 4:2°K 
there is only a small increase in CRSS of 1-7kg/mm?. This may be 
attributed to an increase in G (0-5 kg/mm?) and to the presence of 
dislocations cutting the glide plane. These may be either distributed 
at random or associated with the precipitate. Two types of dislocations. 
associated with 6’ have been observed (Nicholson et al. 1958-59): a 
Van der Merwe grid in the interface and a ring around the platelet and 
parallel to its plane. Thus, 6’ is not truly coherent with the matrix; 
much of the misfit is adjusted by the presence of dislocations. Hardening 
due to misfit strains does not appear to be important here. 

With regard to the stress-strain curves (fig. 7) for both @ and 6” the 
initial work-hardening rate is very large being about G/10 and roughly 
independent of temperature. The amounts of work hardening reach 
maximum values which are many times larger than the yield stresses and 


increase considerably on cooling. 
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At small strains 6 particles are not deformed. This has been shown 
experimentally (Dew-Hughes and Robertson 1960); the precipitates are 
large enough to withstand the stresses on them as the dislocations expand 
between the particles. Following Fisher et al. (1953), Dew-Hughes and 
Robertson (1960) concluded that the initial work hardening is due to the 
residual loops remaining around the particles. These impose a back stress 
on the source. The observed temperature dependence of the initial slope 
is further confirmation since the back stress would only vary as G. 

According to the theory of Fisher, Hart and Pry the metal work hardens 
until the particles are sheared. A plateau in the shear stress—shear strain 
curve is predicted. Hirsch (1957) suggests that the particles are avoided 
by cross slip. Ifso prismatic loops of dislocation will be formed at particles 
and these will give the same mean internal stress calculated by Fisher, 
Hart and Pry, but a reduced stress on the particle. 

At large stresses at room temperature, cutting of some particles is 
observed (Thomas and Nutting 1957-58, Koda and Takeyama 1957-58, 
Nicholson e al. 1960). @ particles. are not cut at room temperature. 
Whether the particles are cut or whether cross slip occurs will depend upon 
the energies involved, and on the size, spacing, and strength of the 
precipitate. With 6’ there is an orientation variable too, since the 
precipitate and matrix lattices are not in registry in all directions. 

Plastic deformation of specimens containing @ and @’ do not show many 
slip limes (Thomas and Nutting 1957-58) indicating that slip is taking 
place on a fine scale. Further, Laue spots develop asterism rather than 
change position. Simultaneous slip on several systems is indicated; the 
dislocations in part appear to be undergoing cross slip to avoid the particles 
as suggested by Hirsch and the occurrence of the maximum in the 
stress-strain curve could correspond to complete replacement of the 
Orowan mechanism by cross slip. Since cross slip in part is aided by 
thermal activation a much higher maximum at 4:2°K could result. One 
thus expects to observe more extensive shearing of particles at low 
temperatures, and the maximum stress may then be reached when the 
particles are sheared. 

Further work is necessary to completely understand with @ and @’ the 
very rapid rate of work hardening, the occurrence of the plateau, and its 
variation with temperature. These are certainly connected with slip on 


many systems, and a rapidly increasing dislocation density with 
deformation. 
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A Dislocation at a Free Surface 
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ABSTRACT 


The stress field of a dislocation meeting a free surface of an elastic body is 
calculated for any angle of incidence and any Burgers vector. 


§ 1. IyTRopucTIOoN 


Many observations have been made on dislocations lying in thin films 
of crystalline material. It is to be expected that the presence of the two 
surfaces of the film would have an effect on the behaviour of the dislocations, 
and this was confirmed by Eshelby and Stroh (1951) for dislocations normal 
to the surface, and by Head (1953) for dislocations parallel to the surface. 
The general case of dislocations in a semi-infinite elastic medium has been 
studied by Steketee (1958). His solution involves six sets of Green’s 
functions and is not easy to apply in practice. In the present paper a 
straight dislocation meets a free surface of an isotropic elastic body at an 
arbitrary angle, and the resulting stress field is to be calculated. 


§ 2. NoTaTIon 


Rectangular cartesian coordinates (x, y, z) are chosen, and the elastic 
body is supposed to occupy the half space z<0, as shown in fig. 1. The 
dislocation lies in the (y, z) plane at an arbitrary angle «to the inward normal, 
the origin lies at the point where it meets the surface and the x axis is taken 
upwards out of the plane of the paper. 

The stress field of an infinite straight dislocation so situated would con- 
tain stress components %2, 72 and 22 on the planez=0. Since this plane is 
to be stress-free it is necessary to find the ‘image’ field which satisfies 
the elastic equations, compensates the stresses on the surface and introduces 
no new singularities in the region z< 0. It is probable that singularities will 
be introduced in the region z> 0. 

The Burgers vector is taken in turn in the x, y and z directions with 
b,, 6, and b, denoting the magnitude of the vector in each case. A separate 
solution is given for each direction, and the general case, with Burgers 
vector in any direction, may be solved by a linear combination of these 


solutions. 
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Fig sl 


The harmonic function used in determining the‘ image’ field is denoted 


by ¢. 
§ 3. DistocatiIon NORMAL TO THE SURFACE 
If the dislocation meets the surface normally the angle « is zero, and at 


least one of the boundary conditions is automatically satisfied. 
The screw dislocation gives 22 zero but shear stresses on the surface, 


gue tC aeoed 
2r x+y? 
Tipe 


This case was solved by Eshelby and Stroh (1951) who found the ‘image’ 


system : 
nies bas ot GS Bae 
7 PZ Qn r—z” 
(1) 
27 r(r—z) 27 r(r—z) 


where 7? =a? + y? + 2? and yw is the shear modulus. 
The edge dislocation has %2 and ¥2 zero with one remaining surface stress 


b 
lee Bi coh. a enmeeae tay 
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where o is Poisson’s ratio. This case is solved by finding the harmonic 
function » such that 


Op opb y 

=—_—- FS Fer ,—?/. h =v. 

a2 a(lte) akan When! 20 
It is related to the harmonic functions ¢,, ¢,, ¢, and ¢3 of Neuber (1937) 
by the equations 


and is described by Green and Zerna (1954). The solution in this case is 


_ op, Frise, Sih Yes ’ 
p= ae ( ylog(r z) =), Le aera a es.) 


The displacements of the ‘image’ system are then derived from q, i.e. 


> opb., Qxyz i (1— =a ) 
2a(1—«a) \ r(r—z)? (r—2z)? j? 
0p 


See ee OG) 
i <oma 2 


iS sate ((1-20)log(r—2)— (8-20) esa) 


27(1l—o 
(3—2a)y? 2? ) 


ee Sos 
For z= 0 this result is in agreement with Indenborn (1959). It may be seen 
that the normal stress in this case is : 


ae “He ( y -#) 


Sek a(l—o)\r(r—z) 
which cancels that of eqn. (2) at z=0. The image shear stress xy has the — 


surface value 


— _ (l—2c)oub, (% _ ) 
= See ag rn ee 


which, for o~ + is equal to about 10% of the corresponding stress for a 
dislocation in an infinite medium. 
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§ 4, OpLiguE DisLocaTIon, 6, 

To find the stress field of a dislocation meeting the free surface obliquely 
it is convenient to use one of the angular dislocations described in a previous 
paper by the author (Yoffe 1960). This angular dislocation is composed 
of the original dislocation of fig. 1 and its image in the surface plane, as 
shown in fig. 2. The stress field of such a dislocation in an infinite 


Fig. 2 


elastic body is known from the previous paper and the corresponding 
displacements may be written down immediately : 


2uu= pbs (tons 2 _tan-171 +tan-! wr sin 2c ) 
27 c x 


yn’ —x* cos 2a 


ee Bee (. nag OM 
Sat Lea) \r lek) earn 
Sup Mee fy sine Sy Bin ary 
) 


(4) 


i ae , — (1-20) cos «log —D(r-L) 


2000 = 
CS 4e(L 0) 


(= 0) oa ee ert) 
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aeprrnnsd ap oa 2c) sin «log a) 
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where n=Yy cosa—zsina, 


C=ysina+zcosa, 
n’ = —ycosa—zsina, 
¢’=ysina—zcosa. 
It is found that the shear stresses £2 and 72 derived from eqns. (4) vanish 
on the plane z= 0, while 22 has the surface value : 


—  pb,cosa 2o sin « Yy COS? (1—20)y 
2a(1—o)\(r-ysina) (r—ysina)? r(r— ysina) 

To remove this stress it is only necessary to find a suitable harmonic 

function y such that 


2 
ze =22 on z=0. 
2 
Such a function is found to be: 
_ pb, cota 


YP alieec) { Clog (r — £) + 2(1 —o)r + (1 — 20)z log (r —z) 


+2(1—c) cota | log (r—2)—ylog (r—2) 


eeeteran tan) tan : (6) 
x c X* CoS a+ YN 


A stress field is readily derived from y as in §3. Not all components will 
be written in full for reasons of space: 


_ Op = 0p 

HEE "ere + (1-20) an 
_ pb, cota RSE lg POP eta egg be pee 
= Pecote sa plies . ae eR ser 


XZ COS a uz LZ COS a 


rr) rr—-O8 ro 


Molt a. We se ae a1 __ersina 
2(1—a)(1 2a) oot a ( tan 5 tan =v an cacaaiie 


—2(1—o) 


2m 25h + (1-2 eo 
guwa2 et — 21-0) 2, 
ay eae, +12) oe maz SE, 


(7) 


The complete stress field of the dislocation meeting the free surface is 
obtained from eqns. (4) and (7). It may be verified that the function 
of eqn. (6) reduces to that of (3) as a> 0. 
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§ 5. OBLIQUE DisLocaTION, b, 

The dislocation with Burgers vector 6, in the y direction may be treated 
in a similar manner. The angular dislocation again leaves only 22 to be 
cancelled on the surface: 

2 1-2 
22= Fae OO #583 ae = icin Rite daechee ac 
2n(1—«) (r—ysina)® r(r—ysin a) 
The function ¢ in this case is : 
eye Oe 2(1—o)x cos alog (r—z 
S| red 0 ) 
—(1— 20+ cos? «)x log (r—Z) 


Tee 1-2 tan 4 —tan-17 taht (9 
(ycosa-+(1=20)n) (and 4 aaa te 


Oo 


which reduces to 


when a->0 for normal incidence. 
The complete displacements are found by adding those of the angular 


dislocation to the additional ‘image’ terms from eqn. (9) : 


+ 


b 42 2 
Dine ad ipa ae et i eee 
ion de Lesoyy = ( he =e c) ) 
Oy Oy 
hes es es | 
beoj- ores ee nee 
_ pb ut ar sin 2x 
B) = ited G Sb fay =e =f 
[Lv maa an . tan = +tan aes 5) | 


pby c— sinacoSsa  YyCOSa em) 


(toeo) mae r—€ or(r—Q) rr) 
+2 Le + (1-20) op 
Oy dz oy” 
a pb, (—F- COs?a  zC0Sa  — z COSa 
dor( Vea) on 0) Poe res eee) 
Os op | 
+27 2(1— ew 


(10) 


§ 6. OBLIQUE DisLocaTion, b, 


The angular dislocation of fig. 2 with Burgers vector b, does not give zero 
shear stresses on the surface : 


Jeo as ( x sin a »  osina. , (l—2a)y 
2r(1—o)\ r(r—ysine)? r—ysina Gas) ; 
als ares y sin sin? x (1—o) (11) 
2n(1—oc) \r(r—ysina)® (r—ysina)®  r(r—ysina) ) j 


zz=0, 
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on z= 0. These stresses could be compensated by an ‘image’ system 
derived from another harmonic function, # say, if can be found such that 


fail 


Merman aaa ZN oe 


However, in this case a more direct method is to arrange two angular 
dislocations as in fig. 3, thereby reversing the sign of theimage. The surface 
stresses are then given by 


le el OY 
Se een 
aw bbe # 12 
Ye Qo | 7?” 7 
os pb, ; 1 -  _GOs* a ) 
= — 4 — .xsin al ——__—- + ——————__,5 

2n(1—c) r(r—ysina) (r—ysina)® 


The shear stresses are simply those of the normal screw dislocation in § 3 
and may be cancelled by Eshelby and Stroh’s solution, The direct stress 


P.M. 4E 
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Zzis removed by the same procedure as before. The function ¢ is found to be 


pb, 


y 2n(1 0) 


{ -#sinalog (r—C)+z 


Sa Mare i: ae (13) 
x (tan Z tan pont an ear ) 


In this case » itself vanishes as «> 0. 
The complete solution is composed of terms derived from (13), from the 
two angular dislocations, and from Eshelby and Stroh’s solution (1): 


yy — Hossine (1=26)log(r— 0) (r= 0) - > ) 

PT An (Ta) (r—t) arr au) 
py Pp 4 ap 
a(r—2) ah 0xdz mu 2) Ox’ 
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2yno= Fe (tan 5s tan ra a a) 
pb,esina/cosa Cosa 2  — 2 ) 
4n(l—o)\r—-€ r—-f’ wr(r—Z) r(r—Z’) 
ep dy 


§ 7. CONCLUSION 


The expressions for the stress field of a dislocation meeting a free surface 
at any angle show that the ‘image’ system is quite simple in form. It 
consists of the image of the original dislocation, and a few additional 
elastic singularities distributed either on the image or on the outward 
normal from the point of incidence. It should not be impossible therefore 
to repeat the process and attempt to estimate the effect of a second free 
surface parallel to the first. This question requires further work and will 
not be answered here. 

From casual inspection it seems unlikely that the dislocation would 
remain straight in all cases. In particular the straight form may sometimes 
be in unstable equilibrium with a tendency to bow out in the slip plane, in an 


attempt to lie along its image. The forces acting on the dislocation near 
the free surface must be calculated in each case. 
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Erratum. The author wishes to correct a misprint in the previous paper 
(Yoffe 1960). 


The second of eqns. (3) should be: 
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ABSTRACT 


Calculations have been made of activation energies for vacancy pair 
movements in NaCl and KCl using the approximation of a static lattice. 
The numerical results indicate that these energies are comparable with 
or greater than the activation energies for single vacancy movement. 
Qualitative reasons are given showing why this should be go. 


Tue classical Born theory of ionic crystals as developed by Mott and 
Littleton (1938) for lattice imperfections has been used recently by Guecione 
et al. (1959) to discuss the activation energies for ionic movements into 
single vacancies in NaCl and KCl. It was found that observed activation 
energies for single vacancy movement are accurately reproduced by the 
theory provided that the closed-shell interactions between ions are rep- 
resented by the ‘Born—Mayer—Verwey ’ potential. This potential, of the 
form A + Br-?? for interionic separations r less than the equilibrium anion— 
cation separation, had previously been proposed by Verwey (1946); it had 
been shown earlier by Verwey and de Boer (1940) to lead to a good account 
of the interionic distances in free molecules. Calculated activation 
energies were not in agreement with experiment when the ‘ Born—Mayer’ 
repulsion potential was assumed, although this potential is well known to 
account for the low pressure compressibilities of ionic crystals. This is 
not surprising since the distances of approach of ions in the activated state 
are close to those in the molecule and much less than the normal interionic 
spacing in the crystal (figure). 

It is evident that a substantial part of the activation energy for single 
vacancy movement is made up of the closed-shell repulsion between the 
moving ion and its two close neighbours in the activated state (ions a, b in 
the figure (A)). It has therefore often been assumed that the activation. 
energies for ionic movements into a vacancy pair should be less than the 
corresponding quantities for single vacancies (figure (C)). This idea has 
been supported by a direct calculation for the Cl- jump into a vacancy pair 
in KCl by Dienes (1948) who obtained the value 0-38ev. The concept of 
mobile vacancy pairs has found favour in the rationalization of the processes 
of colour centre aggregation and the annealing of radiation damage (Seitz 
1946, 1954). However, their direct appearance in diffusion studies has 


OO 
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remained unconfirmed until recently (Schamp and Katz 1954, Laurance 
1960). Laurance studied the tracer diffusion of Cl- in NaCl containing 
variable additions of CaCl, and did find a contribution which could be 
assigned to vacancy pairs. ‘This contribution was governed by a total 
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A. The activated state for a negative ion jumping into a single vacancy. B. 
The ground state of the vacancy pair, showing assumed relaxation of 
neighbours. C. The activated state showing the assumed relaxation of 
neighbours. 


activation energy of 2-49ev. Similar inferences have been made by 
Barr and Morrison (1960) from results obtained with NaCl containing 
CdCl, additions. Now, theoretical estimates} of the energy of formation 


+ Tosi and Fumi (1958). An earlier calculation for NaCl was made by 
Reitz and Gammel and a more recent calculation has been made ona different 
basis by Johnson and Parker (1960). 
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of vacancy pairs, E, using the Born—Mayer potential yield the values 
1-3-l-4ev. Hence Laurance inferred that the energy of movement is 
about 1-1-1-2ey, i.e. as large as for the Cl- jump into the single vacancy. 

On account of this conflict with the intuitive idea that vacancy pairs are 
highly mobile we have re-examined this question theoretically. Firstly 
we noted that the previous calculation by Dienes used the Born—Mayer 
form for the closed-shell repulsions and that this form was too soft to 


Activation energies AH for ionic movements into vacancy pairs compared with 
the values for single vacancies computed by Guccione et al. (in brackets). 
The calculation of AZ has been performed at intervals of é of 0-02 and 
the components of AZ refer to the value of é closest to the minimum; the 
values of these components are therefore approximate since although the 
minimum in A# is shallow the components vary more rapidly with €. 
For the first four entries the Born—Mayer—Verwey potential has been 
used to represent closed-shell interactions. The corresponding values 
for the energies of formation of vacancy pairs are Hy=1-27 ev (NaCl) 
and 1-28 ev (KCl). For the fifth entry the Born—Mayer form has been 
used with constants proposed recently by Tosi and Fumi: this potential 
is less hard than the Born—Mayer—Verwey potential and gives an energy of 
formation #y=1-24ev. For the sixth entry the Born—Mayer potential 
used by Dienes (1948) has been employed; this is the softest of the three 
potentials and gives H;=1:19 ev. We believe the results predicted by 

the Born—Mayer—Verwey model to be the more reliable. 


Moving 


oe AE (ev) | AE (ev) | AE (ev) | AZ (ev) | Dipole 


correction 


Nat in 


NaCl 1-46 6 — 0-02 0-02 0-00 
. 2 


CF in 
NaCl 1-27 1-67 0-24 — (0-87 0-23 
‘ll 1- ).¢ 


K* in 
KCl 1-30 


Cr in 
KCl 
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account for the observed activation energies for single vacancy movement. 
Initially therefore we carried through the calculations using the harder 
Born—Mayer—Verwey potential, but for comparison we have repeated them 
using the Born—Mayer form. The activation energies have been calculated 
by finding the difference in energy of the activated configuration (figure 
(C)) and the ground configuration (figure (B)). All 10 neighbours of the 
two vacancies in the ground configuration have been allowed to relax out- 
wards by an amount A which has been determined by minimizing the energy 
of the crystal. Similarly the neighbours of the three vacancies present in 
the activated complex have been relaxed outwards by an amount 7, 
except for the ion at (010) which has been relaxed by an amount +/2¢ in 
the direction (—1, +1, 0), ie. to the point (—¢, 1+@, 0). The position 
of the moving ion in the activated complex has been taken as ($+ €, $— &€, 0). 
The energy of the crystal in this activated configuration has been minimized 
with respect to 7, ¢ for a range of values of €. The value of € requiring the 
lowest activation energy was then found; the minimum so obtained was 
rather shallow. As in previous calculations of this kind the energy is 
regarded as being composed of Coulomb interactions (monopole—monopole), 
E., closed-shell repulsions, H,, and polarization terms. These latter 
include the self-energies of polarization of the ions together with monopole— 
dipole and dipole-dipole interactions. In the present calculations we are 
interested in the difference between the polarization energy of the activated 
state and of the ground state. The most important part of this is the 
energy of polarization of the moving ion and the neighbours to the vacancy 
complex in the monopole fields acting onthem. Accordingly a polarization 
energy Ly = — 4) '0,F;7, in which q; is the polarizability ion i (Tessman et al. 
1953) and F, is the monopole field acting on it, has been included in the 
energy minimization programmes. A correction for dipole interactions is 
made after the relaxed configuration has been determined ; this correction 
is small. 

The results of these calculations are given in the table. Several points 
emerge. 

(i) AE for a jump into a vacancy pair is greater than or close to AE for a 
jump into a single vacancy. 

(ii) The contribution to AH from the closed-shell repulsions AZ, is 
relatively unimportant for the vacancy pair; as expected however it is 
larger for the large anions than for the cations. 

(iii) The Coulomb contribution AH, to AE is positive for the vacancy 
pair although it is negative for single vacancies; this is due to the absence 
of one of the two barrier ions, a and b, which exert a strong Coulomb 
attraction on the ion moving into the single vacancy. 

(iv) The moving ion is under the influence of a large electric field in the 
activated state and the associated polarization energy is an important 
part of the total. Since the anions are more polarizable than the cations 
the activation energies for anions are less: this contrasts with single vacancy 
movement. 


Mobility of Vacancy Pairs in Ionic Crystals 1161 


(v) Tracer diffusion via vacancy pairs will be limited by the less frequent 
of the two jumps. We must therefore compare the experimental value of 
1-1-1-2ev for Cl~ diffusion in NaCl with the higher of the two values 
obtained for NaCl, i.e. 1-46 ev, which actually governs the Na+ jump into 
the pair. Theory agrees with experiment that the pair is not highly mobile 
and the quantitative agreement is fair. We have evaluated the formation 
energy of vacancy pairs in our model and for this we obtained 1-27 ev which 
is only slightly less than the 1-3-1-4ev used above. The total activation 
energy for Cl- diffusion is then predicted to be 2:73 ev; the difference from 
the experimental value of 2-49 ev is perhaps within the combined error of 
the experiments and the theory. 

(vi) A repetition of our calculations for Cl- in KCl using the Born— 
Mayer potential employed by Dienes gives AZ=0-89ev. The compo- 
sition of this quantity shows the same features (i)—(iv) as the Born—Mayer-— 
Verwey potential. A value of AH =0-91ev was obtained using a recently 
proposed modification of the Born—Mayer constants which gave a potential 
of intermediate hardness (Tosi and Fumi 1961). The effect of increasing 
the hardness of the potential is to push the saddle point of the moving ion 
out to larger € into a region of unfavourable Madelung potential and this 
increases the Coulomb contribution to AE at the expense of AF, (table). 

(vii) It is difficult to analyse the disagreement of these results with those 
of Dienes since the two calculations were carried out differently. The 
energy of ground-state relaxation which we obtain agrees with that of 
Tosi and Fumi (1958) and disagrees with Dienes. A correction for this 
fact would however make Dienes’ predicted activation energy about zero. 
We also find a disagreement with his value for the polarization energy of 
the moving ion (his #,” for the rigid lattice). A correction here of about 
0-7 ev seems necessary and is then such as to bring the results more into. 
line with our own. 

(viii) In conclusion, it seems clear that the mobility of vacancy pairs in. 
alkali halides is less than or comparable with that of single vacancies. 
This implies a revision of the theories of colour centre aggregation at low 
temperature. 
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ABSTRACT 


The specimen independent infra-red absorption spectrum of diamond in 
the 2 1 to 6 region has been examined using a grating spectrometer. The 
results have been analysed in terms of multi-phonon processes. In particular 
the two-phonon absorption has been examined in detail in an attempt to 
deduce phonon energies for the various branches of the frequency spectrum 
at symmetry points in the reduced zone. Sets of lattice wave dispersion 
curves have been constructed from these results and are very unlike those 
for silicon and germanium. 


§ 1. InTRODUCTION 


THERE has recently been a considerable revival of interest in the lattice 
dynamics of diamond structure crystals. This was mainly stimulated 
by the work of Brockhouse and Lyengar (1958) who used inelastic neutron 
scattering techniques to determine directly the frequency (w) wave 
vector (q) dispersion relations for the plane wave normal modes of 
germanium. They obtained complete sets of dispersion curves only for 
the (100) and <{111) symmetry directions, but these were sufficient to 
demonstrate conclusively that the interatomic forces in this crystal 
contained an important long-range component. 

An analysis by Herman (1959) showed that one must include inter- 
actions between atoms at least as far apart as fifth neighbours to fit these 
results, and Cochran (1959) developed a theoretical model which demon- 
strated how such interactions could arise from coupling between dipole 
moments induced on the atoms when the lattice is distorted by a plane. 
wave. Subsequently Brockhouse (1959) determined the dispersion 
curves for lattice waves propagating along a (100) direction in silicon. 
These proved to be very similar to those for germanium, and Cochran 
(1960, private communication) has shown that his model can be used to 
fit the silicon data. 

It is evident that any completely satisfactory theory must reproduce 
the corresponding (w vs. q) dispersion curves for lattice waves in diamond 
itself. Unfortunately neutron scattering work is not at present possible 
owing to the large size of the crystal required, thus the actual dispersion 
curves are unknown for this crystal. One possible source of such 
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information was suggested to us by the work of Johnson (1959) on the 
infra-red lattice absorption bands of silicon. This was based on the 
theory of Lax and Burstein (1955) who showed that, although the direct 
absorption of radiation as single quanta of lattice vibrational energy 
(phonons) is forbidden, it may take place in processes involving the 
creation of two or more phonons. Using this theory he was able to fit 
the energies of all the peaks in his absorption curves using combinations 
of only four distinct phonon energies which he designated: transverse 
optical (T.O.), longitudinal optical (L.O.), longitudinal acoustic (L.A.), 
and transverse acoustic (T.A.). He argued that these energies were 
those at which the maximum density of modes occurred in each branch 
of the frequency spectrum which in turn would be expected to correspond 
to phonon wave vectors in the region of the surface of the reduced zone. 
Comparison with the neutron scattering data indicates that this argument 
is valid. Strictly there should be two T.O. and two T.A. energies since 
the two-fold degeneracy of these branches is only present for directions 
of high symmetry, but failure to resolve these implies that for most q 
vectors this splitting is not large. Thus it seemed that a similar analysis 
of the infra-red lattice absorption in diamond would provide correspond- 
ing information which, together with the known Raman energy 
(corresponding t0 wag=a,, at q=0) and elastic constants 
(McSkimin and Bond 1957), could be used to construct a set of dispersion 
curves. It was evident that earlier data (see Lax and Burstein 1955) were 
much too imprecise to be interpreted in this way with any degree of 
certainty and a new set of measurements was therefore undertaken. 
These were carried out at room temperature and no attempt was made 
to measure the variation of absorption with temperature since Collins 
and Fan (1954) have shown this to be small (a consequence of the high 
Debye Temperature) and therefore not likely to help in assigning the 
various peaks unless one is prepared to use very high temperatures where 
one will encounter difficulties arising from the onset of graphitization of 
the specimen. 


§ 2. EXPERIMENTAL 


Absorption measurements were made in the wavelength range 2 to 
6 on a 7500 lines per inch grating spectrometer with a Golay cell 
detector. No predispersing system was used; instead, second-order 
spectra were suppressed by a series of anti-reflected semi-conductor 
filters; germanium, lead sulphide and indium arsenide were used in the 
ranges 1-8—3-0, 3-0-3-8 and 3-8-6 respectively. The filters show good. 
transparency above their absorption edges (80% transmission) and very 
strong rejection below, thus stray light was negligible. Two diamonds. 
with good optical surfaces were measured—one Type I and one Type II a 
—giving identical results within experimental error. Reflection correc- 
tions were checked by measurements above the three-phonon cut-off 
(0-495 ev) where absorption in both specimens was negligible. Absorption. 
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coefficients were calculated, correcting for multiple reflections, from: 
(1— #)? exp (— Kd) 

where K denotes the absorption coefficient, 7’ the transmission, R the 
reflectivity, and d the specimen thickness. 


Tt — 


§ 3. ANALYSIS OF THE RzEsuLts (First STaGce) 
The absorption spectrum is shown in fig. 1 and can be divided into 
two regions: (a) the relatively strong bands lying below 0-330 ev, with 
a sharp cut-off at this energy, and (b) the weaker, more diffuse bands 
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Infra-red lattice absorption of diamond. (The broken line on (qa) indicates 
how the 0-251 ev peak and the 0-225 ev kink may be associated with the 
same absorption process.) 


above this energy with a rather less well-defined high-energy cut-off at 
0-495 ev. Since the known Raman energy is 0-165 ev this provides 
conclusive evidence that the maximum lattice frequency is that of optical 
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modes at q=0 and that regions (a) and (b) are due to two- and three- 
phonon processes respectively, the first commencing at exactly 2 x 0-165 
ev and the second at 3x 0-165 ev. The sharp cut-off of the two-phonon 
absorption is particularly striking. 

We shall confine our analysis to the two-phonon region since the 
broad three-phonon peaks are difficult to assign unambiguously. Further 
more, we shall assume that all bands are combination bands; i.e., 
hw =h(w,+w,) where fw is the energy of the observed peak and hw, 
and fiw, are the characteristic phonon energies involved. This is justified 
since (see Johnson 1959) difference bands (Aw =h(w,—wg) have essentially 
zero intensity when iw, and iw,> KT, as is the case for diamond at 
room temperature. 


Table 1 
Observed Assignment and Calculated 
energy of phonon energies energy 
feature in ev in ev in ev 
0-319-0-315 2 T.O. 0-316 
(shoulder) 2 x (0-158) 
0-302 T.O. +L.0. 0-302 
(peak) 0-158 +0-144 
0-281 T.O. +L.A. 0-281 
(shoulder) 0-158 +0-123 
0-267 L.O. +L.A. 0-267 
(peak) 0-144 +0-123 
0-251 T.O. +T.A. 0-251 
(peak) 0-158 +0-093 
0-244 L.O. +T.A. 0-237 
(peak) 0-144 +0-093 


Two-phonon cut-off at 
0-330 ev =2 x 0-165 ev: or twice the Raman energy. 


The energies of the various features are shown on fig. 1 and in table 1 
the energies of the main peaks are compared with those obtained by 
adding the following four energies in appropriate pairs: 

T.0.=0-158 ev, L.O.=0-144 ev, 

L.A. =0-123 ev, and T.A. =0-092 ev. 
It will be observed that our assignment is similar to that made by 
Johnson (1959) for silicon. There is no clearly defined 2T.O. peak, 
although the absorption curve exhibits a change of slope between 0-315 
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and 0-318 ev in the region of the 2T.0. energy (0-316 ev). This is not 
surprising since the T.O. energy is close to the Raman energy and the 
T.O. branches thus show little dispersion and consequently the selection 
rule (Lax and Burstein 1955), forbidding combinations of two degenerate 
branches, will ensure that the 2T.0. band is weak. Furthermore the 
T.0.+ L.A. band at 0-281 ev is so close to the 0-267 ev L.0.+L.A. band 
that it appears only as a shoulder on the high-frequency side of the latter. 
The most striking difference is the absence of the very strong T.O. + T.A. 
peak present in the silicon spectrum. The analogous peak for diamond 
at 0-251 ev is relatively weak and this fact provides a strong indication 
that the distribution of T.A. modes in diamond differs drastically from 
that in silicon, being considerably more diffuse. 

From table 1 it can be seen that we have used four energies to fit six 
features and have encountered the same difficulty as Johnson in that 
we are unable to deduce a T.A. energy which will fit both the 0-251 and 
0-244 ev peaks if we use the T.O. and L.O. energies determined from the 
other features. Also, since the distribution of T.A. modes is believed to 
be more diffuse, we expect a marked variation of the frequencies of the 
T.A. branches over the surface of the reduced zone and it is not clear 
from which regions the modes responsible for the T.0.+T.A. and 
L.0.+T.A. peaks come. 


§ 4. ANALYSIS (SECOND STAGE) 


Because of these problems, and also in the hope of establishing the 
other assignments more firmly, we have carried out the following more 
detailed analysis which is an extension of Johnson’s method and is aimed 
at associating specific features with specific symmetry points at the 
surface of the reduced zone. That this should be possible can be seen 
by considering eqn. (5.19) of Lax and Burstein (1955). This states that 
the absorption coefficient K(w) at a frequency w in the two phonon 
summation band at absolute zero is given by: 


2s? 
fe 24) — 9) 


K(w) o dla —w(q)—w(-q)] . . - (1) 
where H,,4 is a matrix element which couples modes of wave vectors 
q and —q from branches ¢ and ¢’ (assuming the wave vector of the 
incident photon is zero) and H,4=0 if t=t’. Van Hove (1953) has shown 
that the distribution function of normal modes p(w), which for a three- 
dimensional lattice is given by: 

ds; 


= w ——_ . . . es. (2 


(where the integrations are over surfaces S, in the reduced zone for which 
w,(q)=«), has a discontinuous first derivative whenever | y4(q) |=9. 
These points are known as critical points (¢.p.). 
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If we consider eqns. (1) and (2), the contribution to K(w), AK(a) 
from a ¢.p. at q in branch ¢ when w=a,(q)+,(—q) is given by: 


anv ne ie 
; aaa ee ae Sas ae 
kee aye nal (eo [v.oda) | 
or 
AK (Gy Geen aaa 
w(q)ay(— q) 


with a similar contribution from any singularity in the branch ¢’ at —q. 
Since there is no reason to believe that (0/dw) | H,% |? is discontinuous 
we see that each c¢.p. is associated with a kink in the absorption curve. 
These kinks do not necessarily coincide with absorption maxima and, 
in practice, they will tend to be smoothed to some extent. Phillips 
(1956, 1959) has shown how the minimum number of c.p. implied by 
symmetry may be located for a given lattice and in table 1 of the (1959) 
article has listed the minimal set for germanium. With the exception 
of the T.A. c.p. at the point & close to the zone surface in the [110] 
direction, these must all be present for diamond. The c.p. at & are 
only present because both T.A. energies there are assumed to be about 
twice that at X, the centre of the (100) zone face, which is equivalent to 
a point lying outside the reduced zone along the [110] direction. This 
implies that | vy ,,(q)| will vanish for both branches at some point 
along this direction. For diamond such c.p. may well be absent from 
one or both T.A. branches. We shall assume that the only c.p. present 
are those implied by symmetry, an assumption which the (1959) article 
indicates to be valid for all diamond structure crystals. Thus we see, 
from table 1 of this reference, that all branches have c.p. at the origin I, 
the centre of the (111) zone face L, the centre of the (100) zone face X, 
and the (140) zone corner W. The longitudinal branches are degenerate 
at the last two points, the transverse branches are two-fold degenerate 
at all c.p., and the optical and acoustic modes are separately degenerate 
at q=0, the former at the Raman frequency and the latter at w=0. 
For the present we do not know whether there are any T.A. (2) singulari- 
ties or not. 

In table 1 it now appears reasonable to assign the T.0.+L.0., 
T.0.+L.A., and L.O.+L.A. features to singularities associated with 
c.p. at L. Combinations involving two longitudinal modes from W or X 
are forbidden by the degeneracy at these points. Consequently we 
would expect the absorption curve to show dips in the regions of twice 
the L.(W) and L.(X) energies which probably correspond to those observed 
at 0-270 ev and 0-262 ev. It is significant that the latter lies at twice 
0-131 ev where there is a pronounced minimum in the defect activated 
single-phonon absorption (Smith and Hardy 1960) indicating a similar 
minimum in the density of longitudinal modes which should produce a 
0-262 ev minimum in the two-phonon absorption. It seems probable 
that this corresponds to 2hwL.(W) as |q| is larger for this point. The 
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0-270 ev dip thus lies close to 2haLl.(X), but, as the modes are degenerate 
for all q’s of the form (1, q, 0) |g |<4, one cannot say that this dip defines 
this energy precisely. However, it is unlikely to be more than 2 or 8 
greater than 0-270 ev. 

To analyse the T.0.+T.A. and L.0.+T.A. peaks we use Phillips’ 
(1959) analysis of the diamond specific heat data of Desnoyers and 
Morrison (1958) which indicates that the T.A. mode energy for the c.p. 
at Lis 0-063 ev. It is at once evident that the absorption peaks involving 
T.A. modes cannot be associated with the T.0.(L)+T.A.(L) and L.0.(L) 
+L.A.(L) energies which are 0-221 ev and 0-206 ev respectively, as the 
first of these lies close to the kink observed at 0-225 ev and the second is 
outside the range of our measurements. However, the presence of the 
0-225 ev feature does provide additional confirmation of the value of the 
T.O.(L) energy used previously and of the T.A.(L) energy derived by 
Phillips. 

Our main object now is to deduce the T.O. and T.A. energies at X 
and, if possible, at W, bearing in mind the fact that there may also be 
a T.A. cp. at &. There are three features which we can use for this 
purpose: the sharp peak at 0-244 ev, the peak at 0-251 ev, and the kink 
at 0-258 ev. ‘Two possible assignments are shown in table 2. In both 
cases the T.A. c.p. at X has been associated with the 0-258 ev kink and 
the 0-244 ev peak. The T.O. (X) energy has been assigned the maximum 
value consistent with that of 0-139 ev for the energy of the degenerate 
L.(X) modes. This last value is consistent with the condition 
2hwL.(X) > 0-270 ev derived earlier and is favoured by the distinctive 
features at this energy in imperfection activated single-phonon absorption 
(Smith and Hardy 1960). The T.O.(X) and T.O. (L) energies are now 
consistent with the somewhat greater dispersion of these modes along 
{100) as compared with 111) directions, which is characteristic of 
theoretical force constant models of the diamond structure (Herman 
1959). We cannot really decide in favour of either 2 (a) or 2 (b) on the 
basis of the present evidence. The former appears consistent, and 
the value of fiw(W) chosen for the T.O. modes is further supported by 
the fact that the T.0.(W)+L.(W) energy lies at 0-289 ev, close to the 
absorption minimum at about 0-290 ev, and one may thus see this 
feature as reflecting the minimum in the density of longitudinal modes 
believed to lie at 0-13lev. Unfortunately there is one serious draw- ~ 
back; for the T.A. modes fiw(W)<fw(X) which is implausible as 
| q(W) | > | q(X) |. It is this fact which leads us to favour the assign- 
ment 2b. This has been made using the same values for the T.O. and 
T.A. mode energies at = and W, an assumption supported by Phillips’ 
(1959) work, and the same energy (0-158 ev) for the T.O. (2, W) modes 
as that used in assignment 2a. The T.A.(x, W) mode energy is then 
chosen so that the L.(=, W)+T.A.(Z, W) combination lies at 0-247 ev 
and can be identified with the commencement of the 0-244 ev peak. The 
0-225 ev kink is now unambiguously associated with the T.O.(L) + T.A.(L) 


P.M. 4F 
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energy, and the T.A.(W)+T.0.(W) process now lies at about 0-274 ev 
in which region there is some slight evidence of a kink. (The 0-247 ev 
minimum, which appears as part of this assignment, is also to be expected 
on the grounds that the 2 L.A. combination is forbidden everywhere in 
the reduced zone and, in particular, for the 2 L.A.(L) energy or 0-246 ev.) 


Table 2 


Calculated 
energy in ev 


Assignment and 
phonon energies in ev 


(a) Keature 
(energy in ev) 


0-258 
(Kink) 


0-251 
(Peak) 


0-244 
(Peak) 


0-225 
(Kink) 


TO Pro 
0-153 +0-105 


T.0.(W) +T.A.(W) 
0-158 +0-093 


Lik) era 
0-139 +0-105 


L.(W) +T.A.(W) 
0-131 +0-093 
or 
T.O(L) een by 
0-158 +0-063 


No T.A.(X) e.p. 


(6) 0-258 
(Kink) 


0-251 
(Peak) 
0-244 
(Peak) 
0-247 


(Minimum) 


(0-274 


Kink ?) 


aNeH Doe Pear Se 
0:153 +0-105 


(not associated 
with c.p.) 


Li ( here TA 
0-139 40-105 


LW)  +'T.A.(W) 
0-131 +0-116 


T.0.(W) +T.A.(W) 
0-158 +0-116 


0-258 


Both these last assignments are uncertain. The corre- 
sponding T.A.(x) singularities probably occur at almost 


the same energies. 


It is evident that neither assignment can be sustained in detail by the 
present data. If one relaxes the condition fwL.(W)< hwL.(X), 
assignment 2a with W and X interchanged is still possible. However 
if this were correct, the resultant dispersion curves would agree with those 
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in fig. 2. to within the limits of error. However, one can say with 
reasonable certainty that the T.0.(X)+T.A.(X) energy lies close to 
0-258 ev, the sharp rise in absorption at this point reflecting a similar 
rise in the density of T.A. normal modes to be expected at the T.A.(X) 
energy. Also the use of this energy to fit the 0-244 ev peak to 
T.A.(X)+1L.(X) leads to an L.(X) energy (0-139 ev) which is consistent 
with the previous results. 


§ 5. CONCLUSIONS 


The results of this analysis have been used to construct the set of 
phonon dispersion curves which are shown in fig. 2. The slopes of the 
acoustic branches at w=0 have been derived from the elastic constants. 
measured by McSkimin and Bond (1957). The curves have been inter- 
polated using the zone boundary energies derived from our results, and 
the probable errors in these are indicated by the short vertical lines. 
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Lattice wave dispersion curves for diamond constructed from the zone 
boundary phonon energies derived from the absorption data. The short 
vertical lines indicate the probable errors in these energies. 


Also shown are the T.A. branches obtained by scaling those of 
germanium according to the homology rule (Brockhouse 1959): 
wa(Ma?)-12, where M is the atomic mass and a the lattice parameter. 
(In the region of w=0 the slopes have been fitted to those predicted by 
the elastic constants of diamond.) The discrepancy between these and 
our results is obviously very large and demonstrates conclusively that 
the interatomic forces in the two crystals are very different. The 
dispersion curves for diamond imply that long-range forces are much less 
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important, since it is these which are responsible for the anomalously 
high dispersion of the T.A. branches in germanium. 

Our results are obviously somewhat tentative, and require con- 
firmation by the more direct neutron scattering technique. Such 
confirmation could well be provided by very low resolution scattering 
work such as may be possible with crystals of a readily available size. 
This would define the phonon energies at ¢.p. to within 10 or 20°% which 
is all one needs to check the correctness of our assignments. If these are 
valid, then the energies we have derived from infra-red absorption are 
considerably more accurate than the neutron scattering results. 
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ABSTRACT 


The pressure and temperature gradients accompanying heat flow in liquid 
helium II have been measured simultaneously. The experiments were 
made in a glass tube of 107-6 ~ diameter and 10:2 cm length between 1-3°x 
and 1-8°K. Significant deviations from Allen and Reekie’s rule were 
observed and are provisionally attributed to an eddy viscosity, ye, due to tur- 
bulence in the superfluid component. An equation, based on dimensional 
arguments, is proposed for 7 which satisfactorily represents the experimental 
data available at present. 


§ 1. [yTRODUCTION 


ACCORDING to an empirical rule of Allen and Reekie (1939) the pressure 
gradient observed in a steady-state heat conduction experiment in helium II 
is proportional to the heat flow. The proportionality holds good, at least 
approximately, even at relatively high values of the heat current density 
despite the fact that the temperature gradient is then found to be a strongly 
non-linear function of the heat flow. 

An explanation of the Allen and Reekie rule was given by Gorter and 
Mellink (1949) who postulated that the rapid increase in the observed 
thermal resistance at high heat current densities is caused by a non-linear 
mutual friction force acting equally and oppositely on the superfluid and 
normal components. With such a mutual friction term included in the 
hydrodynamic equations of both superfluid and normal component it is. 
found that its effect on the pressure gradient cancels, and the pressure 
difference Ap observed in a cylindrical flow channel, for all values of the 
heat flow W, should be given by 

Va 8ynlW 
— raApST 


where 7p is the normal viscosity, pS7T the heat of transport per unit volume, 
and a and / the radius and length of the tube. The Gorter—Mellink cal- 
culation assumes that the mutual friction depends only on the temperature 


(1) 
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and relative velocity v between the two fluids, that there is no restriction 
on the value of the superfluid velocity vs at the wall of the channel, and that 
vs may be rotational without producing any additional dissipative or 
viscous effects. 

Recent work on non-linear processes in helium II has indicated that the 
Gorter—Mellink mutual friction is almost certainly due to scattering of the 
normal fluid excitations by turbulence in the superfluid, the turbulence 
consisting of a random distribution of the quantised vortex lines proposed 
by Feynman (1955). On this picture it seems likely that other dissipative 
processes may occur in the superfluid besides the mutual friction, for 
example it may be possible for the superfluid to interact with the channel 
walls by the attachment of vortex lines, or, as was suggested by Vinen 
(1955, and private communication), there may occur some radial transfer 
of momentum in the superfluid by an eddy viscosity analogous to that found 
in ordinary turbulence. The effect of such processes may be detected and 
separated from the mutual friction by testing the validity of eqn. (1) when 
the helium is in a turbulent condition. 

Small deviations from the equation at large heat currents have already 
been reported in a re-analysis of the earlier data of Mellink (1947) by Vinen 
(1955) and Lifshitz and Andronikashvili (1959), and in some experiments 
on narrow channels by Keller and Hammel (1960). Some of the present 
measurements, which were made on the heat flow in a glass tube of 107-6 
microns diameter, have already been briefly described by us (Brewer and 
Edwards 1958). Another series of experiments of the same kind, but using 
a range of tube sizes, has since been carried out in this laboratory by 
Bhagat and Critchlow (1961). 

The present paper is one of a series describing heat conduction by 
liquid helium II in capillary tubes. Previous papers, following a prelim- 
inary report (Brewer, Edwards and Mendelssohn 1956), have dealt with the 
subcritical region (Brewer and Edwards 1959) and with the transition to 
non-linear conduction (Brewer and Edwards 1961a). In a forthcoming 
paper we shall describe and analyse measurements of the thermal resistance 
in the non-linear region. 


§ 2. APPARATUS 


A diagram of the apparatus is shown in fig. 1. The glass capillary tube 
C communicates at its upper end through a wide metal tube with the 
constant-temperature liquid helium bath. The other end is joined with 
platinum seals through a metal bridge tube to another glass tube S of 
diameter 0-8 mm in which the liquid height can be measured with a catheto- 
meter. The arrangement of the tubes C and S shown in the diagram was 
chosen to allow the maximum possible difference AH ops between levels in S 
and the bath (which may be positive or negative), with the restricted length 
of cryostat available. The tubes are thermally isolated from the bath bya 
metal vacuum jacket V and a glass jacket W, joined together with a 
copper-glass seal. The vacuum space is continuously pumped during 
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experiments through a metal tube which has been omitted from the diagram 
for simplicity. A temperature gradient across C can be established by the 
heater h, and measured by the carbon resistance thermometer ¢t, which is 
glued in thermal contact with the metal tube connecting C and 8. The 
capillary tube C has a diameter 2a = 107-6 microns and length J/=10-2cm; 
its Poiseuille constant, 7a4/8/, was also determined by gas flow measurements 
as previously described (Brewer and Edwards 1959). On applying small 


Fig. 1 


The apparatus. 


heat inputs, the observed level difference, AH ons, increases until an equili- 
brium is reached, when the mass flow of superfluid down the capillary is 
balanced by the upward flow ofnormal fluid. The total pressure difference 
AF is then given by the level difference AH ops, together with the difference 
in vapour pressure between the liquid in S and the bath. The latter was 
calculated from the measured temperature difference between S and the 
bath, using the T;;, vapour pressure scale. At high heat currents the 
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vapour pressure correction becomes so large that the observed level 
difference AH ops is negative. 

To measure the bath level accurately, a second tube of the same diameter 
as 8 was used, dipping directly into the bath. Appropriate corrections for 
surface tension rise were made from separate experiments in water in which 
the rise was measured at different parts of the tube. 


§ 3. RESULTS AND DISCUSSION 
Figure 2 shows the experimental results for the temperature 1-563°K, 
and illustrates the importance of the correction to the observed pressure 
gradient for the increase in vapour pressure. The directly observed 
quantity AHops/W decreases sharply at higher heat currents, whereas the 
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Experimental results at 1-563°K, showing the directly observed quantity 
NE opal W, the corrected values AH/W, and the thermal resistance AT'/W, 
as functions of heat flow W. The arrow indicates the critical heat current 


HA ° tale determined for the same tube (Brewer and Edwards 
a). 


corrected quantity AH/W increases. Values of the temperature difference 
across the tube have been divided by heat current to give the thermal resis - 
tance AT'/W which is also plotted in fig. 2. The thermal resistance is 
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independent of W at low heat currents, and then shows a transition to 
strongly current-dependent values where mutual friction forces occur. 
These observations of A7'/W are in complete agreement with those pre- 
viously obtained with the same tube in a different expermmental arrange- 
ment (Brewer and Edwards 1959, 1961a,b). The arrow on the curve gives 
the critical heat current W, as determined in the previous experiments, and 
is seen to be consistent with the present measurements. 

During these experiments some fluctuations in the temperature and 
pressure were observed at intermediate heat currents just above critical, 
which are probably related to those found in our earlier thermal resistance 
measurements. None of the hysteresis observed earlier was found in the 
present experiments, possibly because there was appreciable vibration 
present, or because of small obstructions or roughness in the tube, or: 
because in this experiment the tube was connected directly to the bath. 

If Allen and Reekie’s rule were followed, the graph of AH/W would be a 
line roughly parallel to the W-axis, but bending down slightly when the 
temperature gradient becomes large. This is seen to be true up to heat 
currents in excess of the critical current W. for the appearance of mutual 
friction. Large deviations then occur quite sharply, and at still higher: 
currents AH/W begins to level off once more. 

We define an effective viscosity in these experiments given by 


map ST Ap @) 
Tei ee 


YHett = 


which may be considered to include dissipative effects occurring in both. 
normal and superfluid components, but which, by the theory of Gorter and 
Mellink, excludes the mutual friction simultaneously present. This. 
effective viscosity is plotted against the relative velocity of the two fluids, v, 
for three temperatures in fig. 3. At low heat currents, yer: is accurately 
equal to the viscosity of the normal fluid, yn, as previously determined in 
the same tube (Brewer and Edwards 1959), so that in this region eqn. (1) is. 
obeyed. The slight decrease of yer; with increasing velocity is due to the 
variation of 7» and pST with temperature along the tube, and agrees with 
the theoretical expectation, shown by the broken curves in the figure. 
These broken curves have been calculated from the known temperature. 
variation of pST and ny and the observed temperature gradients in each 
experiment, by a method described in our previous paper quoted above. 
The quantity represented by the broken curves, which we term 7°", 
therefore represents the expected result of an experiment in which eqn. (1) 
is exactly true and in which the temperature gradient along the tube is the. 
same as in the actual experiment. Arrows in fig. 3 show the critical velo- 
cities for appearance of mutual friction already determined in our previous. 
experiments (Brewer and Edwards 1961 a), and are found to be consistent 
with the present observations of A7/W,as was illustrated in fig. 2 for 1-563°K. 
Thus fig. 3 shows that large contributions to 7ert from causes other than the- 
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viscosity of the normal component or mutual friction occur at velocities 
near but generally rather larger than the critical. - 

To interpret the data we assume provisionally that the additional 
contributions to yer are due to an eddy viscosity of the superfluid, 7e, since 
this seems reasonable both from theory and from other experiments (see, 
for example, Bhagat 1960, Bhagat and Mendelssohn 1961). If we make 


Fig. 3 
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Effective viscosity neg plotted against relative velocity v of superfluid and 
normal fluid. O, 7’=1-306°K; A, 7 =1-563°K; [], 7 =1-806°x. The 
broken curves represent values of 7,°° (see text). Arrows indicate the 
critical velocities for appearance of mutual friction previously determined 
in the same tube (Brewer and Edwards 1961 a). 


the further, simplifying assumption that there is no restriction on the 
value of vs at the wall of the channel, then it can be shown from the two 


fluid equations of motion that ye is given by the difference between the 
broken and solid curves in fig. 3, that is 


Ne=Nett—n"- ofthe Lie meta) 
As mentioned in the Introduction, the eddy viscosity is understood to 
arise from the motions of the random distribution of vortex lines in the 
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superfluid component, and it may be expected to depend rather directly on 
the length of vortex line per unit volume present in the liquid. A plausible 
equation for 7 from dimensional arguments is 


ne=p =, (DL!) SMe eerie oh enti ace 5 (4) 


where 27(f/m) is the quantum of circulation around a vortex line, p the 
liquid density and f(DZ,?) a dimensionless but possibly temperature 
dependent function of vortex line density, LZ), and tube dimension, D. 
(The length D will be defined as the ratio of tube cross section to perimeter 


Fig. 4 
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Eddy viscosity 7. plotted against the dimensionless number DL,'” (see eqns. 
(4) and (5) and text). O, 7=1-306°K; A, T =1-563°K; 0, 7 =1-806°K. 
Closed circles @ represent results calculated from the measurements of 
Bhagat and Critchlow (1961). 


so that it may be applied to tubes of various shapes.) There is no need to 
includeany dependence on the fluid velocity v in eqn. (4)since there is another 
relation between L, and v determining the magnitude of the mutual friction. 
Values of L, have been derived from the present measurements using the 
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theory of mutual friction developed by Hall and Vinen (1956) and Vinen 
(1957) and the experimental values of the temperature gradient. In fig. 4 
we plot ne against the dimensionless number (DL,"*), in order to determine 
the form of the functionf. Itisinteresting that the value of ye is apparently 
zero up to (DL,12)~1 and then increases almost linearly with DE 
Moreover, the magnitude of ye at a given value of (DL,"*) does not depend. 
strongly on temperature. Since Ly? is of the order of the spacing 
between vortex lines, values of (DZ,!/?) equal to or less than unity represent 
situations where the average line spacing is of the order of or greater than 
the dimensions of the channel. It is natural to expect that the behaviour 
of both the mutual friction and eddy viscosity should change in this region. 
Measurements on the mutual friction at low heat current densities (Brewer 
and Edwards 1961b) have shown that the vortex line theory of Vinen, 
based on the assumption of a homogeneous, random distribution of vortex 
lines, breaks down (as may be expected) when (DL,') = 1, and the variation 
of L, with v is then completely different from that at higher vortex line 
densities. 

The present experiments can therefore be explained by postulating an 
eddy viscosity in the superfluid, with no constraint on the value of vs at 
the channel wall. The magnitude of this proposed eddy viscosity is given 
approximately by 


h nee 
ne® & p — [(Dio)— Cal, DLS Cen) Se. ee 


where ¢, and ¢, are dimensionless numbers (¢, + 0-05, €, + 1), only slightly 
dependent on temperature. As can be seen from fig. 4, there are systematic 
deviations from this linear equation for small values of (DZ,"?) at the two. 
higher temperatures. 

Equations (4) or (5) predict a definite relation between the proposed 
eddy viscosity y, and tube size. We are unable to verify this prediction 
from our own measurements since only one tube of 107-6 micron diameter 
was investigated. However, Bhagat and Critchlow (1961) give simultane- 
ous values of the mutual friction and pressure gradient for a tube of 297 
microns diameter at 1-41°K, which have been plotted as closed circles in 
fig. 4. Itis satisfactory that these independent points also lie on a straight 
line and are in good agreement with the values at our closest temperature, 
1-306°K. The only other available data which might be used to test eqn. 
(5) or the more exact eqn. (4) are the results of Mellink (1947) re-analysed by 
Lifshitz and Andronikashvili (1959); these are in rather rough agreement 
with our proposed equations, but since they refer to radial flow in a disk- 
shaped slit, they are not exactly comparable with our own. The measure- 
ments of Keller and Hammel (1960) were made with very large temperature 
differences across the tube, and are again rather difficult to in terpret. 


ACKNOWLEDGMENTS 
The authors are grateful to Dr. K. Mendelssohn, F.R..S., for his interest. 
in the work, to Dr. W. F. Vinen for some stimulating discussions and. 


Heat Conduction by Liquid Helium IL in Capillary Tubes 1181 


particularly for pointing out the importance of the Allen and Reekie 
rule, and to Dr. 8. M. Bhagat for a copy of his paper with Dr. P. R. Critchlow 
before publication. The work has been supported by Research Fellowships 


from the Nuffield Foundation (D.F.B.) and the Pressed Steel Company 
(D.0.E.). 


REFERENCES 


ALLEN, J. F., and Renxts, J., 1939, Proc. Camb. phil. Soc., 35, 114. 

Buaeat, 8. M., 1960, Proc. phys. Soc., Lond., 75, 303. 

Buaaat, 8. M., and Crircntow, P. R., 1961, Cryogenics, 2, 39. 

Buragat, 8. M., and Menpexssonn, K., 1961, Cryogenics, 2, 34. 

Brewer, D. F., and Epwarps, D. O., 1958, Low Temperature Physics and 
Chemistry (Madison: University of Wisconsin Press), p. 12. 

Brewer, D. F., and Epwarps, D. O., 1959, Proc. roy. Soc. A, 251, 247; 1961 a, 
Phil. Mag., 6, 775; 1961 b, Phil. Mag. (in print). 

Brewer, D. F., Epwarps, D. O., and MenpELssoun, K., 1956, Phil. Mag., 1, 
1130. 

FrynMan, R. P., 1955, Progress in Low Temperature Physics (Amsterdam: 
North-Holland Publishing Co.), 1, 36. 

GorTER, C. J., and Metrix, J. H., 1949, Physica, 15, 285. 

Hat, H. E., and Vinen, W. F., 1956, Proc. roy. Soc. A, 238, 215. 

Keer, W. E., and Hammet, E. F., 1960, Ann. Phys. (U.S.A), 10, 202. 

Lirsuirz, E. M., and ANDRONIKASHVILTI, E. L., 1959, 4 Supplement to Helium 
(Chapman and Hall), p. 121. 

MeEtuink, J. H., 1947, Physica, 13, 180. 

Vinen, W. F., 1955, Cong. de Phys. des Basses Temp., Paris, p. 60; 1957, Proc. 
roy. Soc. A, 242,493; Ibid. A, 243, 400. 


Femi soin| 


CORRESPONDENCE 


The Determination of Dislocation Densities in Thin Films 
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I wish to describe (1) a new method of estimating the density of 
dislocations observed in a thin metal film by electron microscopy, and 
(2) an investigation of the sources of random error in estimates of the 
mean density corresponding to a particular treatment. 

(1) Bailey and Hirsch (1960) have described a method of estimating 
the density of dislocations. They measure the total projected length R., 
of dislocation line in a given area A on a typical micrograph. Then, on 
the assumption that the dislocation segments are randomly orientated 
with respect to the plane of the film, the dislocation density is 
p=(4/7)R,/At, where t is the thickness of the film. 

A disadvantage of this method is the labour of measuring R&,. 
However, an estimate of R, may be obtained easily by making use of 
the results of Smith and Guttman (1953). They showed that if a set 
of random lines with a total length L is marked on A, and the number 
of intersections V which dislocations make with the grid lines is measured, 
then R,=7NA/2L, giving p=2N/Lt, provided N is large enough. 

The two methods of finding R,, have been compared for 20 micrographs 
taken from thin films of a specimen of aluminium cold-rolled to about 
9% reduction. The new method proved to be about ten to twenty times 
faster in measuring time alone. Five lines drawn in random directions 
on a plate taken at x20000 gave enough intersections (about 50), 
and it was not necessary to enlarge the micrographs, as it was for the 
direct measurement of R,. The mean difference between the dislocation 
densities (of about 2x 10°/cm?) estimated by the two methods was 
~4-1x107/em?, the standard deviation of the mean difference being 
~ 3-6 x 107/em?. Hence there was no significant difference between the 
two methods. 

(2) As part of another programme, the new method has been used in a 
preliminary investigation of the sources of random error in determinations 
of dislocation density. Three specimens of aluminium were prepared at 
each of three levels of deformation, namely 18, 45 and 88% reduction, by 
cold rolling. A thin film was prepared from each specimen and two 
micrographs taken from each of two grains. The specimens were treated 
in three groups, all levels being represented in each group. In this way 
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systematic variation during the measurement was shown to be absent. 
The dislocation density was determined for each micrograph, and an 
analysis of variance made of these data. In what follows, references 
have been given to specific sections of Davies (1958) which describe 
the procedure. For each level of deformation the specimens, grains, and 
micrographs form a hierarchic classification with three sources of 
variation (§ 6.3). The additive property of variance (§ 6.11) was used 
in the separation and estimation of components of variance associated 
with specimens, grains, and micrographs. These estimates are given in 
the table. The variance of the mean dislocation density (§ 6.11) based 


Components of variance 
Level Between Between Between 
specimens, Ss? grains, Sg” micrographs, Sm? 
18 0 0-87 0-68 
45 0 6-03 2-14 
88 1-07 0-17 1-23 


on Ng Specimens, Ng grains per specimen and nm micrographs per grain is 
S2/ngt+S,?/Nsngt+Sm?/nsrgnm. Hence it is seen from this formula and 
the estimates in the table that at levels 18 and 45 the most efficient way 
to reduce the error in the mean is to take more than two grains per 
specimen ; two specimens and two micrographs per grain are enough. 
At level 88, however, the most efficient procedure is to take more 
specimens. The variance estimates are fully efficient, but since they are 
based on few observations these suggestions should be regarded as a guide 
to the design of further investigations. 
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§ 1. INTRODUCTION 


Ir has been reported on several occasions (Ogilvie 1959, Gillam 1959 and 
Ogilvie and Thomson 1961) that the bombardment of a metal surface with 
positive ions of low energy (from 4000 to 10 ev) results in the formation of 
crystallites which are considerably misoriented with respect to the crystal 
structure before bombardment. This was inferred from experiments in 
which flat specimens were electropolished to form holes, and transmission 
electron diffraction patterns were taken through thin regions of the metal 
round the edge of a hole. The pattern from a polished but unbombarded 
specimen consisted of sharp spots, if it was large grained or a single crystal. 
In the diffraction pattern of the same specimen after ion bombardment, the 
spots had become arced, and the pattern was similar to those obtained by 
Hirsch e¢ al. (1955) in beaten and etched gold foil. A misorientation of 35° 
has recently been reported when 130ev argon ions bombarded a silver 
single erystal (Ogilvie and Thomson 1961). 

More detailed investigations of these effects can be carried out by 
transmission electron microscopy. The technique of observing disloca- 
tions and other defects used by Whelan et al. (1957), Bailey (1960) and many 
others, seems a promising way of finding out how ion bombardment causes 
these changes in the crystalline texture. Dislocation loops have been 
shown to arise from the neutron bombardment of copper (Silcox and Hirsch 
1959) and from the bombardment of copper with high energy «-particles 
(Barnes and Mazey 1960) in this way. In all the experiments mentioned _ 
so far, ‘electropolishing to holes’ was used to produce the thin films 
required in the microscope. 

In the experiments to be described, evaporated single crystal gold films 
were used instead of electropolished specimens. These had the advantage 
of being of uniform thickness, and a large area of such specimens was thus 
available for inspection ; with electropolished specimens the part examined 
is wedge-shaped, and only small sections near the holes are available. 
Extinction contours and impurities from the polishing bath can also be a 
nuisance with polished specimens. Evaporated films are frequently 
objected to because they have different properties from bulk specimens, but 
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the misorientation effects due to ion bombardment should also appear in 
films, as the effects are thought to occur in a layer of about 100 A thickness 
and this should be observable in a film some 500 4 thick after bombardment. 


§ 2, EXPERIMENTAL METHOD 


Single crystal gold films were prepared by the method of Bassett and 
Pashley (1959). A layer of silver at least 1000 A thick was evaporated on 
the (100) face of a freshly cleaved rock-salt crystal held at 250-300°c, and 
the gold was immediately evaporated on the silver. The rock-salt was 
dissolved in water and the silver in concentrated nitric acid, and a coherent 
single crystal gold film remained. Gold films made in this way were 
examined in a Siemens Elmiskop 1 electron microscope and found to be of 
the same general appearance as that described by Bassett and Pashley. 

The ion bombardment of the film took place in a low pressure argon 
discharge similar to the type used in the experiments mentioned previously. 
A cylindrical anode of aluminium with open ends was used with a magnetic 
field of about 500 gauss directed along its axis. A hot tungsten filament 
at one end of the cylinder and the specimen at the other acted as cathodes. 
The greased cone and socket joins in the glassware used in earlier experi- 
ments to make the system demountable were abandoned, as evidence of 
carbon contamination appeared in the first micrographs. Thus the system 
could be baked and the gaseous impurities in the glass and electrodes 
pumped away before the argon was introduced. The baking temperature 
was 500°C. The pressure during the ion bombardment was 0-6 x 10-3 torr 
which was sufficiently low to prevent sputtered material being scattered 
back to the specimen. 

Backus (1949) discusses this discharge system, in which most of the ions 
bombard the specimen normal to its surface with energies slightly less 
than that given by the anode voltage. The ion bombardment with a 
discharge voltage of 300 volts lasted for 1 hour and the current density 
was 0:15ma/sq.cm. The film which was initially 2000 4 thick appeared to 
be much thinner in the microscope and was probably about 500A thick. 
This would correspond to an approximate sputtering ratio of the order of 
0-1 atoms removed per incident positive ion, under these conditions. 


§ 3. Rusuts 


Figure 1} is a micrograph showing the appearance of a film after the 
above treatment. The whole of the field is covered with specks and there 
are a few dislocation rings. Barnes and Mazey (1960) found similar spots 
after a-particle bombardment of copper. Figure 2 shows the diffraction 
pattern taken from the same part of the specimen, and the noticeable 
feature is that the arcing of the spots which had been expected did not 
appear. There were no edges in the sections of the film examined in this 
case. 


{ All figures are shown as plates. 


Correspondence 1187 


Some of the films, however, lay across the grids in such a way that an 
edge of the film was seeninagridhole. Figure 3 shows a diffraction pattern. 
taken from an area including such an edge, and misorientations of the kind 
found previously in diffraction patterns are seen. It seems therefore that. 
these are a feature, not of the ion bombardment, but of the twisting and 
bending of the specimen in the very thin fragile region of anedge. Moving 
to an adjacent grid hole on the same specimen which did not contain an 
edge, diffraction patterns were obtained similar to fig. 2. For a gold film 
bombarded in a 130 volt discharge, the observations were similar. 

Further support for the view that edge effects have been misleading was 
obtained by making an evaporated gold film in which a wire placed in front 
of the rock-salt prior to the evaporation resulted in a thin edge being formed 
on the unbombarded film. Once again, well away from the edge, the diffrac- 
tion pattern was that of a single crystal, but fig. 4 shows the ‘misorientation’ 
effect occurring near the edge. 

It remains to be explained why the misorientations have not been found 
before ion bombardment in diffraction patterns of films made by electro- 
polishing, but have been observed after ion bombardment. This is 
possibly because the thin areas round the electropolished hole were in fact 
not thin enough, and the edges were not at this stage bent. Jon bombard- 
ment caused a further thinning of the edge and may then have also been 
instrumental in bending it. 

The nature of the defects from which the specks in fig. 1 arise has not 
been studied in detail. Ion bombardment is expected to produce vacancies 
and interstitials in equal numbers, some of which can combine to form 
clusters. Collapse of the lattice and the formation of dislocation loops 
can then occur. Annealing at 350°c for 10 min resulted in a disappearance 
of the specks and the formation of complex dislocation networks. Figure 
5 shows a resulting micrograph. 
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THERE is considerable interest in the problem of the electrical resistance 
due to stacking faults in pure metals (for references see, for example, 
Howie 1960) but at present the experimental evidence for the magnitude 
of this resistance is scanty. We have already noted that in pure sodium 
the dominant residual electron scattering is not due to stacking faults 
(Basinski ef al. 1959) ; a consideration of the electrical resistivity of the 
low-temperature phase of lithium leads us to conclude that stacking-fault 
scattering is small in lithium also. 

When lithium is cooled below about 80°K a hexagonal close-packed 
phase begins to form by the martensitic mechanism (Barrett 1947, 1951, 
1956, Barrett and Trautz 1948). According to Barrett ‘it would be 
appropriate to describe the structure of the spontaneously formed 
low-temperature phase as a nearly random placement of each successive 
layer of atoms in the hollows of the layer just below, or alternatively as 
an f.c.c. structure with very imperfect stacking of (111) planes” (1951, 
p. 354). This implies that the stacking-fault density may be nearly as 
high as one plane in three. At 4°K approximately 90° of the material 
has transformed to this ‘ hexagonal ’ phase, the rest being untransformed 
b.c.c. material (Hull and Rosenberg 1960, Barrett 1951, Martin 1960). 

Typical values for the resistivity of such ‘ pure’ lithium specimens 
at 4°x, where the scattering is essentially all residual scattering, are 
between 1 and 1:5x10-8ohmecm. If we assume that the residual 
resistivities of the two phases are equal and attribute all the resistivity 
of the close-packed phase to stacking faults then, since the mean free 
path of electrons at room temperature is about 10-* em, we deduce that 
at 4°K it is about 10-*cm in such samples. In the extreme case that _ 
every third plane is a faulted plane the distance between stacking faults 
is about 10-7cm and the probability of reflection at an individual 
stacking fault is approximately 10-4t. In this estimate we have taken 
the free electron Fermi velocity for the electrons; this should not 
introduce a serious error. We are also assuming the scattering from 
one stacking fault to be independent of that from another ; according 
to Howie (1960) this should be true if they are further apart than one 
lattice spacing. 


+ Now at the H. H. Wills Physics Laboratory, University of Bristol. 
+ We are assuming that the stacking faults are effectively infinitely wide. 
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It might be argued that perhaps in the resistive specimens less than 
90%, of the b.c.c. phase transforms and, being free from stacking faults, 
short-circuits the h.c.p. phase. If this were a large effect it would show 
itself in the temperature region around 80°K where we can compare the 
resistivity of the two-phase mixture with that of the pure b.c.c. phase. 
We deduce from such measurements (cf. Dugdale and Gugan 1961) 
that even if only 50°% transforms the reflection coefficient can hardly 
exceed 10-3. In fact, however, there seems no reason to doubt the figure 
of about 90%. The. amount transformed appears to be insensitive to 
the amount of strain in the specimen before transformation (Martin 1960) 
and, moreover, our resistance measurements were made on a variety of 
specimens (bare wires with diameters between 0-5 and 5mm and with 
different heat treatments) which all gave similar results. 

As noted above, we concluded earlier that stacking-fault scattering 
in sodium was small but this was hardly surprising because even in the 
hexagonal phase this metal approximates well to a free-electron metal 
(Dugdale and Gugan 1960). On the other hand, it has been argued that 
in a metal with a distorted Fermi surface the reflection coefficient should 
be fairly high (cf. Ziman 1960), and for copper a value of about 0-2 
has been calculated (cf. Howie 1960). The Fermi surface in lithium is 
thought to be highly distorted in the b.c.e. phase (cf. Cohen and Heine 
1958) and in the close-packed phase the situation may well be similar. 
It is therefore interesting that the reflection coefficient for the scattering 
of electrons by stacking faults is apparently so small in lithium. However, 
as the next letter shows, it is perhaps possible to account for the small 
value on the basis of existing theories of the scattering of electrons by 
stacking faults. 
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A Possible Explanation of the Results of the Preceding Note 
Regarding Scattering of Electrons by Stacking Faults in Lithium 


By A. Howrs 


Cavendish Laboratory, Cambridge 
[Received July 5, 1961] 


In the preceding letter it is reported that samples of h.c.p. lithium, 
containing an estimated density of 107 cm—! of independently scattering 
stacking faults, have a resistivity of about 1x10-8ohm-cm. The 
resistivity to be associated with unit stacking-fault density is therefore 
1 x 10-* ohm-cm and may be compared with a value of 2 x 10-8 ohm-em 
recently measured for stacking faults in gold (Cotterill 1962). This latter 
figure is in quite good agreement with the lower of two estimates made for 
copper by Howie (1960) on the basis of a wave-matching calculation and 
is thought to be typical of copper, gold and silver. The apparent difference 
of about two orders of magnitude between the resistivity of stacking faults 
in lithium and of stacking faults in the noble metals may result from several 
causes. 

(1) Lithium, as the authors of the preceding note point out, is regarded 
as being less free-electron-like than sodium. It is however almost certainly 
more free-electron-like than the noble metals. Estimates of the size of 
the band gap at the (110) Brillouin zone boundaries in b.c.c. lithium 
(Ham 1960) seem to be only about half to one-third of the value for the 
gap at the (111) zone boundary in copper (Pippard 1957). This would 
imply that in lithium the value (suitably averaged over the Fermi surface) 
of the most important diffracted wave intensity is between one-quarter 
and one-ninth of what it is in copper and the resistivity of stacking faults 
which is essentially a diffraction effect might be expected to fall by a 
similar factor (Ziman 1958). 

(2) In the noble metals, the electrons most strongly scattered by a 
stacking fault on a (111) plane were found (Howie 1960) to be those 
diffracting from (111), (111) and (111) planes, i.e. electrons travelling at 
about 70° to the direction of normal incidence on the fault. In the h.c.p. 
structure however the situation is quite different. Electrons travelling 
normal to the fault and diffracting from planes parallel to the fault would 
be unaffected by the fault (since these planes are merely sheared at the 
fault). Electrons diffracted by (1010) pyramidal planes would be affected 
by the fault but would generally be travelling almost parallel to the faults 
and hence would not encounter them so often. It would be necessary 
to consider diffraction from higher order inclined planes such as (1011) 
for which the associated Fourier component of the lattice potential may 
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be very small. It does not seem impossible that further reduction im 
scattering by as much as a factor of 5 or so could occur purely as a result 
of the different geometry in the h.c.p. system. 

(3) The estimate of a stacking-fault density of 107 cm™ (equivalent 
to a fault on every third plane) with the faults scattering independently 
seems dangerously large. Even in highly deformed sodium, where the 
f.c.c. and h.c.p. phases appear to be relatively closer in energy than in 
lithium, Barrett (1956) estimated from x-ray studies that only 1 plane in 5: 
was faulted while in unworked samples the faulting was less. A faulting 
density involving | plane in 10 might therefore be a more reasonable figure 
to assume for the case in question. 

The suggestion (Howie 1960) that faults more than one layer apart 
should scatter independently was made for copper and it is likely that 
faults would interact over greater distances in more free-electron-like 
metals. The assumption of genuinely independent scattering from arrays 
of parallel faults should moreover evidently be applied with the greatest 
caution since any approach to a regular arrangement of faults must lead 
to cooperative effects. The h.c.p. structure could for instance be regarded. 
as f.c.c. with a fault on every second plane and a very large resistivity 
would be expected for perfect crystals of this structure if the faults were 
imagined to scatter independently. 

It seems possible therefore that a figure for the resistivity associated with 
unit density of stacking fault in h.c.p. lithium as low as 4 x 10- ohm-em 
could be understood without any modification to our picture of the way 
in which stacking faults scatter electrons and that the remaining 
discrepancy between this figure and the one quoted earlier may result 
from an over-estimation of the effective stacking fault density. 
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A Systematic Error in the Determination of Dislocation 
Densities in Thin Films 


By R. K. Ham and N. G. SuarPE 


Division of Tribophysics, Commonwealth Scientific and Industrial Research 
Organization, University of Melbourne, Victoria, Australia 


[Received July 17, 1961] 


THE methods of estimating the density of dislocations in thin films des- 
eribed by Bailey and Hirsch (1960) and Ham (1961) both depend on the 
assumption that the dislocation segments are randomly orientated with 
respect to the plane of the film. Bailey and Hirsch have suggested that 
when a bulk specimen is thinned dislocations may shorten their lengths by 
re-orientating themselves so as to lie more nearly normal to the plane of the 
film. A. K. Head (private communication) has pointed out that this 
departure from randomness in the distribution of dislocation orientations 
might be detected by counting the number of intersections which disloca- 
tions make with the surfaces of the film. — If the dislocations are randomly 
orientated with respect to the plane of the film, then the dislocation 
density (i.e. length of line per unit volume) obtained is p’ = 2N’/A, where N’ 
is the number of intersections with both surfaces whose total area is A 
(Frank 1957). Shortening of the dislocations in the way suggested by 
Bailey and Hirsch will cause p’ to be greater than the density p obtained by 
the random line method described by Ham, which counts intersections of 
dislocations with planes normal to the film. 

The dislocation density has been determined both by surface intersection 
counts and by the random line method for the same areas on 20 micrographs 
taken from thin films of a specimen of aluminium cold-rolled to give about 
9° reduction in thickness. The mean density obtained with random lines 
was p=1-78 x 10®°cm/cm?, while that obtained with surface counts was 
p’=2-18x 10®cm/cm?. Thus the mean difference between the densities 
determined by the two methods was 4x 10%cm/cm?. The standard 
deviation of the mean difference was 1-75 x 108 em/cm', so that the ditference | 
between the two results is significant (at the 3-5°% level) and p’ is greater 
than p. Hence the distribution of dislocation orientations is biased in 
favour of the direction normal to the film, as suggested by Bailey and 
Hirsch. 

This means that the true dislocation density in the film (in the absence 
of other systematic errors, e.g. invisibility of certain dislocations) must lie 
between p and p’. Since the difference between p and p’ is about 20% 
then the error in either p or p’ is about 10%. 

It may also be noted that the dislocation re-orientations which can 
oceur during the preparation of a thin film will not affect the number of 
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surface intersections. Thus the result of the surface count p’ gives an 
estimate of the true dislocation density in the bulk (again in the absence of 
other systematic errors, e.g. loss of dislocations during thinning from the 
bulk) if the dislocations are randomly orientated in the bulk, or if the 
orientations of the surfaces examined are selected at random. Hence the 
random line result underestimates the bulk density by about 20%. 

Surface counting has a further advantage over the other methods in 
that it does not depend on a determination of the thickness of the film. 
Unfortunately it becomes very difficult to count surface intersections for 
dislocation densities higher than about 3 x 10° cm/cm?. 
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REVIEWS OF BOOKS 


Elements of Maser Theory. By ArtHurR A. VUYLSTEKE. (Van Nostrand, 1961.) 
[Pp. 363.] Price £3 11s. 6d. 


‘Tuts is a remarkable book, and both the author and his employers, General 
Motors Corporation, are to be congratulated on their foresight in anticipating 
‘that there would be a need for a comprehensive account of the general back- 
ground of theory needed to understand masers. The author (and your reviewer) 
believes that a new branch of technology, quantum electronics, is about to 
develop, and that unlike semi-conductors, for example, classical analogues are 
unlikely to appear. Thus those who are going to be responsible for designing 
the new devices will be faced with the need to have a really good understanding 
of quantum mechanics and statistical mechanics. The attempt is made, in this 
book, to satisfy this need and, in consequence, many pages are devoted to what 
might be called the standard book-work of modern theoretical physics. The 
readers’ requirements are always kept in mind and although I feel that the aver- 
age British physics graduate will be hard pushed to absorb all the material pre- 
sented, he really hasn’t any option if he proposed to contribute to the 
‘development of masers and quantum electronics. Further, once he has 
mastered the subject, he will have achieved something worthwhile. I would 
‘certainly recommend this book to anyone who is prepared to work hard to 
understand masers. Ke Were s: 


‘Computing Methods and the Phase Problem in X-ray Crystal Analysis. By R. 
Peprysky, J. M. RoBertson and J.C. SPEAKMAN. (Pergamon Press, 1961.) 
[Pp. 326.] Price £3 3s. Od. 

A CONFERENCE of crystallographers was held in Glasgow in August of last year 

to discuss computing methods and methods of solving the ‘ phase problem ’—in 

other words, how to reconstitute the image of an object when the phases of the 
waves diffracted by it are unknown. The present book is a summary of the 
proceedings. 

Usually such books are not very successful: they tend to record work at an 
inappropriate stage of development; the papers often do not hang together ; 
-and the discussions, illuminating as they may have seemed to the participants 
at the time, may look rather trivial in the cold light of print. This book, how- 
-ever, emerges as better than the average. It is extremely useful to know what 
various people are doing in the field of computation, for their methods are not 
usually published in the ordinary journals. ; 

Of the 28 chapters in the book, nine are devoted to the phase problem. These 

-are chiefly concerned with so-called ‘ direct methods ’—that is, finding the 

phases by utilizing relations that are likely to exist between the various orders 

of diffracted beams. Other methods discussed are the ‘ heavy-atom ’ method, 
in which the phases are taken to be those of a predominant atom in the crystal; 
the isomorphous-replacement method, in which comparison is made between 

‘two isomorphous crystals; and anomalous dispersion, in which use is made of the 

change in phase produced when the radiation has a wavelength near to an 

-absorption edge of one of the atoms. Nothing at all is included about recent 

uses of Fourier transforms. 

The book will form a useful addition to the library of the practising Bay 
HL: 


‘erystallographer. 
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The Theory of Neutral and Ionized Gases. Edited by C. Dx Wirt and J. F 
Detorur. (New York: John Wiley & Sons, Inc., 1960.) [Pp. 469.] Price 
£7 Os. Od. 


Tuer 1958 and 1959 Les Houches Summer Schools in Theoretical Physics were 
devoted to the theory of many-body systems. The 1958 Summer School 
considered the quantum mechanics of such systems and its proceedings have by 
now become a recognized vade-mecum. The 1959 Summer School was devoted 
to the discussion of the statistical mechanics of systems of interacting particles 
and the volume reviewed here constitutes its proceedings. Unfortunately these 
proceedings have taken about half a year longer to be published than the 1958 
ones, but now that they have appeared, they should also soon become part of the 
all too scarce number of compendia on the theory of many-body systems. 

Once again the directors of the Les Houches Summer Schools succeeded in 
attracting a number of theorists who themselves contributed greatly to the 
subject matter under discussion. The longest contribution is the first one by 
Montroll in which he discusses ensemble theory, cluster integral theory, trans- 
port processes, random walk, phase transitions, and the Ising problem. This is 
followed by Van Hove’s lectures on non-equilibrium processes, and a series of 
lectures by Deleroix, Kruskal, Kaufman, and Denisse on various aspects of 
plasma theory. The volume is completed by a discussion by Schatzman of 
plasma-physics as applied to astrophysics. 

Unfortunately the editors have not supplied an index which would have 
appreciably enhanced the already considerable value of this book. Surely, if 
these lectures are worth publishing—which they certainly are—they are worth 
publishing with an index. D. TER Haar. 


Balthasar van der Pol: Selected Scientific Papers. (2 volumes.) Edited by 
H. BremMMER and C.J. Bouwkamp. (Amsterdam: North Holland Publishing 
Company, 1960.) [Pp. xvi+1339.] Price £6 13s. Od. 


BALTHASAR VAN DER Por’s scientific papers range over a wide field in the border 
region between Physics, Electrical Engineering and Mathematics. They 
include papers of fundamental importance to which workers in the appropriate 
fields need to make frequent reference and there must be many who will want to 
have these volumes on their shelves. 

The topics covered may be divided broadly into three groups. First, there 
are the papers on the propagation of radio waves including the famous series on 
propagation over a spherical earth with finite conductivity. Second, there are 
the papers on thermionic valves including studies of relaxation oscillations and 
the two well known papers written in collaboration with Sir Edward Appleton. 
Third, there are the mathematical papers which cover such topics as operational 
calculus, theory of numbers, Tchebycheff polynomials and Jacobi functions. 

These volumes are very well reproduced from the original journals by a 
photographic process. About 71% of the material is in English, 13° in Dutch 
and the remander about equally in French and German. Van der Pol’s doctoral 
thesis, in Dutch, is included, but the editors point out that most of the material 
is repeated elsewhere in English. Two plates show photographic reproductions. 
from van der Pol’s manuscript notes. KRiGuB: 


Relativity. The General Theory. By J. L. Syner. (North-Holland Publishing: 
Company, 1960.) [Pp. 505.] Price £5 10s. Od. 


THIs is the companion volume to Professor Synge’s Relativity: the special theory 


(1956) but it can be read independently. Like the earlier book, it is a major: 
addition to the literature. 
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In a general way, the ground covered is about what we should expect of a 
comprehensive treatise. However, the approach and development and the 
selection of particular topics within the general scope are all highly individual- 
istic. Synge exploits more recondite resources of differential geometry than those 
that are at the command of most other writers. He combines the physics and 
the geometry almost axiomatically and without appeal to the outworn traditional 
heuristic arguments. He gives a many-sided analysis of the field equations and 
their consequences. He considers many applications of the general theory but 
mainly for the sake of the principles or methods they illustrate rather than the 
immediate particular results. While he expressly does not aim at completeness 
in his treatment of the ground covered, he provides an extensive bibliography 
(64 pages) that includes references to numerous topics not dealt with in the text. 

The work is the outcome of Synge’s unrivalled experience of relativity in 
research, teaching and collaboration with younger colleagues. As we have 
come to expect, the presentation is spiced by Synge’s mordant wit. His book 
will be of immense value to all workers in relativity and it will be a special boon 
to those beginning to do research in the field. The only readers on whom its 
merits would probable by wasted are those who are quite without previous 
knowledge of the theory. W. H. McC. 


Resonance Absorption in Nuclear Reactors. By LAWRENCE DRESNER. (Perga- 
mon Press, 1960.) [Pp. 131.] Price £2 Os. 0d. 


Tuts book is extremely useful, containing a comprehensive summary of the 
theory and a comparison of theory with experiment. It is divided into eight 
chapters, commencing with a short historical review and two chapters on the 
theory of resonance absorption in an homogeneous medium. The interesting 
part of the book, from the point of view of people interested in reactor design 
physics, is contained in the subsequent chapters on heterogeneous resonance 
integrals. These chapters contain a description of the approximate analytic 
theory for estimating the resonance absorption in a lattice structure. It is 
shown how, with certain assumptions, the calculation in the heterogeneous 
case can be reduced to that of the homogeneous case. The argument is more 
satisfying from a logical point of point, than the theory previously used to 
justify the form of the empirical relations fitted to the experimental results. 

The use of the purely numerical Monte Carlo method is not discussed in detail, 
although some results are quoted. It would have been interesting to compare 
the Monte Carlo results with the approximate analytic ones for a greater variety 
of cases. The comparison with theory and experiment illustrate the usefulness 
of Hellstrand’s series of experiments. 

This book should be welcomed, it will be very useful to reactor design 
physicists. Jats 


Radio Waves in the Ionosphere. By K.G. BuppEN. (Cambridge: University 
Press, 1961.) [Pp. 542.] Price £4 15s. Od. 


THE propagation of radio waves through the inhomogeneous, anisotropic 
ionosphere presents problems of considerable complication, especially when the 
electron concentration changes appreciably within the space of one wavelength. 
The results of work on these problem are spread over a wide range of papers so 
that it is not easy to obtain a clear picture of what has been accomplished. 
Moreover, the mathematical apparatus involved is complicated, so that many 
physicists and engineers find difficulty in understanding clearly how the prob- 
lems have been approached and. what the results are. by 

Dr. Budden’s book provides, in a very clear form, an exposition of these 
matters and a clear statement of what has been accomplished. His ability as a 
teacher is evident from the way in which he introduces each part of the subject 
from first principles, and explains any less familiar mathematical concepts. 
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The book can be unreservedly recommended to all who are concerned with the 
theory of radio waves in the ionosphere. For those familiar with the subject it 
will suffice to list some of the section headings to indicate its value: 

“ Slowly varying media”, ““ W.K.B. solutions”, “ Ray theory for oblique 


incidence ”’, ‘‘ The Booker quartic ”’, ‘‘ The ray velocity and the ray surface ”’, 
“The general problem of ray tracing”, “ Coupled wave equations’, “ The 
Z-trace ’’, ‘“ The coupling echo ”’, “ Coupling between upgoing and downgoing 
waves ”’, ‘‘ The phase integral method ”’, “‘ Reciprocity.” 


It will also form a valuable text-book for anyone setting out to study some of 
the more complicated aspects of electro-magnetic wave theory, and for them 
the examples included at the ends of the separate chapters should prove useful. 

The dedication of the book “‘ To the memory of Douglas R. Hartree ” is a 
happy and most appropriate one. J.A.R. 


Electrical Noise. Fundamentals and Physical Mechanism. By D. A. BELL. 
(D. van Nostrand, 1960.) [Pp. 342.] Price £2 10s. Od. 


Turs is an excellent book which will undoubtedly be warmly welcomed by 
engineers, research workers and teachers in a wide variety of fields. The 
author’s intention was to give a complete and coherent account of the present 
state of knowledge and he has amply succeeded in this scholarly, closely 
reasoned and yet readable volume. 

In the first few chapters a very thorough account is given of the fundamentals 
of noise theory, namely equipartition, statistical mechanics and Nyquist’s 
theorem. Succeeding chapters deal with noise in various types of valve, photo- 
multipliers and more recent devices including masers and parametric amplifiers. 
This is followed by a detailed discussion of current noise in semi-conductors, noise 
in thin films and Barkhausen noise. Some mention is also made of radiation 
detectors and measuring techniques although space forbids other than an out- 
line of the principles involved. Throughout, the text is liberally dressed with 
references to original papers which provide a full, if not exhaustive, bibliography. 

It is difficult to criticize this book without niggling at a few details and it 
deserves far better treatment than this. A. H. 


Quantum Theory of Atomic Structure, Vol. 2. By J.C. Stater. (McGraw-Hill 
Publishing Company Ltd., 1960.) [Pp. 4389.] Price £5 ls. Od. 


Tuts book is one of the best attempts to lead the advanced student easily into 
the sophisticated theory of the quantum structure of atoms. At the same time 
it provides a very useful collection of reference data and applications of 
advanced techniques. It might be regarded as for the specialist and for reference 
but it is not the compendium of theorems that this might seem to imply. There 
is an air of teaching the reader to do it for himself. Professor Slater writes as 
if he thinks you can do quite a lot of advanced theory on an odd piece of paper 
and this book will show you how. The content, however, is really quite heavy 
and, for example, one finds a very informative review of Racah’s methods for 
complex spectra but without all the details of his papers. 

The main contents are the mathematical theory and an account of quantita- 
tive predictions of atomic multiplet levels, together with many reference tables, 
the theoreties of magnetic properties and hyperfine structure and the treatment 
of the fractional parentage theory of an configurations by projection operators. 
Finally there is a complete bibliography of all theoretical investigations on 
atomic structure. 

It will be a most valuable book for those who research or lecture on atomic 
structure and it should be in all science libraries. coe ia 3S 
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Introduction to Mechanics, Matter, and Waves. By U. Inaarp and W. L. 
KravusHaar. (Addison-Wesley Publishing Company Inc., 1960.) [Pp. 672.] 
Price £2 13s. Od. 

Tuts textbook is based upon a year’s course in general physics at MIT, and it 

covers mechanics, dynamics, wave motion, elementary properties of matter, 

the kinetic theory of gasses, and a little thermodynamics, starting essentially 
from scratch. As is common with such American texts it assumes a concurrent 
course in mathematics, so that the mathematical standard of the book becomes 
more advanced as it progresses; calculus is used freely throughout the latter 
half. The basic concepts, such as force, energy and so on, are discussed very 
carefully, with the aid of many excellent diagrams illustrating actual or con- 
ceptual experiments, and the argument is developed from these as far as such 
relatively sophisticated topics as Rutherford scattering, coriolis force, gyroscopic 
motion, the Maxwell—Boltzman distribution, entropy and probability, and the 

Bernoulli effect. There are a number of problems after each chapter. In 

general the presentation is good, and the book seems likely to be a useful one. 

de Hie Be 


The Theory of Transition-Metal Jon. By J. S. Grirriru. (Cambridge: 
University Press, 1961.) [Pp. x+453.] Price £4 15s. Od. 


TRANSITION metal compounds are interesting because of the great variety of 
their properties, chemical, optical and magnetic. Basically these all derive 
from an incomplete electron shell. There are many energy levels of an ion 
within a relatively small energy range. The electron charge distribution is 
therefore flexible and can adapt itself to its environment in such a way as to 
lower the energy, thereby making a particularly stable complex. There are 
many levels separated from the ground state by energy differences in the range 
afforded by quanta of visible light. Light is therefore selectively absorbed, 
giving transition metal compounds their colour, and on a more sophisticated 
level, their interesting and complex absorption spectra. The magnetic moments 
of individual electrons in the ion are not forced by the exclusion principle to cancel 
each other out, and the ion therefore has a magnetic moment. 

Since the pioneer work of Van Vleck in the 1930s, more and more attention 
has been devoted to those interesting phenomena and to the theories which 
unify them. Professor Griffith’s book is both timely and admirable. It 
brings together the various physical and chemical properties of transition 
metal ions in a comprehensive manner and relates them to a well-developed 
theoretical framework. The exposition is admirably lucid. General principles 
are first outlined and then specified problems are worked out in fair detail. 
Nowhere is the reader bogged down by a surfeit of generality, though by the end 
of an argument he may have been led to a very general point of view. Professor. 
Griffith prefers to treat a problem first approximately, and later to point out 
what modifications are required to incorporate the features which were neglected. 

Five chapters are devoted to theoretical ground work. This covers ground 
which is to be found elsewhere, and makes the book rather bulky. However, 
this approach has the merit of making the book self-contained and giving it 
unity. Some of the results presented, particularly in the chapter on groups, are 
new. The general theory is always economically set out for the purpose in. 
hand. The remainder of the book is devoted to a more detailed discussion of 
various aspects of transition metal ions as follows: complex ions; crystal-field 
theory and the weak-field coupling scheme; the strong-field coupling scheme; 
paramagnetic susceptibilities; optical spectra and thermodynamic properties; 
paramagnetic resonance. MPH. LP. 
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An Introduction to Applied Anisotropic Elasticity. By R. F. 8. Hearmon. 
(Clarendon Press; Oxford University Press, 1961.) [Pp. 136.] Price 
£1 15s. Od: 


Arrer assuming that the reader has a knowledge of the fundamentals of 
theoretical elasticity, the first part of this book gives an authoritative discussion 
of the stress-strain relations in aelotropic media. The effects of symmetry are 
described in full, and there is a chapter on single crystals, sheet materials and 
polycrystalline aggregates. 

In later chapters, the author has collected, from many sources, solutions to a 
large number of problems in anisotropic elasticity. In a few cases solutions are 
given out of context and are, therefore, hard to follow. However, the com- 
prehensive list of references makes it easy to refer to the original work. Ina 
chapter on wave propagation, recent research work is described, but one gains 
the impression that the presentation could be simplified if a more general 
approach were adopted. 

The book fills a gap and it will be found very useful, not only by the engineer 
or physicist, but also by the mathematician working in this field. Week bs: 


Dispersion Relations. Edited by G. R. Screaton. (Oliver & Boyd, 1961.) 
[Pp. 290.] Price £2 15s. Od. 


THis volume contains the lecture notes of the Scottish Universities Summer 
School on Dispersion Relations arranged by N. Kemmer in August 1960. 

The first chapter is by J. D. Jackson and gives an introduction to dispersion 
relation techniques that is suitable for a reader having only an elementary 
knowledge of quantum field theory. In the second chapter J. C. Polkinghorne 
gives an account of the basic techniques used to establish analyticity properties 
in perturbation theory. This is balanced by a discussion of analyticity properties 
of scattering amplitudes based on axiomatic field theory given by W. Thirring. 
The mathematical foundations of quantum field theory are discussed in a short 
chapter by J. M. Jauch. 

The four remaining chapters give an excellent introduction to the problems 
of applying dispersion theory to physical processes. M. J. Moravesik discusses 
the use of the ‘ nearest singularity ’ approximation with particular reference to 
nucleon-nucleon scattering. G. F. Chew gives a very readable account of the 
ambitious programme directed towards basing a theory of strong interactions 
on double dispersion relations and unitarity; this includes a discussion of the 
Mandelstam representation and partial wave dispersion relations. W. R. 
Frazer gives an account of the more limited problem of the nucleon form-factor, 
based on the two-pion approximation. 8. Fubini describes the use of single 
variable dispersion relations in pion physics. 

This book gives an excellent introduction to the important subject of analytic 
structure and elementary particle physics. It is suitable for reading at graduate 
level following an elementary course on quantum field theory. R. J. E. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 


JACEE I Phil. Mag. Ser. 8, Vol. 6, Pl. 130. 
Fig. 2 


A large crystal (top) showing changes of hue as it grows ACKOSS & SETICS of shallow 
steps of height 300 to 8004. The lower crystal, thickening by layers 
spreading over its surface, is becoming polygonal. 


i A ‘tal surfaces after initia- 
f thickness 300 to 12004 spreading over crys ti 
Be ei at the edge. Temperature —16-5°c. Excess vapour density 


0-4 x 10-* g em-3. 


J. HALLETT Phil. Mago Ser. 8;-Voloonhlel3 ie 


Fig. 8 


Layers of thickness 12004 growing around the crystal edges, after initiation by 
contact with thicker crystals. 


Fig. 9 


An irregularly shaped layer growing over an ice surface. 


ARALLET Pala Magy ser. 8, Vol.16,Plal'32: 


Fig. 10 


A dendritic layer with growth maxima in <1120) directions. Note that all 
crystals are growing on steps in the substrate. 


Fig. 15 


Dendritic growth of ice crystals on a freshly cleaved mica surface. 


EILEEN M. WILKS Phil. Mag. Ser. 8, Vol. 6, Pl. 133. 


Photomicrograph (x 140) of a natural octahedron surface showing a low-angle 
pyramidal trigon X and a steep-sided flat-based trigon Y. 


Fizeau fringes across a steep-sided flat-based trigon ( x 400). 


EILEEN M. WILKS Phil. Mag. Ser. 8, Vol. 6, Pl. 134. 


Fig. 3 


Fizeau fringes across a low-angle pyramidal trigon (x 400). 


Fig. 5 


Photomicrograph (x 28) showing low-angle pyramidal trigons and steep-sided 
trigons. 
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(a) Initial precipitation in Al4+4% Cu alloy aged for a few minutes at ~270°o, 
The majority of the precipitates are 6’ and are believed to be at the foil 
surfaces. (b) Same area at ~340°c. Note the coarsening. Large 
precipitates A, B and C are 6. Foil surface near (112). x30 000. 
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G. THOMAS and M. J. WHELAN Phil. Mag. Ser. 8, Vol. 6, PI. 136. 
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(a) Precipitation at ~270°C characteristic of a foil with near (100) surface. 
The large 6’ precipitates viewed edge-on are known to have nucleated in the 
interior on dislocations. The irregular precipitates of lower visibility 
are 6’, and are thought to lie at the (100) surfaces. (b) Same grain at 
~330°c showing extended precipitate patterns of 9 phase. The extended 
patterns are believed to be at the surfaces. 30000. 


G. THOMAS and M. J. WHELAN Phil. Mag. Ser. 8. Vol. 6, eal, We 


Region with near (110) surface showing 6 and @’ precipitates at ~320°c. 
x 30 000. 


Fig. 4 
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(110) area at ~320°o, showing @ and @’ precipitates. x 30000. 
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Examples of precipitates at ~300°C near a triple grain boundary junction. 
A, inclined grain boundary with preferential precipitation at the surfaces 
B and C. Note the regions denuded of precipitates near the grain 
boundaries; (a) refers to the second precipitation, (5) refers to the third 


precipitation. x 15 000. 


G. THOMAS and M. J. WHELAN Phil. Mag. Ser. 8, Vol. 6, Pl. 139. 


Ip (1) 0-5 (d) 
| een SS 


(4) @ precipitates in a thin region. The precipitates run right through the foil. 
Temperature ~430°c. x 10000. —(b) 6 precipitate running right through 
the foil. Note the inclined boundary showing thickness interference fringes. 


x 30 000. 
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(a) Region showing initial precipitates after 30 min. at ~250°c. The precipi- 
tates are probably 6’. x 20000. (b) Same region after a small tempera- 
ture rise, showing growth of large precipitates surrounded by a denuded 
zone. ‘The precipitates A and B are probably @. x 10 000. 


G. THOMAS and M. J. WHELAN Phil. Mag. Ser. 8, Vol. 6, PI. 140. 


Fig. 8 Fig. 9 


Thin foil aged for 10 min. at 270°c in high Thick region of the same speci- 
vacuum prior to examination, Compare men as fig. 8 after the same 
with fig. 1. x 30 000. ageing treatment and a 


further electropolish. The 0’ 
precipitate patterns are 
probably the result of 
nucleation at helical dis- 
locations. x 26 000. 


Fig. 10 


Examples of precipitation at 
helical dislocations. Aged 
in bulk for 2 min at 270°c. 
Note the precipitates A and 
the dislocation B. x 26000. 
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G. THOMAS and M. J. WHELAN Phil. Mag. Ser. 8, Vol. 6, PI. 141. 
Fig. 11 


Selected area electron diffraction 
pattern from a thin foil aged at 
~270°c, showing  cross-grating 
pattern from 6’ marked with 
arrows. A, 200 matrix and 6’ 
reflections. Foil surface near (011). 


Fig. 12 


(a) (b) 


(a) 0’ reflections of type 110 (marked with arrows) from precipitates in a foil of 
(001) orientation similar to fig. 2 (a). Ageing temperature ~250°c 
(b) Diffraction pattern from the same grain at ~330°c, showing 6 cross- 
grating pattern. a-axes of @ phase are outlined. Note the 200 type 
matrix reflection A, and the 220 type 6 reflection B. 


G. THOMAS and M. J. WHELAN Phil. Mag. Ser. 8, Vol. 6, Pl. 142. 


0-5 
ae 
Sequence of micrographs showing dissolution of # precipitates as the temperature 
is raised between 400 and 500°c. Note the rounding-off of the angular- 
shaped precipitates. Foil surface is (100). x 30 000. 


J Gy By RNE@eteai: Phil. Mag. Ser. 8, Vol. 6, PI. 143. 


Vig. 1 


Apparatus for tensile testing at 4:2°xK. A, vacuum jacketed helium dewar ; 
B,, compression tube ; By, flange for attaching helium dewar; Bs, 
attachment for helium transfer tube ; B,, exhaust opening ; B,, flange 
for attaching to Instron ; C, tension tube; D, and D,, fixtures for 
holding specimen grips; E, and E,, grips for holding specimen ; 
F, adaptor for attaching tension tube to Instrom ; G, neoprene gasket ; 
H, nitrogen dewar ; J, helium gas inlet for purging ; C fits into opening 
in J; K, typical specimen ; L, bubbler for helium gas outlet. 


E. GILLAM and M. PHILLIPS Phil. Mag. Ser. 8, Vol. 6, Pl. 144. 
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Single crystal gold film bombarded with 300 ev argon ions. x40 000. 


As fig. 1, selected area 5 microns in diameter. 


E. GILLAM and M. PHILLIPS Phil. Mag. Ser. 8, Vol. 6, Pl. 145. 
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Diffraction pattern at edge of bombarded specimen. 5 microns. 


Diffraction pattern at edge of unbombarded specimen. 5 microns. 


E, GILLAM and M. PHILLIPS Phil. Mag. Ser. 8, Vol. 6, Pl. 146. 


Fig. 5 


Gold film annealed at 350°c after ion bombardment. x 100 000. 


